Changeset 617 for LMDZ4/trunk/libf/dyn3d/guide.F

Timestamp:

Apr 13, 2005, 3:53:20 PM (19 years ago)

Author:

lmdzadmin

Message:

corrections pour le mode guide :

• mode non zoome
• guidage de certains champs seulement

MAF

File:

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

LMDZ4/trunk/libf/dyn3d/guide.F

@@ -81 +81 @@
       INTEGER step_rea,count_no_rea
 
-      real aire_min,aire_max
       integer ilon,ilat
       real factt,ztau(ip1jmp1)
 
       INTEGER itau,ij,l,i,j
-      integer ncidt,varidpl,nlev,status
+      integer ncidpl,varidpl,nlev,status
       integer rcod,rid 
       real ditau,tau,a
@@ -208 +207 @@
          step_rea=1
          count_no_rea=0
+         ncidpl=-99
 
 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 (guide_u) then
+           if (ncidpl.eq.-99) ncidpl=NCOPN('u.nc',NCNOWRIT,rcod)
+         endif
+c
+         if (guide_v) then
+           if (ncidpl.eq.-99) ncidpl=NCOPN('v.nc',NCNOWRIT,rcod)
+         endif
+c
+         if (guide_T) then
+           if (ncidpl.eq.-99) ncidpl=NCOPN('T.nc',NCNOWRIT,rcod)
+         endif
+c
+         if (guide_Q) then
+           if (ncidpl.eq.-99) ncidpl=NCOPN('hur.nc',NCNOWRIT,rcod)
+         endif
+c
          if (ncep) then
-          status=NF_INQ_DIMID(ncidt,'LEVEL',rid)
+          status=NF_INQ_DIMID(ncidpl,'LEVEL',rid)
          else
-          status=NF_INQ_DIMID(ncidt,'PRESSURE',rid)
-         endif
-          status=NF_INQ_DIMLEN(ncidt,rid,nlev)
+          status=NF_INQ_DIMID(ncidpl,'PRESSURE',rid)
+         endif
+          status=NF_INQ_DIMLEN(ncidpl,rid,nlev)
          print *,'nlev', nlev 
-          call ncclos(ncidt,rcod)
+          call ncclos(ncidpl,rcod)
 c   Lecture du premier etat des reanalyses.
          call read_reanalyse(1,ps
@@ -296 +311 @@
       tau=tau-aint(tau)
 
-      print*,'ATTENTION !!!! ON NE GUIDE QUE JUSQU A 15N'
-
 c  ucov
       if (guide_u) then
@@ -485 +498 @@
                zlat=rlatv(j)*180./pi
             endif
+         if (abs(grossismx-1.).lt.0.1.or.abs(grossismy-1.).lt.0.1) then
+c  pour une grille reguliere, xi=xxx**0=1 -> alpha=alphamin
+             alpha(i,j)=alphamin
+         else
             xi=((dxdy_max-dxdy_)/(dxdy_max-dxdy_min))**gamma
-c  pour une grille reguliere, xi=xxx**0=1 -> alpha=alphamin
             xi=min(xi,1.)
             if(lat_min_guide.le.zlat .and. zlat.le.lat_max_guide) then
@@ -493 +509 @@
                alpha(i,j)=0.
             endif
+         endif
          enddo
       enddo