Changeset 782 for LMDZ4/trunk/libf/phylmd/phyetat0.F

Timestamp:

Jun 11, 2007, 4:50:43 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Adaptation du code a la nouvelle interface avec les surface de Josefine
LF

File:

• LMDZ4/trunk/libf/phylmd/phyetat0.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -5 +5 @@
 c
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
-     .            rlat_p,rlon_p, pctsrf_p, tsol_p,tsoil_p,
-cIM "slab" ocean
-     .           ocean, tslab_p,seaice_p,
-     .           qsurf_p,qsol_p,snow_p,albe_p, alblw_p, evap_p,
+     .           rlat_p, rlon_p, pctsrf_p, tsol_p,
+     .           ocean_in, ok_veget_in,
+     .           albe_p, alblw_p,
      .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
-     .           fder_p,radsol_p,frugs_p,agesno_p,clesphy0,
+     .           radsol_p,clesphy0,
      .           zmea_p,zstd_p,zsig_p,zgam_p,zthe_p,zpic_p,zval_p,
      .           rugsrel_p,tabcntr0,
      .           t_ancien_p,q_ancien_p,ancien_ok_p, rnebcon_p, ratqs_p,
-     .           clwcon_p,run_off_lic_0_p)
+     .           clwcon_p)
+
       USE dimphy
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_para
       USE iophy
+      USE ocean_slab_mod,   ONLY : ocean_slab_init
+      USE ocean_cpl_mod,    ONLY : ocean_cpl_init
+      USE ocean_forced_mod, ONLY : ocean_forced_init
+      USE fonte_neige_mod,  ONLY : fonte_neige_init
+      USE pbl_surface_mod,  ONLY : pbl_surface_init
+      USE surface_data,     ONLY : ocean, ok_veget
+
       IMPLICIT none
 c======================================================================
@@ -25 +32 @@
 c======================================================================
 #include "dimensions.h"
-cym#include "dimphy.h"
 #include "netcdf.inc"
 #include "indicesol.h"
@@ -111 +117 @@
       real ratqs(klon_glo,klev)
 
-      CHARACTER*6 ocean
+      CHARACTER*6 ocean_in
+      LOGICAL ok_veget_in
 
       INTEGER        longcles
@@ -515 +522 @@
 c Lecture de tslab (pour slab ocean seulement):      
 c
-      IF (ocean .eq. 'slab  ') then
+      IF (ocean_in .eq. 'slab  ') then
         ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
         IF (ierr.NE.NF_NOERR) THEN
@@ -1543 +1550 @@
       call Scatter( zmasq_glo,zmasq)
 
+c
+c Initilalize variables in module surface_data
+c
+      ok_veget = ok_veget_in
+      ocean    = ocean_in
+c
+c Initialize module pbl_surface_mod 
+c
+      CALL pbl_surface_init(qsol_p, fder_p, snow_p, qsurf_p,
+     $     evap_p, frugs_p, agesno_p, tsoil_p)
+
+c Initialize ocean module according to ocean type
+      IF ( ocean == 'slab' ) THEN
+c        initilalize module ocean_slab_init
+         CALL ocean_slab_init(dtime, tslab_p, seaice_p, pctsrf_p)
+      ELSEIF ( ocean == 'couple' ) THEN
+c        initilalize module ocean_cpl_init
+         CALL ocean_cpl_init(dtime, rlon_p, rlat_p)
+      ELSE
+c        initilalize module ocean_forced_init
+         CALL ocean_forced_init
+      ENDIF
+c
+c Initilialize module fonte_neige_mod      
+c
+      CALL fonte_neige_init(run_off_lic_0_p)
+
+
       RETURN
       END