not being able to store "PRECCFRQ"

[sandeep_sahany@indiatimes_com]

Dear Cam-users,
I am trying to store the following variable but not being able to do so.

VARIABLE = PRECCFRQ, PRECLFRQ, PRECCINT, PRECLINT

APPROACH USED = As the above mentioned variable are in the master field list, I added them in namelist as follows
FINCL1 = " PRECCFRQ" , "PRECLFRQ" , "PRECCINT" , "PRECLINT"

Problem facing = The model is running successfully. But it stores the above all variable as zero.

Then I tried to look into the code and seen that, the following are the only locations where, let say " PRECCFRQ" is involved :

models/atm/cam/src/control/comctl.h: real(r8) precc_thresh ! Precipitation threshold for PRECCINT and PRECCFRQ
models/atm/cam/src/control/history.F90: call addfld ('PRECCFRQ','fraction',1, 'A',&
models/atm/cam/src/control/history.F90: 'mm/hr is set to zero -- to get intensity divide by PRECCFRQ)',phys_decomp)
models/atm/cam/src/control/runtime_opts.F90:! precc_thresh Precipitation threshold to use for PRECCINT and PRECCFRQ (mm/hr)
models/atm/cam/src/control/runtime_opts.F90:! Precipitation threshold for PRECCINT, PRECLINT, PRECCFRQ, and PRECLFRQ output fields


From the above it seems that it is not calling outfld, which might be the reason.
If any body can point me the location where "PRECCFRQ" is calculated, which I am not being able to find out, will be a great help.

Thanks and regards
Saroj

 

[rneale]

You have it correct. Every field you include in the fincl1 namelist statement needs an addfld call and outfld call to be included in the history file.

 

[sandeep_sahany@indiatimes_com]

Dear Neale,
Thanks for your constant help towards solving the problems encountered in my research work. Yes, I am aware that outfld needs to be called to output a variable in the history file. But, my problem is to find out the exact location where the outfld should be called, for let's say "PRECCFRQ" .

Warm regards,
Saroj

 

[rneale]

There is currently no PREECFRW outfld call as far as I can tell, so the documentation is not quite right for this variable. But it's pretty easy to include this yourself.

I would have an outfld call for PRECCFRQ at the same point as the outfld call for PRECC which I think is in tphysbc.F90. You just need to output a variable that is = 1 if PRECC is above the minimum threshold (could be 0.1 mm/hr as in the documentaion, but it could be anything you want it to be) or 0 if not. Then the history file mechanism would calculate the appropriate average of this (I think you wanted monthly).

 

[sandeep_sahany@indiatimes_com]

Dear Neale,
Thanks for your guidance. Yes, now I am being able to store it.

Regards,
Saroj