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

Timestamp:

Jul 28, 2000, 2:38:04 PM (24 years ago)

Author:

lmdzadmin

Message:

Mise au point de l'interface en force, ca tourne sur un pas de temps
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (5 diffs)

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

@@ -70 +70 @@
       logical rnpb
       parameter(rnpb=.true.)
-      PARAMETER (npas=1440)
-      PARAMETER (nexca=48)
-      PARAMETER (itimestep=1800)
+c      PARAMETER (npas=1440)
+c      PARAMETER (nexca=48)
+c      PARAMETER (itimestep=1800)
       EXTERNAL fromcpl, intocpl, inicma
       REAL cpl_sst(iim,jjmp1), cpl_sic(iim,jjmp1)
       REAL cpl_alb_sst(iim,jjmp1), cpl_alb_sic(iim,jjmp1)
+      character *6 ocean
+      parameter (ocean = 'force ')
 c======================================================================
 c ok_ocean indique l'utilisation du modele oceanique "slab ocean",
@@ -669 +671 @@
          PRINT*, 'La frequence de sortie region est de ', ecrit_reg
          ENDIF
+
+c
+c Initialiser le couplage si necessaire
+c
+      npas = 0
+      nexca = 0
+      if (ocean == 'couple') then
+        npas = itaufin/ iphysiq
+        nexca = 86400 / dtime
+        write(*,*)' ##### Ocean couple #####'
+        write(*,*)' Valeurs des pas de temps'
+        write(*,*)' npas = ', npas
+        write(*,*)' nexca = ', nexca
+      endif        
 c
 c
@@ -1536 +1552 @@
       DO i = 1, klon
         if (.not. ok_veget) then
-         frugs(i,is_ter) = SQRT(frugs(i,is_ter)**2+rugoro(i)**2)
+          frugs(i,is_ter) = SQRT(frugs(i,is_ter)**2+rugoro(i)**2)
+        endif 
          frugs(i,is_lic) = rugoro(i)
          frugs(i,is_oce) = rugmer(i)
@@ -1563 +1580 @@
       ENDIF
 
-      CALL clmain(dtime,pctsrf,
+      CALL clmain(dtime,itap,pctsrf,
      e            t_seri,q_seri,u_seri,v_seri,
      e            julien, rmu0, 
-     e            ok_veget, ftsol,
+     e            ok_veget, ocean, npas, nexca, ftsol,
      e            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,
      e            rain_fall, snow_fall, solsw, sollw, fder,
      e            rlon, rlat, frugs,
-     e            debut, lafin,
+     e            debut, lafin, agesno,
      s            d_t_vdf,d_q_vdf,d_u_vdf,d_v_vdf,d_ts,
      s            fluxt,fluxq,fluxu,fluxv,cdragh,cdragm,rugmer,
@@ -2793 +2810 @@
      .      falbe, fevap, rain_fall, snow_fall,
      .      solsw, sollw,
-     .      radsol,rugmer,agesno,
+     .      radsol,frugs,agesno,
      .      zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro,
      .      t_ancien, q_ancien)