Changeset 1150 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/conf_phys.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/conf_phys.F90 (modified) (12 diffs)

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

@@ -12 +12 @@
  &                     iflag_cldcon, &
  &                     iflag_ratqs,ratqsbas,ratqshaut, &
- &                     ok_ade, ok_aie, aerosol_couple, &
+ &                     ok_ade, ok_aie, aerosol_couple, &
+ &                     flag_aerosol, new_aod, &
  &                     bl95_b0, bl95_b1,&
  &                     iflag_thermals,nsplit_thermals,tau_thermals, &
@@ -60 +61 @@
   logical              :: ok_LES
   LOGICAL              :: ok_ade, ok_aie, aerosol_couple
+  INTEGER              :: flag_aerosol
+  LOGICAL              :: new_aod
   REAL                 :: bl95_b0, bl95_b1
   real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
@@ -71 +74 @@
   logical,SAVE        :: ok_LES_omp   
   LOGICAL,SAVE        :: ok_ade_omp, ok_aie_omp, aerosol_couple_omp
+  INTEGER, SAVE       :: flag_aerosol_omp
+  LOGICAL, SAVE       :: new_aod_omp
   REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
   REAL,SAVE           :: freq_ISCCP_omp, ecrit_ISCCP_omp
@@ -239 +244 @@
   CALL getin('aerosol_couple',aerosol_couple_omp)
 
+!
+!Config Key  = flag_aerosol
+!Config Desc = which aerosol is use for coupled model
+!Config Def  = 1
+!Config Help = Used in physiq.F
+!
+! - flag_aerosol=1 => so4 only (defaut) 
+! - flag_aerosol=2 => bc  only 
+! - flag_aerosol=3 => pom only
+! - flag_aerosol=4 => all aerosol 
+! - flag_aerosol=5 => bcpom
+! - flag_aerosol=6 => pomsulf
+
+  flag_aerosol_omp = 1
+  CALL getin('flag_aerosol',flag_aerosol_omp)
+
+! Temporary variable for testing purpose!!
+!Config Key  = new_aod
+!Config Desc = which calcul of aeropt
+!Config Def  = false
+!Config Help = Used in physiq.F
+!
+  new_aod_omp = .true.
+  CALL getin('new_aod',new_aod_omp)
+
 ! 
 !Config Key  = aer_type 
 !Config Desc = Use a constant field for the aerosols 
 !Config Def  = scenario 
-!Config Help = Used in readsulfate.F 
+!Config Help = Used in readaerosol.F90 
 ! 
   aer_type_omp = 'scenario' 
@@ -918 +948 @@
 !Config Help =
 !
-  iflag_coupl = 0
+  iflag_coupl_omp = 0
   call getin('iflag_coupl',iflag_coupl_omp)
 
@@ -927 +957 @@
 !Config Help = 
 !
-  iflag_clos = 1
+  iflag_clos_omp = 1
   call getin('iflag_clos',iflag_clos_omp)
 !
@@ -935 +965 @@
 !Config Help = 
 !
-  iflag_cvl_sigd = 0
+  iflag_cvl_sigd_omp = 0
   call getin('iflag_cvl_sigd',iflag_cvl_sigd_omp)
 
@@ -943 +973 @@
 !Config Help = 
 !
-  iflag_wake = 0
+  iflag_wake_omp = 0
   call getin('iflag_wake',iflag_wake_omp)
 
@@ -1266 +1296 @@
     ok_aie = ok_aie_omp
     aerosol_couple = aerosol_couple_omp
+    flag_aerosol=flag_aerosol_omp
+    new_aod=new_aod_omp
     aer_type = aer_type_omp
     bl95_b0 = bl95_b0_omp
@@ -1326 +1358 @@
        WRITE(numout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
        CALL abort_gcm('conf_phys','version_ocean not valid',1)
+    END IF
+
+! Test sur new_aod. Ce flag permet de retrouver les resultats de l'AR4
+! il n'est utilisable que lors du couplage avec le SO4 seul 
+    IF (ok_ade .OR. ok_aie) THEN 
+       IF ( .NOT. new_aod .AND.  flag_aerosol .NE. 1) THEN
+          CALL abort_gcm('conf_phys','new_aod=.FALSE. not compatible avec flag_aerosol=1',1)
+       END IF
     END IF
 
@@ -1396 +1436 @@
   write(numout,*)' ok_aie = ',ok_aie
   write(numout,*)' aerosol_couple = ', aerosol_couple
+  write(numout,*)' flag_aerosol = ', flag_aerosol
+  write(numout,*)' new_aod = ', new_aod
   write(numout,*)' aer_type = ',aer_type
   write(numout,*)' bl95_b0 = ',bl95_b0
@@ -1409 +1451 @@
   write(numout,*)' type_run = ',type_run 
   write(numout,*)' ok_isccp = ',ok_isccp 
-  WRITE(numout,*)' solarlong0 = ', solarlong0
+  write(numout,*)' solarlong0 = ', solarlong0
   write(numout,*)' qsol0 = ', qsol0
   write(numout,*)' inertie_sol = ', inertie_sol