Changeset 431

Timestamp:

Dec 19, 2002, 3:52:05 PM (22 years ago)

Author:

lmdzadmin

Message:

Convergence avec la version de Ionela Musat

bilan_dyn.F : correction bug dimension (PF)
conf_gcm.F : remize a zero de la date initiale (LF)
control.h : remize a zero de la date initiale (LF)
diagedyn.F : correction bug dimension (IM)
etat0_netcdf.F : lecture de co2_ppm & solaire ds physiq.def
gcm.F : remize a zero de la date initiale (LF)

Location:

LMDZ.3.3/branches/rel-LF/libf/dyn3d

Files:

• bilan_dyn.F (modified) (1 diff)
• conf_gcm.F (modified) (1 diff)
• control.h (modified) (1 diff)
• diagedyn.F (modified) (1 diff)
• etat0_netcdf.F (modified) (2 diffs)
• gcm.F (modified) (1 diff)

LMDZ.3.3/branches/rel-LF/libf/dyn3d/bilan_dyn.F

@@ -41 +41 @@
       real phi(iip1,jjp1,llm)
       real ucov(iip1,jjp1,llm)
-      real vcov(iip1,jjp1,llm)
+      real vcov(iip1,jjm,llm)
       real trac(iip1,jjp1,llm,ntrac)
 

LMDZ.3.3/branches/rel-LF/libf/dyn3d/conf_gcm.F

@@ -94 +94 @@
       anneeref = 1998
       CALL getin('anneeref',anneeref)
+
+!Config  Key  = raz_date
+!Config  Desc = Remise a zero de la date initiale
+!Config  Def  = 0 (pas de remise a zero)
+!Config  Help = Remise a zero de la date initiale 
+!Config         0 pas de remise a zero, on garde la date du fichier restart
+!Config         1 prise en compte de la date de gcm.def avec remise a zero
+!Config         des compteurs de pas de temps
+      raz_date = 0
+      CALL getin('raz_date', raz_date)
 
 !Config  Key  = nday

LMDZ.3.3/branches/rel-LF/libf/dyn3d/control.h

@@ -4 +4 @@
       COMMON/control/nday,day_step,
      $              iperiod,iconser,iecri,idissip,iphysiq ,
-     $              periodav,ecritphy,iecrimoy,dayref,anneeref
+     $              periodav,ecritphy,iecrimoy,dayref,anneeref,
+     $              raz_date
 
       INTEGER   nday,day_step,iperiod,iconser,iecri,
-     $          idissip,iphysiq,iecrimoy,dayref,anneeref
+     $          idissip,iphysiq,iecrimoy,dayref,anneeref, raz_date
       REAL periodav, ecritphy
 

LMDZ.3.3/branches/rel-LF/libf/dyn3d/diagedyn.F

@@ -67 +67 @@
       REAL ps(ip1jmp1)                       ! pression  au sol
       REAL p (ip1jmp1,llmp1  )  ! pression aux interfac.des couches
-      REAL pk (ip1jmp1,llmp1  )  ! = (p/Pref)**kappa
+      REAL pk (ip1jmp1,llm  )  ! = (p/Pref)**kappa
       REAL teta(ip1jmp1,llm)                 ! temperature potentielle 
       REAL q(ip1jmp1,llm)               ! champs eau vapeur

LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F

@@ -92 +92 @@
       REAL :: pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
       REAL :: w(ip1jmp1,llm)
-      REAL ::phystep,co2_ppm,solaire
+cIM   REAL ::phystep,co2_ppm,solaire
+      REAL ::phystep
       INTEGER :: radpas
        real zrel(iip1*jjp1),chmin,chmax
@@ -613 +614 @@
       radpas    = NINT (86400./phystep/ FLOAT(nbapp_rad) )
       write(*,*)'phystep =', phystep, radpas
-      co2_ppm   = 330.0
-      solaire   = 1370.0
+cIM : lecture de co2_ppm & solaire ds physiq.def
+c     co2_ppm   = 348.0
+c     solaire   = 1365.0
 
 c

LMDZ.3.3/branches/rel-LF/libf/dyn3d/gcm.F

@@ -272 +272 @@
          PRINT*,'WARNING!!! changement de pas de temps',dtvr,'>',zdtvr
         ENDIF
+C
+C on remet le calendrier à zero
+c
+      if (annee_ref .ne. anneeref .or. day_ref .ne. dayref) then
+        write(*,*)' Attention les dates initiales lues dans le fichier'
+        write(*,*)' restart ne correspondent pas a celles lues dans '
+        write(*,*)' gcm.def'
+        if (raz_date .ne. 1) then
+          write(*,*)' On garde les dates du fichier restart'
+        else
+          annee_ref = anneeref
+          day_ref = dayref
+          itau_dyn = 0
+          itau_phy = 0
+          write(*,*)' On reinitialise a la date lue dans gcm.def'
+        endif
+      endif
 
 c  nombre d'etats dans les fichiers demarrage et histoire