ComputeMungsmTemperaturePrecipitation.cpp File Reference

#include "./elements.h"
#include "../Numerics/numerics.h"
#include <cmath>

Go to the source code of this file.

Functions
void	ComputeMungsmTemperaturePrecipitation (IssmDouble TdiffTime, IssmDouble PfacTime, IssmDouble *PrecipitationsLgm, IssmDouble *PrecipitationsPresentday, IssmDouble *TemperaturesLgm, IssmDouble *TemperaturesPresentday, IssmDouble *monthlytemperaturesout, IssmDouble *monthlyprecout)

Function Documentation

ComputeMungsmTemperaturePrecipitation()

void ComputeMungsmTemperaturePrecipitation	(	IssmDouble	TdiffTime,
		IssmDouble	PfacTime,
		IssmDouble *	PrecipitationsLgm,
		IssmDouble *	PrecipitationsPresentday,
		IssmDouble *	TemperaturesLgm,
		IssmDouble *	TemperaturesPresentday,
		IssmDouble *	monthlytemperaturesout,
		IssmDouble *	monthlyprecout
	)

Definition at line 11 of file ComputeMungsmTemperaturePrecipitation.cpp.

                                                                                 { 
 
  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/yr]
 
    /*Assign output pointer*/
    *(monthlytemperaturesout+imonth) = monthlytemperaturestmp[imonth];
    *(monthlyprecout+imonth) = monthlyprectmp[imonth];
  }
  // printf(" tempera %f\n",monthlytemperaturestmp[1]);
}