Changeset 618 for LMDZ4/trunk/libf/dyn3d

Timestamp:

Apr 15, 2005, 2:31:21 PM (20 years ago)

Author:

lmdzadmin

Message:

Sauvegarde de ncidpl MAF+AnneC

File:

• LMDZ4/trunk/libf/dyn3d/read_reanalyse.F (modified) (5 diffs)

LMDZ4/trunk/libf/dyn3d/read_reanalyse.F

@@ -43 +43 @@
       integer varidpl,ncidQ,varidQ
       save ncidu,varidu,ncidv,varidv,ncidt,varidt,ncidps,varidps
+      save ncidpl
       save varidpl,ncidQ,varidQ
 
@@ -48 +49 @@
       real*4 tnc(iip1,jjp1,nlevnc),psnc(iip1,jjp1)
       real*4 Qnc(iip1,jjp1,nlevnc)
-      real*4 pl(nlevnc),presnc(iip1,jjp1,nlevnc)
+      real*4 pl(nlevnc)
 
       integer start(4),count(4),status
@@ -106 +107 @@
             endif
 
-c Coordonnee vertcale
+c Coordonnee verticale
             if (ncep) then
                print*,'Vous etes entrain de lire des donnees NCEP'
@@ -115 +116 @@
             endif
             print*,'ncidu,varidpl',ncidu,varidpl
-
-      endif
-
+      endif
       print*,'ok1'
 
-c -----------------------------------------------------------------
-c   lecture des champs u, v, T, ps
-c -----------------------------------------------------------------
-
-c  niveaux de pression
-c  -------------------
-
+c Niveaux de pression
       print*,'WARNING!!! Il n y a pas de test de coherence'
       print*,'sur le nombre de niveaux verticaux dans le fichier nc'
@@ -136 +129 @@
 c  passage en pascal
       pl(:)=100.*pl(:)
-c  passage des presions sur une grille 3D
-      do l=1,nlevnc
-         presnc(:,:,l)=pl(l)
-      enddo
-      if(first) then
-         do l=1,nlevnc
-            print*,'PL(',l,')=',pl(l)
-         enddo
-      endif
-
+      if (first) then
+       do l=1,nlevnc
+          print*,'PL(',l,')=',pl(l)
+       enddo
+      endif
+
+c -----------------------------------------------------------------
+c   lecture des champs u, v, T, ps
+c -----------------------------------------------------------------
 
 c  dimensions pour les champs scalaires et le vent zonal