ComputeD18OTemperaturePrecipitationFromPD.cpp File Reference

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

Go to the source code of this file.

Functions
void	ComputeD18OTemperaturePrecipitationFromPD (IssmDouble d018, IssmDouble dpermil, bool isTemperatureScaled, bool isPrecipScaled, IssmDouble f, IssmDouble *PrecipitationPresentday, IssmDouble *TemperaturePresentday, IssmDouble *PrecipitationReconstructed, IssmDouble *TemperatureReconstructed, IssmDouble *monthlytemperaturesout, IssmDouble *monthlyprecout)

Function Documentation

ComputeD18OTemperaturePrecipitationFromPD()

void ComputeD18OTemperaturePrecipitationFromPD	(	IssmDouble	d018,
		IssmDouble	dpermil,
		bool	isTemperatureScaled,
		bool	isPrecipScaled,
		IssmDouble	f,
		IssmDouble *	PrecipitationPresentday,
		IssmDouble *	TemperaturePresentday,
		IssmDouble *	PrecipitationReconstructed,
		IssmDouble *	TemperatureReconstructed,
		IssmDouble *	monthlytemperaturesout,
		IssmDouble *	monthlyprecout
	)

Definition at line 10 of file ComputeD18OTemperaturePrecipitationFromPD.cpp.

                                    {
 
  IssmDouble monthlytemperaturestmp[12],monthlyprectmp[12];
  IssmDouble deltaTemp;
 
  /* Constants */
  // dpermil = 2.4;/*degrees C per mil*/
 
  /*Create Delta Temp to be applied to monthly temps and used in precip scaling*/
  deltaTemp = dpermil * (d018+34.83);   
 
  for(int imonth = 0; imonth<12; imonth++){
 
    if(isTemperatureScaled)monthlytemperaturestmp[imonth] = TemperaturePresentday[imonth] + deltaTemp;
    else{
       monthlytemperaturestmp[imonth] = TemperatureReconstructed[imonth];
       deltaTemp=TemperatureReconstructed[imonth]-TemperaturePresentday[imonth];
    }
 
    if (isPrecipScaled)monthlyprectmp[imonth] = PrecipitationPresentday[imonth]*exp((f/dpermil)*deltaTemp);
    else monthlyprectmp[imonth] = PrecipitationReconstructed[imonth];
 
    /*Assign output pointer*/
    *(monthlytemperaturesout+imonth) = monthlytemperaturestmp[imonth];
    *(monthlyprecout+imonth) = monthlyprectmp[imonth];
  }
}