Index: /issm/trunk-jpl/src/c/shared/Elements/ComputeMungsmTemperaturePrecipitation.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/ComputeMungsmTemperaturePrecipitation.cpp	(revision 18969)
+++ /issm/trunk-jpl/src/c/shared/Elements/ComputeMungsmTemperaturePrecipitation.cpp	(revision 18969)
@@ -0,0 +1,31 @@
+/* file:  ComputeMungsmTemperaturePrecipitation.cpp
+   Compute the temperature and precipitation at time t from 
+   the data at present day and lgm.
+   The interpolation is done from some factors extracted from the MUNGSM
+ */
+
+#include "./elements.h"
+#include "../Numerics/numerics.h"
+
+void ComputeMungsmTemperaturePrecipitation(IssmDouble TdiffTime, IssmDouble PfacTime,
+					   IssmDouble* PrecipitationsLgm, IssmDouble* PrecipitationsPresentday,
+					   IssmDouble* TemperaturesLgm, IssmDouble* TemperaturesPresentday,
+					   IssmDouble* monthlytemperaturesout, IssmDouble* monthlyprecout){ 
+
+  IssmDouble monthlytemperaturestmp[12],monthlyprectmp[12];
+  IssmDouble tdiffh;  
+
+
+  for (int imonth = 0; imonth<12; imonth++){
+    tdiffh = TdiffTime*( TemperaturesLgm[imonth] - TemperaturesPresentday[imonth] );
+    monthlytemperaturestmp[imonth] = tdiffh + TemperaturesPresentday[imonth] ;
+
+    monthlyprectmp[imonth] =min(1.5, PrecipitationsPresentday[imonth] * pow(PrecipitationsLgm[imonth],PfacTime));   // [m/month]
+
+    /*Assign output pointer*/
+    *(monthlytemperaturesout+imonth) = monthlytemperaturestmp[imonth];
+    *(monthlyprecout+imonth) = monthlyprectmp[imonth];
+  }
+  // printf(" tempera %f\n",monthlytemperaturestmp[1]);
+}
+ 
Index: /issm/trunk-jpl/src/m/enum/SurfaceforcingsIsmungsmEnum.m
===================================================================
--- /issm/trunk-jpl/src/m/enum/SurfaceforcingsIsmungsmEnum.m	(revision 18969)
+++ /issm/trunk-jpl/src/m/enum/SurfaceforcingsIsmungsmEnum.m	(revision 18969)
@@ -0,0 +1,11 @@
+function macro=SurfaceforcingsIsmungsmEnum()
+%SURFACEFORCINGSISMUNGSMENUM - Enum of SurfaceforcingsIsmungsm
+%
+%   WARNING: DO NOT MODIFY THIS FILE
+%            this file has been automatically generated by src/c/shared/Enum/Synchronize.sh
+%            Please read src/c/shared/Enum/README for more information
+%
+%   Usage:
+%      macro=SurfaceforcingsIsmungsmEnum()
+
+macro=StringToEnum('SurfaceforcingsIsmungsm');
Index: /issm/trunk-jpl/src/m/enum/SurfaceforcingsPrecipitationsLgmEnum.m
===================================================================
--- /issm/trunk-jpl/src/m/enum/SurfaceforcingsPrecipitationsLgmEnum.m	(revision 18969)
+++ /issm/trunk-jpl/src/m/enum/SurfaceforcingsPrecipitationsLgmEnum.m	(revision 18969)
@@ -0,0 +1,11 @@
+function macro=SurfaceforcingsPrecipitationsLgmEnum()
+%SURFACEFORCINGSPRECIPITATIONSLGMENUM - Enum of SurfaceforcingsPrecipitationsLgm
+%
+%   WARNING: DO NOT MODIFY THIS FILE
+%            this file has been automatically generated by src/c/shared/Enum/Synchronize.sh
+%            Please read src/c/shared/Enum/README for more information
+%
+%   Usage:
+%      macro=SurfaceforcingsPrecipitationsLgmEnum()
+
+macro=StringToEnum('SurfaceforcingsPrecipitationsLgm');