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

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/dyn3dpar

Files:

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

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

@@ +1 @@
+! $Id$
+!
 MODULE infotrac
 
 ! nqtot : total number of tracers and higher order of moment, water vapor and liquid included
   INTEGER, SAVE :: nqtot
-!!$OMP THREADPRIVATE(nqtot)   
 
 ! nbtr : number of tracers not including higher order of moment or water vapor or liquid
 !        number of tracers used in the physics
   INTEGER, SAVE :: nbtr
-!!$OMP THREADPRIVATE(nbtr)   
 
 ! Name variables
   CHARACTER(len=20), ALLOCATABLE, DIMENSION(:), SAVE :: tname ! tracer short name for restart and diagnostics
   CHARACTER(len=23), ALLOCATABLE, DIMENSION(:), SAVE :: ttext ! tracer long name for diagnostics
-!!$OMP THREADPRIVATE(tname,ttext)   
 
 ! iadv  : index of trasport schema for each tracer
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE    :: iadv
-!!$OMP THREADPRIVATE(iadv)   
 
 ! niadv : vector keeping the coorspondance between all tracers(nqtot) treated in the 
 !         dynamic part of the code and the tracers (nbtr+2) used in the physics part of the code. 
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE    :: niadv ! equivalent dyn / physique
-!!$OMP THREADPRIVATE(niadv)   
 
 ! Variables for INCA
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: conv_flg
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE  :: pbl_flg
-!!$OMP THREADPRIVATE(conv_flg, pbl_flg)   
 
 CONTAINS
@@ -305 +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/dyn3dpar/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