Changeset 356 for LMDZ.3.3/trunk

Timestamp:

Apr 2, 2002, 4:17:16 PM (23 years ago)

Author:

lmdz

Message:

Lecture du fichier netcdf pour determination du nombre de niveaux verticaux MAF
LF

Location:

LMDZ.3.3/trunk/libf/dyn3d

Files:

• nudge.F (modified) (4 diffs)
• read_reanalyse.F (modified) (2 diffs)

LMDZ.3.3/trunk/libf/dyn3d/nudge.F

@@ -26 +26 @@
 
 c   ...  Possibilite de choisir le shema de Van-leer pour l'advection de
-c         q  , en faisant iadv = 3  dans   traceur  (29/04/97) .
+c         q  , en faisant iadv = 10  dans   traceur  (29/04/97) .
 c
 c-----------------------------------------------------------------------
@@ -70 +70 @@
 c
       INTEGER itau,ij,l
+      integer ncidt,varidpl,nlev,status
+      integer rcod,rid 
       real ditau,tau,a
-
-      LOGICAL first
+      save nlev
+
+      LOGICAL first,ncep
       integer online
-      save first,online
+      save first,online,ncep
       data first,online/.true.,1/
+      data ncep/.false./
 
       save ucovrea1,vcovrea1,tetarea1,masserea1,psrea1
@@ -139 +143 @@
 
 c    itau_test    montre si l'importation a deja ete faite au rang itau
-
-
+c lecture d'un fichier netcdf pour determiner le nombre de niveaux
+         ncidt=NCOPN('T.nc',NCNOWRIT,rcod)
+         if (ncep) then
+          status=NF_INQ_DIMID(ncidt,'LEVEL',rid)
+         else
+          status=NF_INQ_DIMID(ncidt,'PRESSURE',rid)
+         endif
+          status=NF_INQ_DIMLEN(ncidt,rid,nlev)
+         print *,'nlev', nlev 
+          call ncclos(ncidt,rcod)
 c   Lecture du premier etat des reanalyses.
          call read_reanalyse(1
-     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
+     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1,nlev)
 
 c-----------------------------------------------------------------------
@@ -180 +192 @@
            itau_test=itau
            call read_reanalyse(step_rea
-     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
+     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1,nlev)
         endif
       else

LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F

@@ -2 +2 @@
 c $Header$
 c
-      subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode)
+      subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode,nlevnc)
 
 c   mode=0 variables naturelles
@@ -15 +15 @@
 
       integer nlevnc
-      parameter (nlevnc=15)
-C pour annee 2000      parameter (nlevnc=21)
+c      parameter (nlevnc=15)
+C pour annee 2000      
+c      parameter (nlevnc=21)
       integer timestep,mode,l