Changeset 1246 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw.F90

Timestamp:

Sep 23, 2009, 4:47:40 PM (15 years ago)

Author:

Laurent Fairhead

Message:

• En deconnectant les aérosols (ok_ade=ok_aie=n) on a les mêmes

résultats avant et après les modifs.

• preindustrial readin fields are used to compute natural aerosol fields

to allow for clean double calls to radiation

• full forcing diagnostics (NAT, ANT, ZERO, Cloud forcing, CS,AS) are

activated with lev_histmth 4, If lev_histmth is not 4, the call to the
radiation is minimized, for efficiency, but ade and aie are computed and
applied (however for species wise forcing one would need to do
difference runs) (still quite a bit new forcing info, requires probably

some more explanation)

• there is a hardcoded key in sw_aeroAR4.F90 which lets you choose to use the zero aerosol, or natural aerosol perturbation acting on the meteorology, but still would put out the full forcing diagnostics.
• aod fields from offline aerosol fields are also output in histmth for

all aerosol tracers read in and available for evaluation

• aeropt contains the ss humidity correction from nicolas&yves

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw.F90 (modified) (5 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw.F90

@@ -18 +18 @@
    topswad0_aero, solswad0_aero,&
    topsw_aero, topsw0_aero,&
-   solsw_aero, solsw0_aero)
+   solsw_aero, solsw0_aero, &
+   topswcf_aero, solswcf_aero)
 
 
@@ -135 +136 @@
   REAL, DIMENSION(kdlon,9), INTENT(out) :: solsw_aero
   REAL, DIMENSION(kdlon,9), INTENT(out) :: solsw0_aero
+  REAL, DIMENSION(kdlon,3), INTENT(out) :: topswcf_aero
+  REAL, DIMENSION(kdlon,3), INTENT(out) :: solswcf_aero
 
 ! Local variables
@@ -181 +184 @@
   REAL(KIND=8) ztopsw_aero(kdlon,9), ztopsw0_aero(kdlon,9)
   REAL(KIND=8) zsolsw_aero(kdlon,9), zsolsw0_aero(kdlon,9)
+  REAL(KIND=8) ztopswcf_aero(kdlon,3), zsolswcf_aero(kdlon,3)     
   real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
 
@@ -338 +342 @@
                ztopsw_aero,ztopsw0_aero,&
                zsolsw_aero,zsolsw0_aero,&
+               ztopswcf_aero,zsolswcf_aero, & 
                ok_ade, ok_aie) 
           
@@ -384 +389 @@
           solswad_aero(iof+i) = zsolswadaero(i)
           solswad0_aero(iof+i) = zsolswad0aero(i)
-          topsw_aero(iof+i,:) = ztopsw_aero(iof+i,:)
-          topsw0_aero(iof+i,:) = ztopsw0_aero(iof+i,:)
-          solsw_aero(iof+i,:) = zsolsw_aero(iof+i,:)
-          solsw0_aero(iof+i,:) = zsolsw0_aero(iof+i,:)
-          
+! MS the following lines seem to be wrong, why is iof on right hand side???
+!          topsw_aero(iof+i,:) = ztopsw_aero(iof+i,:)
+!          topsw0_aero(iof+i,:) = ztopsw0_aero(iof+i,:)
+!          solsw_aero(iof+i,:) = zsolsw_aero(iof+i,:)
+!          solsw0_aero(iof+i,:) = zsolsw0_aero(iof+i,:)
+          topsw_aero(iof+i,:) = ztopsw_aero(i,:)
+          topsw0_aero(iof+i,:) = ztopsw0_aero(i,:)
+          solsw_aero(iof+i,:) = zsolsw_aero(i,:)
+          solsw0_aero(iof+i,:) = zsolsw0_aero(i,:)
+          topswcf_aero(iof+i,:) = ztopswcf_aero(i,:)
+          solswcf_aero(iof+i,:) = zsolswcf_aero(i,:)          
         ENDDO
     ELSE