Changeset 101 for LMDZ.3.3/branches/rel-LF

Timestamp:

Jul 5, 2000, 6:15:11 PM (24 years ago)

Author:

lmdzadmin

Message:

Rajout interfoce_cpl. LF

Location:

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

Files:

• clmain.F (modified) (8 diffs)
• interface_surf.F90 (modified) (8 diffs)
• oasis.F (modified) (1 diff)
• physiq.F (modified) (1 diff)

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

@@ -2 +2 @@
      .                  ok_veget,ts,
      .                  paprs,pplay,radsol,snow,qsol,evap,albe,
-     .                  rain_f, snow_f, solsw, sollw,
+     .                  rain_f, snow_f, solsw, sollw, fder,
      .                  rlon, rlat, rugos,
      .                  debut, lafin,
@@ -99 +99 @@
       REAL albe(klon,nbsrf)
       real rain_f(klon), snow_f(klon)
+      REAL fder(klon)
       REAL sollw(klon), solsw(klon)
       REAL rugos(klon,nbsrf)
@@ -115 +116 @@
       real yrain_f(klon), ysnow_f(klon)
       real ysollw(klon), ysolsw(klon), ysolswnet(klon)
+      real yfder(klon), ytaux(klon), ytauy(klon)
       REAL yrugm(klon), yrads(klon)
       REAL y_d_ts(klon)
@@ -258 +260 @@
         yrain_f(j) = rain_f(i)
         ysnow_f(j) = snow_f(i)
+        yfder(j) = fder(i)
+        ytaux(j) = flux_u(i,nsrf)
+        ytauy(j) = flux_v(i,nsrf)
         ysolsw(j) = solsw(i)
         ysollw(j) = sollw(i)
@@ -302 +307 @@
      e          ycoefh,yt,yq,yts,ypaprs,ypplay,ydelp,yrads,
      e          yevap,yalb, ysnow, yqsol, yrain_f, ysnow_f,
+     e          yfder, ytaux, ytauy,
      e          ysollw, ysolsw,
      s          y_d_t, y_d_q, y_d_ts, 
@@ -423 +429 @@
      e                t,q,ts,paprs,pplay,
      e                delp,radsol,evap,albedo,snow,qsol, 
-     e                precip_rain, precip_snow,
+     e                precip_rain, precip_snow, fder, taux, tauy,
      e                lwdown, swdown,
      s                d_t, d_q, d_ts, flux_t, flux_q,dflux_s,dflux_l)
@@ -532 +538 @@
       logical debut, lafin, ok_veget
       real zlev1(klon)
+      real fder(klon), taux(klon), tauy(klon)
       real temp_air(klon), spechum(klon)
       real hum_air(klon), ccanopy(klon)
@@ -679 +686 @@
      . tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef,
      . precip_rain, precip_snow, lwdown, swnet, swdown,
+     . fder, taux, tauy,
      . albedo, snow, qsol,
      . ts, p1lay, coef1lay, psref, radsol,

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

@@ -50 +50 @@
       & tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
       & precip_rain, precip_snow, lwdown, swnet, swdown, &
+      & fder, taux, tauy, &
       & albedo, snow, qsol, &
       & tsurf, p1lay, coef1lay, ps, radsol, &
@@ -71 +72 @@
 !   dtime        pas de temps de la physique (en s)
 !   jour         jour dans l'annee en cours
+!   nexca        pas de temps couplage
 !   nisurf       index de la surface a traiter (1 = sol continental)
 !   knon         nombre de points de la surface a traiter
@@ -105 +107 @@
 !   radsol       rayonnement net aus sol (LW + SW)
 !   ocean        type d'ocean utilise (force, slab, couple)
+!   fder         derivee des flux (pour le couplage)
+!   taux, tauy   tension de vents
 !
 ! output:
@@ -143 +147 @@
   real, dimension(knon), intent(IN) :: radsol
   real, dimension(klon), intent(IN) :: zmasq
