Changeset 204 for LMDZ.3.3/trunk/libf/dyn3d/nudge.F

Timestamp:

Apr 13, 2001, 12:44:53 PM (23 years ago)

Author:

lmdz

Message:

Debogage du guidage et de la version debranchee et abandon de la version
debranchee non-netcdf FH/MAF
LF

File:

• LMDZ.3.3/trunk/libf/dyn3d/nudge.F (modified) (6 diffs)

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

@@ -99 +99 @@
      s   ,'1: en-ligne, 0: hors-ligne (x=x_rea), -1: climat (x=x_gcm)'
 cnec         read(*,*) online
-          online=-1
-         print*,'Entrer les constantes de temps de rappel en jours'
-         print*,'alphaT,alphau,alphav,alphaP'
+cnec      online=-1
+cnec     print*,'Entrer les constantes de temps de rappel en jours'
 cnec         read(*,*)alphaT
 cnec         read(*,*)alphau
@@ -110 +109 @@
              alphav=0.1
              alphaP=1.e10
+         print*,'alphaT,alphau,alphav,alphaP'
+     s          ,alphaT,alphau,alphav,alphaP
          if(online.eq.-1) return
          print*,'alpha rappel pour T, u, v, P ',
@@ -142 +143 @@
 c   Lecture du premier etat des reanalyses.
          call read_reanalyse(1
-     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,psrea2,1)
+     s   ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
 
 c-----------------------------------------------------------------------
@@ -171 +172 @@
            CALL SCOPY( ijp1llm,ucovrea2, 1, ucovrea1 , 1 )
            CALL SCOPY( ijp1llm,tetarea2,1,tetarea1   , 1 )
-           CALL SCOPY( ijp1llm,masserea2,1,masserea1   , 1 )
-           CALL SCOPY( ip1jmp1,psrea2, 1,  psrea1, 1 )
+c          CALL SCOPY( ijp1llm,masserea2,1,masserea1   , 1 )
+c          CALL SCOPY( ip1jmp1,psrea2, 1,  psrea1, 1 )
 
           print*,'LECTURE REANALYSES, pas ',step_rea
@@ -179 +180 @@
            itau_test=itau
            call read_reanalyse(step_rea
-     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,psrea2,1)
+     s     ,ucovrea2,vcovrea2,tetarea2,masserea2,ps,1)
         endif
       else
@@ -205 +206 @@
             a=(1.-tau)*tetarea1(ij,l)+tau*tetarea2(ij,l)
             teta(ij,l)=alphaT*teta(ij,l)+(1-alphaT)*a
-            a=(1.-tau)*masserea1(ij,l)+tau*masserea2(ij,l)
-            masse(ij,l)=alphaP*masse(ij,l)+(1-alphaP)*a
+c           a=(1.-tau)*masserea1(ij,l)+tau*masserea2(ij,l)
+c           masse(ij,l)=alphaP*masse(ij,l)+(1-alphaP)*a
          enddo
          do ij=1,ip1jm