problem about the modification of radiation.F90

[arthur_panch@gmail_com]

Hi,
Now I want to turn off the radiative effect of ozone over the Tibet Plateau (TP). I define the TP region in the area of 25°~40°N, 75°~105°E. I add several lines after the line “call radinp(ncol, state%pmid, state%pint, pbr, pnm, eccf) ” in the radiation.F90 file. The modification is as follows:
!==================================================
pi=pi/180.0_r8
do k=1, pver
do i=1, ncol
if ( state%lat(i) >= 25.0_r8*pi .and. state%lat(i) = 75.0_r8*pi .and. state%lon(i)

 

[mmills]

Arthur,There are 2 copies of radiation.F90. The one in physics/cam is used for CAM-RT radiation, which is used in CAM4 and WACCM4. The one in physics/rrtmg is used in CAM5 and WACCM5. If you are modifying the wrong one, you will see no effect.If you are not interested in chemistry interactions, you may want to run WACCM with specified chemistry. We have several compsets for SC-WACCM. This would allow you to prescribe O3 from an input file based on a previous WACCM run to acheive what you wish to do.

 

[arthur_panch@gmail_com]

Mills,Thank you very much. I have modified the wrong file.