+  real, dimension(klon), intent(IN) :: fder, taux, tauy
   character (len = 6)  :: ocean
   real, dimension(knon), intent(INOUT) :: evap, snow, qsol
@@ -162 +167 @@
   real, parameter      :: calice=1.0/(5.1444e+06*0.15), tau_gl=1./86400.*5.
   real, parameter      :: calsno=1./(2.3867e+06*.15)
+  integer              :: nexca !pas de temps couplage
+  real, dimension(knon):: alb_ice
 
 #include "YOMCST.inc"
@@ -253 +260 @@
 ! Surface "ocean" appel a l'interface avec l'ocean
 !
-!    if (ocean == 'couple') then
-!      call interfoce(nisurf, ocean)
+    if (ocean == 'couple') then
+      nexca = 0
+      if (nexca == 0) then
+        abort_message='nexca = 0 dans interfoce_cpl'
+        call abort_gcm(modname,abort_message,1)
+      endif
+
+      call interfoce(itime, dtime, &
+      & klon, iim, jjm, nisurf, pctsrf, knon, knindex, rlon, rlat, &
+      & ocean, nexca, debut, lafin, &
+      & swdown, lwdown, precip_rain, precip_snow, evap, tsurf, &
+      & fder, albedo, taux, tauy, &
+      & tsurf_new, alb_new, alb_ice, pctsrf_new)
+
 !    else if (ocean == 'slab  ') then
 !      call interfoce(nisurf)
@@ -262 +281 @@
 !     &  debut, &
 !     &  tsurf_new, alb_new, z0_new, pctsrf_new)
-!    endif
+    endif
 !
     cal = 0.
@@ -556 +575 @@
 
 ! Cette routine sert d'interface entre le modele atmospherique et un 
-! coupleur avec un modele d'ocean 'complet'
+! coupleur avec un modele d'ocean 'complet' derriere
 !
 ! L. Fairhead 02/2000

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

@@ -865 +865 @@
 
 
+      SUBROUTINE sipc_init_model
+      print *, 'Attention dans oasis.F, sipc_init_model est non defini'
+      RETURN
+      END
+
+      SUBROUTINE svipc_write
+      print *, 'Attention dans oasis.F, svipc_write est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_export
+      print *, 'Attention dans oasis.F, clim_export est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_init
+      print *, 'Attention dans oasis.F, clim_init est non defini'
+      RETURN
+      END
+
+      SUBROUTINE sipc_write_model
+      print *, 'Attention dans oasis.F, sipc_write_model est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_start
+      print *, 'Attention dans oasis.F, clim_start est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_define
+      print *, 'Attention dans oasis.F, clim_define est non defini'
+      RETURN
+      END
+
+      SUBROUTINE sipc_attach
+      print *, 'Attention dans oasis.F, sipc_attach est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_import
+      print *, 'Attention dans oasis.F, clim_import est non defini'
+      RETURN
+      END
+
+      SUBROUTINE svipc_read
+      print *, 'Attention dans oasis.F, svipc_read est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_stepi
+      print *, 'Attention dans oasis.F, clim_stepi est non defini'
+      RETURN
+      END
+
+      SUBROUTINE sipc_read_model
+      print *, 'Attention dans oasis.F, sipc_read_model est non defini'
+      RETURN
+      END
+
+      SUBROUTINE svipc_close
+      print *, 'Attention dans oasis.F, svipc_close est non defini'
+      RETURN
+      END
+
+      SUBROUTINE clim_quit
+      print *, 'Attention dans oasis.F, clim_quit est non defini'
+      RETURN
+      END
+

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

@@ -1556 +1556 @@
      e            ftsol,paprs,pplay,radsol,
      e            fsnow,fqsol,fevap,falbe,
-     e            rain_fall, snow_fall, solsw, sollw,
+     e            rain_fall, snow_fall, solsw, sollw, fder,
      e            rlon, rlat, frugs,
      e            debut, lafin,