Changeset 1179 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d

Timestamp:

Jun 11, 2009, 4:18:47 PM (16 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.

Location:

LMDZ4/branches/LMDZ4-dev/libf/dyn3d

Files:

• infotrac.F90 (modified) (2 diffs)
• pres2lev.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/infotrac.F90

@@ +1 @@
+! $Id$
+!
 MODULE infotrac
 
@@ -299 +301 @@
     END DO
 
+!
+! Test for advection schema. 
+! This version of LMDZ only garantees iadv=10 and iadv=14 (14 only for water vapour) .
+!
+    DO iq=1,nqtot
+       IF (iadv(iq)/=10 .AND. iadv(iq)/=14) THEN
+          WRITE(lunout,*)'STOP : The option iadv=',iadv(iq),' is not tested in this version of LMDZ'
+          CALL abort_gcm('infotrac_init','In this version only iadv=10 and iadv=14 is tested!',1)
+       ELSE IF (iadv(iq)==14 .AND. iq/=1) THEN
+          WRITE(lunout,*)'STOP : The option iadv=',iadv(iq),' is not tested in this version of LMDZ'
+          CALL abort_gcm('infotrac_init','In this version iadv=14 is only permitted for water vapour!',1)
+       END IF
+    END DO
+
 !-----------------------------------------------------------------------
 ! Finalize :

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/pres2lev.F

@@ -1 @@
-!
-! $Header$
+! $Id$
 !
 c******************************************************
@@ -21 +20 @@
 c  ARGUMENTS
 c  """""""""
-       LOGICAL ok_invertp
-       INTEGER lmo ! dimensions ancienne couches (input)
-       INTEGER lmn ! dimensions nouvelle couches (input)
-       INTEGER lmomx ! dimensions ancienne couches (input)
-       INTEGER lmnmx ! dimensions nouvelle couches (input)
+       LOGICAL, INTENT(IN) :: ok_invertp
+       INTEGER, INTENT(IN) :: lmo ! dimensions ancienne couches
+       INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches
+       INTEGER lmomx ! dimensions ancienne couches
+       INTEGER lmnmx ! dimensions nouvelle couches
 
        parameter(lmomx=10000,lmnmx=10000)
 
-        real po(ni,nj,lmo)! niveau de pression ancienne grille
-        real pn(ni,nj,lmn) ! niveau de pression nouvelle grille
+        real, INTENT(IN) :: po(ni,nj,lmo) ! niveau de pression ancienne grille
+        real, INTENT(IN) :: pn(ni,nj,lmn) ! niveau de pression nouvelle grille
 
-       INTEGER i,j,Nhoriz,ni,nj ! nombre de point horizontale (input)
+       INTEGER, INTENT(IN) :: ni,nj ! nombre de point horizontale
 
-       REAL varo(ni,nj,lmo) ! var dans l'ancienne grille (input)
-       REAL varn(ni,nj,lmn) ! var dans la nouvelle grille (output)
+       REAL, INTENT(IN)  :: varo(ni,nj,lmo) ! var dans l'ancienne grille
+       REAL, INTENT(OUT) :: varn(ni,nj,lmn) ! var dans la nouvelle grille
 
        real zvaro(lmomx),zpo(lmomx)
@@ -41 +40 @@
 c Autres variables
 c """"""""""""""""
-       INTEGER n, ln ,lo 
+       INTEGER n, ln ,lo, i, j, Nhoriz
        REAL coef