problems in output of FUS,FDS,FUL,FDL

[boyle5@llnl_gov]

The above variables appear in the master field list and so if cited in the namelist they will appear in the output file without any error.
But all the values are zero - the reason is that the appropriate outfld call is only included for the SCAM compilation.
Thus, if I want the actual values of these variables I must modify the code in radlw.F90 and radsw.F90. I deal with a number of collaborators and thus many source codes for the CAM, it is very difficult to always remember to make the changes necessary to get the radiation fluxes.

In addition, to get UBOT, VBOT, QBOT, ZBOT one must also change the source code. These variables are only valid if COUP_CSM is defined.

I have two problems (1) that the code is allowed to run when putting out meaningless output, that is only discovered after the run - (2) that to get variables on the master list the source code needs to be modified.

I would like the above variable to be available without source code modification.

 

[rneale]

I don't think any of those variables are actually on the 'default' master field list so you have to explictly add them to the namelist. CAM is designed to run under a number of non-default configurations (coupled, SCAM) and it is certainly not difficult to break the model. In a vanilla/default run only the default diagnostic list is really guaranteed. Just as a note the *BOT variables are just the lowest level values of the 3D equivalent e.g. UBOT is the lowest level of U.

 

[boyle5@llnl_gov]

Richard,

Thanks for the reply - I guess I was just venting frustration. A lot of cycles wasted just producing fields of zeros.
I would like to advocate that the code be changed such that the radiative fluxes are available without changing the source code.

I understand that UBOT etc are the lowest level of winds. For looking at surface processes at high frequency, it is really nice to just use 1/26 of the storage for the lowest level winds, temp and moisture.
In the old days, we use to get 2 meter T and Q and ten meter winds, now the model only seems to provide the 2 meter T.

--Jim