Using the RadApp tool

Overview

The RadApp tool is a lightweight application running thru MAPL_Cap that runs the Radiation gridded component as well as the PCHEM gridded component to provide aerosols to the radiation. Since it is run through MAPL_Cap the application can be run for multiple time steps and output any export of the Radiation components through the History, just as in the full GEOS5 model. The application uses ExtData to drive the tool with external data on disc. The user should consult with the documentation on ExtData for details in how this component works.

Setup

This is the current checkpoint procedure. Since you are running the Radiation code you must have a working copy of the full GEOS5 model. This page will not detail how to do this. The current RadApp has been developed and tested with the Ganymed-3_0 and Ganymed-4_BETA(1-3) tags.

1. Checkout and built the model tag
2. Checkout the application by issuing: cvs co -r tag GEOSdiag_App
3. Make sure the environment variable ESMADIR is setup probably for the model build you want to link the application to. Typically this is /path_to_model/GEOSagcm. If this is correct there should be a src and linux directory under ESMADIR. Make you source $ESMADIR/Linx/bin/g5_modules to make sure the modules you need to build are added.
4. Now simply issue: make install, and RadApp.x should be built if everything is ok.

Resource File

The following are keywords that can be used in the resource file for RadApp. The name of the resource file for RadApp is specified in the CAP.rc file with the ROOT_CF keyword. Many of these keywords are the same as in the AGCM.rc for the full model.

NX: integer layout in first dimension

NY: integer layout in the second dimension

IM: integer gridsize in first dimension

JM: integer gridsize in second dimension

LM: integer number of vertical levels

GRIDNAME: name of the grid, follows the same convention as model; for example a 2 degree lat-lon grid would be PC144x91-DC

SOLAR_DT: integer timestep for full radiation calculation

IRRAD_DT: integer timestep for full radiation calculation

IRRAD_CALLED_LAST: 0 if you want full radiation before flux evaluation, otherwise 1

SOLAR_CALLED_LAST: 0 if you want full radiation before flux evaluation, otherwise 1

RADCOUPLE_METHOD:

You can restart RadApp using a restart for the internal states of SOLAR and IRRAD and checkpoint. Or you can just omit these entirely to bootstrap the internal state. The same conventions and options apply as in the full model.

SOLAR_INTERNAL_RESTART_FILE: solar_internal_rst

SOLAR_INTERNAL_RESTART_TYPE: pnc4

SOLAR_INTERNAL_CHECKPOINT_FILE: solar_internal_checkpoint

SOLAR_INTERNAL_CHECKPOINT_TYPE: pnc4

IRRAD_INTERNAL_RESTART_FILE: irrad_internal_rst

IRRAD_INTERNAL_RESTART_TYPE: pnc4

IRRAD_INTERNAL_CHECKPOINT_FILE: irrad_internal_checkpoint

IRRAD_INTERNAL_CHECKPOINT_TYPE: pnc4

You must also specify the number of bands and the location of the optics files for example:

DU_OPTICS: ExtData/g5chem/x/opticsBands_DU.v11_2.nc

SS_OPTICS: ExtData/g5chem/x/opticsBands_SS.v3_3.nc

SU_OPTICS: ExtData/g5chem/x/opticsBands_SU.v1_3.nc

OC_OPTICS: ExtData/g5chem/x/opticsBands_OC.v1_3.nc

BC_OPTICS: ExtData/g5chem/x/opticsBands_BC.v1_3.nc

NUM_BANDS: 18

You can also specify the RATS and AERO providers. Currently this is either PCHEM or you can Choose RADENV and fill these through ExtData

RATS_PROVIDER: PCHEM or RADENV AERO_PROVIDER: PCHEM or RADENV

If you are running PCHEM for the aerosols you will need to provide the files, for example:

1. For Climatological Aerosols (or non-interactive GOCART)
2. -------------------------------------------------------

AEROCLIM: ExtData/g5chem/L72/aero_clm/dR_MERRA-AA-r1.aer_Nv.%y4%m2clm.nc4 AEROCLIMDEL: ExtData/g5chem/L72/aero_clm/dR_MERRA-AA-r1.del_aer_Nv.%y4%m2clm.nc4 AEROCLIMYEAR: 2008