Changeset 782 for LMDZ4/trunk/libf/phylmd/phyredem.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/phyredem.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/phyredem.F

@@ -3 +3 @@
 !
 c
-      SUBROUTINE phyredem (fichnom,dtime,radpas,
-     .           rlat_p,rlon_p, pctsrf_p,tsol_p,tsoil_p,
-cIM "slab" ocean
-     .           tslab,seaice,
-     .           qsurf_p,qsol_p,snow_p,albedo_p, alblw_p, evap_p,
-     .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,fder_p,
-     .           radsol_p,frugs_p,agesno_p,zmea_p,zstd_p,zsig_p,
+      SUBROUTINE phyredem (fichnom,dtime,radpas,ocean,
+     .           rlat_p,rlon_p, pctsrf_p,tsol_p,
+     .           albedo_p, alblw_p,
+     .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
+     .           radsol_p,zmea_p,zstd_p,zsig_p,
      .           zgam_p,zthe_p,zpic_p,zval_p,rugsrel_p,
-     .           t_ancien_p, q_ancien_p, rnebcon_p, ratqs_p, clwcon_p,
-     .           run_off_lic_0_p)
+     .           t_ancien_p, q_ancien_p, rnebcon_p, ratqs_p, clwcon_p)
+
       USE dimphy
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_para
+      USE ocean_slab_mod,   ONLY : ocean_slab_final
+      USE fonte_neige_mod,  ONLY : fonte_neige_final
+      USE pbl_surface_mod,  ONLY : pbl_surface_final
+
       IMPLICIT none
 c======================================================================
@@ -21 +23 @@
 c Objet: Ecriture de l'etat de redemarrage pour la physique
 c======================================================================
-cym#include "dimensions.h"
-cym#include "dimphy.h"
 #include "netcdf.inc"
 #include "indicesol.h"
@@ -36 +36 @@
       REAL tsol_p(klon,nbsrf)
       REAL tsoil_p(klon,nsoilmx,nbsrf)
+      CHARACTER*6 ocean
 cIM "slab" ocean
       REAL tslab_p(klon), seaice_p(klon)
@@ -110 +111 @@
       CHARACTER*7 str7
       CHARACTER*2 str2
-c
+
+c======================================================================
+c 
+c Get variables which will be written to restart file from module 
+c pbl_surface_mod
+      CALL pbl_surface_final(qsol_p, fder_p, snow_p, qsurf_p, 
+     $     evap_p, frugs_p, agesno_p, tsoil_p)
+
+c Get a variable calculated in module fonte_neige_mod
+      CALL fonte_neige_final(run_off_lic_0_p)
+
+c If slab ocean then get 2 varaibles from module ocean_slab_mod
+      IF ( ocean == 'slab' ) THEN
+         CALL ocean_slab_final(tslab_p, seaice_p)
+      ELSE 
+         tslab_p(:)  = 0.0
+         seaice_p(:) = 0.0
+      ENDIF      
+
+c======================================================================
+
       call Gather( rlat_p,rlat)
       call Gather( rlon_p,rlon)