Changeset 205 for LMDZ.3.3

Timestamp:

Apr 13, 2001, 4:58:01 PM (23 years ago)

Author:

lmdzadmin

Message:

Rajout date0 dans physiq pour interface avec ORCHIDEE
Fermeture des fichiers restart par restclo dans abort_gcm
LF

Location:

LMDZ.3.3/branches/rel-LF

Files:

• libf/dyn3d/abort_gcm.F (modified) (2 diffs)
• libf/phylmd/clmain.F (modified) (7 diffs)
• libf/phylmd/interface_surf.F90 (modified) (12 diffs)
• libf/phylmd/physiq.F (modified) (8 diffs)
• makegcm (modified) (1 diff)

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

@@ +1 @@
+c
+c $Header
+c
       SUBROUTINE abort_gcm(modname, message, ierr)
      
@@ -16 +19 @@
       write(*,*) 'in abort_gcm'
       call histclo
+      call restclo
 c     call histclo(2)
 c     call histclo(3)

LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

@@ -1 @@
-      SUBROUTINE clmain(dtime,itap,pctsrf,
+      SUBROUTINE clmain(dtime,itap,date0,pctsrf,
      .                  t,q,u,v,
      .                  jour, rmu0,
@@ -36 +36 @@
 c dtime----input-R- interval du temps (secondes)
 c itap-----input-I- numero du pas de temps
+c date0----input-R- jour initial
 c t--------input-R- temperature (K)
 c q--------input-R- vapeur d'eau (kg/kg)
@@ -76 +77 @@
 c
       REAL dtime
+      real date0
       integer itap
       REAL t(klon,klev), q(klon,klev)
@@ -368 +370 @@
 
 c calculer la diffusion de "q" et de "h"
-      CALL clqh(dtime, itap, jour, debut,lafin,
+      CALL clqh(dtime, itap, date0,jour, debut,lafin,
      e          rlon, rlat, cufi, cvfi,
      e          knon, nsrf, ni, pctsrf,
@@ -479 +481 @@
       RETURN
       END
-      SUBROUTINE clqh(dtime,itime, jour,debut,lafin,
+      SUBROUTINE clqh(dtime,itime, date0,jour,debut,lafin,
      e                rlon, rlat, cufi, cvfi, 
      e                knon, nisurf, knindex, pctsrf,
@@ -512 +514 @@
       INTEGER knon
       REAL dtime              ! intervalle du temps (s)
+      real date0
       REAL u1lay(klon)        ! vitesse u de la 1ere couche (m/s)
       REAL v1lay(klon)        ! vitesse v de la 1ere couche (m/s)
@@ -729 +732 @@
       ccanopy = 365.
 
-      CALL interfsurf(itime, dtime, jour, rmu0,
+      CALL interfsurf(itime, dtime, date0, jour, rmu0,
      e klon, iim, jjm, nisurf, knon, knindex, pctsrf, 
      e rlon, rlat, cufi, cvfi, 

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -49 +49 @@
 !############################################################################
 !
-  SUBROUTINE interfsurf_hq(itime, dtime, jour, rmu0, &
+  SUBROUTINE interfsurf_hq(itime, dtime, date0, jour, rmu0, &
       & klon, iim, jjm, nisurf, knon, knindex, pctsrf, &
       & rlon, rlat, cufi, cvfi,&
@@ -77 +77 @@
 !   iim, jjm     nbres de pts de grille
 !   dtime        pas de temps de la physique (en s)
+!   date0        jour initial 
 !   jour         jour dans l'annee en cours,
 !   rmu0         cosinus de l'angle solaire zenithal
@@ -138 +139 @@
   integer, intent(IN) :: klon
   real, intent(IN) :: dtime
+  real, intent(IN) :: date0
   integer, intent(IN) :: jour
   real, intent(IN)    :: rmu0(klon)
@@ -345 +347 @@
       scatter= .true.
       if (.not. scatter) then
-        call interfsol(itime, klon, dtime, nisurf, knon, &
+        call interfsol(itime, klon, dtime, date0, nisurf, knon, &
      &  knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
      &  debut, lafin, ok_veget, &
@@ -355 +357 @@
      &  tsol_rad, tsurf_new, alb_new, emis_new, z0_new, dflux_l, dflux_s)
       else
-        call interfsol_scat(itime, klon, dtime, nisurf, knon, &
+        call interfsol_scat(itime, klon, dtime, date0, nisurf, knon, &
      &  knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
      &  debut, lafin, ok_veget, &
@@ -622 +624 @@
 !#########################################################################
 !
-  SUBROUTINE interfsol(itime, klon, dtime, nisurf, knon, &
+  SUBROUTINE interfsol(itime, klon, dtime, date0, nisurf, knon, &
      & knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
      & debut, lafin, ok_veget, &
@@ -694 +696 @@
   integer, intent(IN) :: klon
   real, intent(IN)    :: dtime
+  real, intent(IN)    :: date0
   integer, intent(IN) :: nisurf
   integer, intent(IN) :: knon
@@ -754 +757 @@
   real, dimension(klon):: qsurf
   real, dimension(klon):: snow, qsol
-  real, save                 :: date0 = 0.
   real, dimension(knon,2) :: albedo_out
 ! Pb de nomenclature
@@ -963 +965 @@
 !#########################################################################
 !
-  SUBROUTINE interfsol_scat(itime, klon, dtime, nisurf, knon, &
+  SUBROUTINE interfsol_scat(itime, klon, dtime, date0, nisurf, knon, &
      & knindex, rlon, rlat, cufi, cvfi, iim, jjm, pctsrf, &
      & debut, lafin, ok_veget, &
@@ -1035 +1037 @@
   integer, intent(IN) :: klon
   real, intent(IN)    :: dtime
+  real, intent(IN)    :: date0
   integer, intent(IN) :: nisurf
   integer, intent(IN) :: knon
@@ -1095 +1098 @@
   real, dimension(klon):: qsurf
   real, dimension(klon):: snow, qsol
-  real, save                 :: date0 = 0.
   real, dimension(knon,2) :: albedo_out
 ! Pb de nomenclature
@@ -1192 +1194 @@
 !
   if (debut) then
+    if (check) write(*,*) 'debut orchidee itime - 1', itime-1,date0
     call intersurf_main (itime-1, iim, jjm+1 , knon, knindex, dtime, &
      & lrestart_read, lrestart_write, lon_sc, lat_sc, &

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -94 +94 @@
 ccc      PARAMETER (soil_model=.FALSE.)
       logical ok_veget
-      parameter (ok_veget = .false.)
-c      parameter (ok_veget = .true.)
+      parameter (ok_veget = .true.)
+c      parameter (ok_veget = .false.)
 c======================================================================
 c Dans les versions precedentes, l'eau liquide nuageuse utilisee dans
@@ -118 +118 @@
 c
       LOGICAL ok_instan ! sortir le fichier instantane
-      PARAMETER (ok_instan=.true.)
+      PARAMETER (ok_instan=.false.)
 c
       LOGICAL ok_region ! sortir le fichier regional
@@ -499 +499 @@
       character*80 abort_message
       logical ok_sync
-
+      real date0
 
 c
@@ -1525 +1525 @@
 c Verifier les temperatures
 c
-
       CALL hgardfou(t_seri,ftsol,'debutphy')
 c
@@ -1593 +1592 @@
 
       fder = 0.
-
-      CALL clmain(dtime,itap,pctsrf,
+      date0 = day_ini
+
+      CALL clmain(dtime,itap,date0,pctsrf,
      e            t_seri,q_seri,u_seri,v_seri,
      e            julien, rmu0, 
@@ -1607 +1607 @@
      s            dsens, devap,
      s            ycoefh,yu1,yv1) 
+
 c
 C§§§ PB
@@ -1633 +1634 @@
 c         evap(i) = - fluxq(i,1) ! flux d'evaporation au sol
          evap(i) = - zxfluxq(i,1) ! flux d'evaporation au sol
+C LF test signe flux
+         sens(i) = zxfluxt(i,1)
+         evap(i) = zxfluxq(i,1)
          fder(i) = dsens(i) + devap(i)
       ENDDO
@@ -1738 +1742 @@
       CALL abort
       ENDIF
+
       CALL homogene(paprs, q_seri, d_q_con, u_seri,v_seri,
      .              d_u_con, d_v_con)
       DO k = 1, klev
-      DO i = 1, klon
+        DO i = 1, klon
          t_seri(i,k) = t_seri(i,k) + d_t_con(i,k)
          q_seri(i,k) = q_seri(i,k) + d_q_con(i,k)
          u_seri(i,k) = u_seri(i,k) + d_u_con(i,k)
          v_seri(i,k) = v_seri(i,k) + d_v_con(i,k)
-      ENDDO
+        ENDDO
       ENDDO
       IF (check) THEN

LMDZ.3.3/branches/rel-LF/makegcm

@@ -77 +77 @@
 #Force/couple
 set COUPLE=0
-set ocean=`grep '^      parameter (ocean' libf/phylmd/physiq.F | awk -F\' '{print $2}'`
+set ocean=`grep '^      parameter (ocean' $LMDGCM/libf/phylmd/physiq.F | awk -F\' '{print $2}'`
 banner $ocean
 if ($ocean == couple) then
   set COUPLE=1
 endif
-
-
+#veget ou non
+set veget=`grep '^      parameter (ok_veget' $LMDGCM/libf/phylmd/physiq.F | awk -F. '{print $2}'`
+if ($veget == true) then
+  banner orchidee
+endif
 ########################################################################
 #  Les differentes platformes reconnues