Changeset 786 for LMDZ4/trunk/libf/dyn3d

Timestamp:

Jun 19, 2007, 12:12:13 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Changement dans le restart suite a la réorganisation de l'interface avec
les surfaces JG
LF

File:

• LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F (modified) (4 diffs)

LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

@@ -9 +9 @@
       USE ioipsl
       USE dimphy
+      USE fonte_neige_mod
+      USE pbl_surface_mod
+      USE ocean_slab_mod
+      USE surface_data
       !
       IMPLICIT NONE
@@ -24 +28 @@
 #include "comconst.h"
 #include "indicesol.h"
-c#include "dimphy.h"
 #include "dimsoil.h"
 #include "temps.h"
@@ -115 +118 @@
       integer :: isst(klon-2)
       real zx_tmp_2d(iim,jjp1)
+
+
+      REAL :: dummy(klon, nbsrf)
       !
       !   Constantes 
@@ -674 +680 @@
       run_off_lic_0 = 0.0 
 
-cIM   call phyredem("startphy.nc",phystep,radpas, co2_ppm, solaire,
-    
-      call phyredem("startphy.nc",phystep,radpas,
-     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, tslab, seaice, 
-     $    qsolsrf, qsol,
-     $    snsrf, 
-     $    albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder,
-     $    radsol, frugs,  agesno, 
+c
+c Avant l'appel a phyredem, on initialize les modules de surface
+c avec les valeurs qui vont etre ecrit dans startphy.nc
+c
+      dummy(:,:) = 1.0
+      ocean = "slab"
+      call fonte_neige_init(run_off_lic_0)
+      call pbl_surface_init(qsol, fder, snsrf, qsolsrf,
+     $     evap, frugs, agesno, tsoil)
+      call ocean_slab_init(dummy(1,1), tslab, seaice, dummy(:,:))
+
+      call phyredem("startphy.nc",phystep,radpas,ocean, 
+     $    latfi, lonfi, pctsrf, tsolsrf,
+     $    albe, alblw, rain_fall, snow_fall, solsw, sollw,
+     $    radsol, 
      $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel, 
-     $    t_ancien, q_ancien, rnebcon, ratqs, clwcon,
-     $    run_off_lic_0)
+     $    t_ancien, q_ancien, rnebcon, ratqs, clwcon )
+
+
 
 C     Sortie Visu pour les champs dynamiques