Changeset 1150 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Apr 17, 2009, 5:34:01 PM (16 years ago)

Author:

jghattas

Message:

Ajoute possiblilite de forcer avec d'autre aersols avec le meme principe que pour les so4.

Par default la convergence est rempu avec flag_aersosol=1 + ok_ade=ok_aie=y. Exactement les memes resultats peuvent etre retrouver avec new_aod=n. Les resultats sont exactement les memes sans prise en compte aucun aerosol, avec ok_ade=ok_aie=n.

• flag_aerosol indique quel aersosol ou combinasion a utiliser (=1 uniquement les SO4 comme avant)
• les fichiers d'entree so4.runXXXX.cdf change du nom pour majescule SO4.runXXXX.cdf

• readsulfate change du nom pour readaerosol qui trait tous les aerosols
• radlwsw change du nom pour radlwsw_aero.
• aeropt_2bands.F90 et aeropt_5wv.F90 correspond a un reecriture de aeropt.F (premier et deuxieme moitie respectivement)
• aeropt.F est gardé pour retrouver la convergence si demande avec new_aod=false (=true default)
• sw_aero est un version evolue de sw_lmdar4 (dans fichier radiation_ar4.F)

ACo + JG

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_state_var_mod.F90 (modified) (7 diffs)

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

@@ -255 +255 @@
 !$OMP THREADPRIVATE(snow_con)
 !
-! sulfate_pi : SO4 aerosol concentration [ug/m3] (pre-industrial value)
-
-      REAL,SAVE,ALLOCATABLE :: sulfate_pi(:, :)
-!$OMP THREADPRIVATE(sulfate_pi)
       REAL,SAVE,ALLOCATABLE :: rlonPOS(:)
 !$OMP THREADPRIVATE(rlonPOS)
@@ -269 +265 @@
 ! ok_aie=T ->
 !       ok_ade=T -AIE=topswai-topswad
-!        ok_ade=F -AIE=topswai-topsw
+!       ok_ade=F -AIE=topswai-topsw
 !
 !topswad, solswad : Aerosol direct effect
@@ -277 +273 @@
       REAL,SAVE,ALLOCATABLE :: topswai(:), solswai(:)
 !$OMP THREADPRIVATE(topswai,solswai)
-      REAL,SAVE,ALLOCATABLE :: tau_ae(:,:,:), piz_ae(:,:,:)
-!$OMP THREADPRIVATE(tau_ae,piz_ae)
-      REAL,SAVE,ALLOCATABLE :: cg_ae(:,:,:)
-!$OMP THREADPRIVATE(cg_ae)
-
-! Les variables suivants uniquement pour un configuration avec INCA
-! topswad_inca, solswad_inca : Aerosol direct effect
-      REAL,SAVE,ALLOCATABLE :: topswad_inca(:), solswad_inca(:)
-!$OMP THREADPRIVATE(topswad_inca,solswad_inca)
-! topswad0_inca, solswad0_inca : Aerosol direct effect
-      REAL,SAVE,ALLOCATABLE :: topswad0_inca(:), solswad0_inca(:)
-!$OMP THREADPRIVATE(topswad0_inca,solswad0_inca)
-! topswai_inca, solswai_inca : Aerosol indirect effect
-      REAL,SAVE,ALLOCATABLE :: topswai_inca(:), solswai_inca(:)
-!$OMP THREADPRIVATE(topswai_inca,solswai_inca)
-      REAL,SAVE,ALLOCATABLE :: topsw_inca(:,:), solsw_inca(:,:)
-!$OMP THREADPRIVATE(topsw_inca,solsw_inca)
-      REAL,SAVE,ALLOCATABLE :: topsw0_inca(:,:), solsw0_inca(:,:)
-!$OMP THREADPRIVATE(topsw0_inca,solsw0_inca)
-      REAL,SAVE,ALLOCATABLE :: tau_inca(:,:,:,:)
-!$OMP THREADPRIVATE(tau_inca)
-      REAL,SAVE,ALLOCATABLE :: piz_inca(:,:,:,:)
-!$OMP THREADPRIVATE(piz_inca)
-      REAL,SAVE,ALLOCATABLE :: cg_inca(:,:,:,:)
-!$OMP THREADPRIVATE(cg_inca)
+
       REAL,SAVE,ALLOCATABLE :: ccm(:,:,:)
 !$OMP THREADPRIVATE(ccm)
@@ -402 +374 @@
       ALLOCATE(ibas_con(klon), itop_con(klon))
       ALLOCATE(rain_con(klon), snow_con(klon))
-!
-      ALLOCATE(sulfate_pi(klon, klev))
       ALLOCATE(rlonPOS(klon))
       ALLOCATE(newsst(klon))
@@ -409 +379 @@
       ALLOCATE(topswad(klon), solswad(klon))
       ALLOCATE(topswai(klon), solswai(klon))
-      ALLOCATE(tau_ae(klon,klev,2), piz_ae(klon,klev,2))
-      ALLOCATE(cg_ae(klon,klev,2))
-
-      IF (config_inca /= 'none') THEN
-         ALLOCATE(topswad_inca(klon), solswad_inca(klon))
-         ALLOCATE(topswad0_inca(klon), solswad0_inca(klon))
-         ALLOCATE(topswai_inca(klon), solswai_inca(klon))
-         ALLOCATE(topsw_inca(klon,9), solsw_inca(klon,9))
-         ALLOCATE(topsw0_inca(klon,9), solsw0_inca(klon,9))
-      END IF
-      ! Following 4 variables are needed only by INCA but must be
-      ! allocated as they exist in the phytrac argument list
-      ALLOCATE(tau_inca(klon,klev,9,2))
-      ALLOCATE(piz_inca(klon,klev,9,2))
-      ALLOCATE(cg_inca(klon,klev,9,2))
+
       ALLOCATE(ccm(klon,klev,2))
 
@@ -505 +461 @@
       deallocate(ibas_con, itop_con)
       deallocate(rain_con, snow_con)
-!
-      deallocate(sulfate_pi)
       deallocate(rlonPOS)
       deallocate(newsst)
@@ -513 +467 @@
       deallocate(topswai, solswai)
 
-      deallocate(tau_ae, piz_ae)
-      deallocate(cg_ae)
-
-      IF (config_inca /= 'none') THEN
-         deallocate(topswad_inca, solswad_inca)
-         deallocate(topswad0_inca, solswad0_inca)
-         deallocate(topswai_inca, solswai_inca)
-         deallocate(topsw_inca, solsw_inca)
-         deallocate(topsw0_inca, solsw0_inca)
-      END IF
-      deallocate(tau_inca)
-      deallocate(piz_inca)
-      deallocate(cg_inca)
       deallocate(ccm)