Changeset 1246 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h

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/phys_output_write.h (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h

@@ -828 +828 @@
        ENDIF
 
+! OD550 per species
+      IF (new_aod) THEN
+      DO naero = 1, naero_spc
+          IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
+             CALL histwrite_phy(nid_files(iff),
+     $            o_tausumaero(naero)%name,itau_w,
+     $            tausum_aero(:,2,naero) )
+          ENDIF
+      END DO
+      ENDIF
+
        IF (ok_ade) THEN
           IF (o_topswad%flag(iff)<=lev_files(iff)) THEN
@@ -837 +848 @@
      $            solswad_aero)
           ENDIF
+
+!====MS forcing diagnostics
+        if (new_aod) then              
+        IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoaas_nat%name,itau_w,
+     $      topsw_aero(:,1))
+        ENDIF
+
+        IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfas_nat%name,itau_w,
+     $      solsw_aero(:,1))
+        ENDIF
+
+        IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoacs_nat%name,itau_w,
+     $      topsw0_aero(:,1))
+        ENDIF
+
+        IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfcs_nat%name,itau_w,
+     $      solsw0_aero(:,1))
+        ENDIF
+  
+!ant
+        IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoaas_ant%name,itau_w,
+     $      topsw_aero(:,2))
+        ENDIF
+
+        IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfas_ant%name,itau_w,
+     $      solsw_aero(:,2))
+        ENDIF
+
+        IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoacs_ant%name,itau_w,
+     $      topsw0_aero(:,2))
+        ENDIF
+
+        IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfcs_ant%name,itau_w,
+     $      solsw0_aero(:,2))
+        ENDIF
+
+!cf
+        IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
+     $      topswcf_aero(:,1))
+        ENDIF
+
+        IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfcf_nat%name,itau_w,
+     $      solswcf_aero(:,1))
+        ENDIF
+
+        IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoacf_ant%name,itau_w,
+     $      topswcf_aero(:,2))
+        ENDIF
+
+        IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfcf_ant%name,itau_w,
+     $      solswcf_aero(:,2))
+        ENDIF
+
+        IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swtoacf_zero%name,itau_w,
+     $      topswcf_aero(:,3))
+        ENDIF
+
+        IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN
+        CALL histwrite_phy(nid_files(iff),o_swsrfcf_zero%name,itau_w,
+     $      solswcf_aero(:,3))
+        ENDIF
+
+        endif ! new_aod
+!====MS forcing diagnostics
+
        ENDIF