Changeset 1179 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

Timestamp:

Jun 11, 2009, 4:18:47 PM (15 years ago)

Author:

jghattas

Message:

• Ajout de l'interpolation vertical pour les nouveaux fichiers de forcage des aerosols. Utilisant les anciennes fichiers de SO4 pas d'interpolation possible. Convergence numerique avec la version precedente en utilisant les anciens fichiers des SO4. aerosol_optic.F90 change du nom pour readaerosol_optic.F90 (lecture d'aerosol + optic) Les fichiers de forcage aerosol doit maintenant avoir le suffix .nc.

• Correction des bugs pour inca et certain diagnostiques optionelles de radlwsw.
• Ajout de test pour le choix advection schema.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F (modified) (9 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

@@ -1 @@
-!
 ! $Id$
 !
@@ -537 +536 @@
 c
 c Variables propres a la physique
-c
-c      INTEGER radpas
-c      SAVE radpas                 ! frequence d'appel rayonnement
-ccccccccc$OMP THREADPRIVATE(radpas)
-c
-cc      INTEGER iflag_con
-c
       INTEGER itap
       SAVE itap                   ! compteur pour la physique
@@ -1074 +1066 @@
       CHARACTER*4, DIMENSION(9)      :: rfname 
       REAL, DIMENSION(klon)          :: aerindex     ! POLDER aerosol index
-      REAL, DIMENSION(klon,klev)     :: maerosol     ! aerosol concentration [ug/m3]
-      REAL, DIMENSION(klon,klev)     :: maerosol_pi  ! aerosol concentration [ug/m3] (pre-industrial value)
-      REAL, DIMENSION(klon,klev,9,2) :: tau_aero, piz_aero, cg_aero
-      REAL, DIMENSION(klon)          :: topswad_aero, solswad_aero   ! diag
-      REAL, DIMENSION(klon)          :: topswai_aero, solswai_aero   ! diag
-      REAL, DIMENSION(klon)          :: topswad0_aero, solswad0_aero ! pas utilise, eventuellment pour diag
-      REAL, DIMENSION(klon,9)        :: topsw_aero, solsw_aero       ! pas utilise
-      REAL, DIMENSION(klon,9)        :: topsw0_aero, solsw0_aero     ! pas utilise
-
+      REAL, DIMENSION(klon,klev)     :: mass_ins_aero! total mass concentration for all indissoluble aerosols[ug/m3]
+      REAL, DIMENSION(klon,klev)     :: mass_ins_aero_pi  ! - " - (pre-industrial value)
 
       ! Parameters
@@ -1228 +1213 @@
          tau_overturning_th(:)=0.
 
-         IF (config_inca /= 'none') ccm(:,:,:) = 0.
+         IF (config_inca /= 'none') THEN
+            ! jg : initialisation jusqu'au ces variables sont dans restart
+            ccm(:,:,:) = 0.
+            tau_aero(:,:,:,:) = 0.
+            piz_aero(:,:,:,:) = 0.
+            cg_aero(:,:,:,:) = 0.
+         END IF
 
          rnebcon0(:,:) = 0.0
@@ -2644 +2635 @@
       IF (ok_ade.OR.ok_aie) THEN
          IF (.NOT. aerosol_couple)
-     &        CALL aerosol_optic(
+     &        CALL readaerosol_optic(
      &        debut, new_aod, flag_aerosol, rjourvrai, pdtphys,
      &        pplay, paprs, t_seri, rhcl, 
-     &        maerosol, maerosol_pi, 
+     &        mass_ins_aero, mass_ins_aero_pi, 
      &        tau_aero, piz_aero, cg_aero )
       ELSE
@@ -2829 +2820 @@
 
       IF (aerosol_couple) THEN 
-         maerosol(:,:)    = ccm(:,:,1) 
-         maerosol_pi(:,:) = ccm(:,:,2) 
+         mass_ins_aero(:,:)    = ccm(:,:,1) 
+         mass_ins_aero_pi(:,:) = ccm(:,:,2) 
       END IF
 
@@ -2839 +2830 @@
      .            flwp, fiwp, flwc, fiwc,
      e            ok_aie,
-     e            maerosol, maerosol_pi,
+     e            mass_ins_aero, mass_ins_aero_pi,
      e            bl95_b0, bl95_b1,
      s            cldtaupi, re, fl)
@@ -2847 +2838 @@
      .            cldh, cldl, cldm, cldt, cldq,
      e            ok_aie,
-     e            maerosol, maerosol_pi,
+     e            mass_ins_aero, mass_ins_aero_pi,
      e            bl95_b0, bl95_b1,
      s            cldtaupi, re, fl)
@@ -2922 +2913 @@
          
 
-      ENDIF
+      ENDIF ! aerosol_couple
       itaprad = 0
-      ENDIF
+      ENDIF ! MOD(itaprad,radpas)
       itaprad = itaprad + 1