Changeset 995 for LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

Timestamp:

Sep 9, 2008, 3:17:08 PM (16 years ago)

Author:

lsce

Message:

• Modifications liées au calcul des nouveau sous-fractions
• Nettoyage de ocean slab : il reste uniquement la version avec glace de mer forcé
• Nouveaux variables pour distiguer la version et type d'ocean : type_ocean=force/slab/couple, version_ocean=opa8/nemo pour couplé ou version_ocean=sicOBS pour slab

JG

File:

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

LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

@@ -11 +11 @@
       USE fonte_neige_mod
       USE pbl_surface_mod
-      USE ocean_slab_mod
-      USE surface_data,     ONLY : ocean, ok_veget
       USE phys_state_var_mod
       !
@@ -46 +44 @@
       REAL :: alblw(klon,nbsrf)
       REAL :: tsoil(klon,nsoilmx,nbsrf) 
-cIM "slab" ocean
-      REAL :: tslab(klon), seaice(klon)
       REAL :: frugs(klon,nbsrf), agesno(klon,nbsrf)
       REAL :: rugmer(klon)
@@ -113 +109 @@
       real zx_tmp_2d(iim,jjp1)
 
-      REAL :: dummy(klon, nbsrf)
+      REAL :: dummy
 
       logical              :: ok_newmicro
@@ -152 +148 @@
 !      CALL defrun_new(99,.TRUE.,clesphy0)
       CALL conf_gcm( 99, .TRUE. , clesphy0 )
-      call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
+      call conf_phys(ok_journe, ok_mensuel, ok_instan,                  &
      &                 ok_hf,                                           &
      &                 solarlong0,seuil_inversion,                      &
@@ -418 +414 @@
      .    jjm, rlonu, rlatv , interbar )
       !
-cIM "slab" ocean
-      varname = 'tslab'
-      tslab(:) = 0.0
-      CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,tslab,0.0,
-     .     jjm, rlonu, rlatv , interbar )
-c
-      varname = 'seaice'
-      seaice(:) = 0.0
-      CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,seaice,0.0,
-     .     jjm, rlonu, rlatv , interbar )
-      !
       varname = 'rugmer'
       rugmer(:) = 0.0
@@ -710 +695 @@
       q_ancien = 0.
       agesno = 0.
-cIM "slab" ocean
-      tslab(1:klon) = ftsol(1:klon,is_oce)
-      seaice = 0. !kg m-2
-c
       frugs(1:klon,is_oce) = rugmer(1:klon)
       frugs(1:klon,is_ter) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
@@ -729 +710 @@
 c avec les valeurs qui vont etre ecrit dans startphy.nc
 c
-      dummy(:,:) = 1.0
-      ocean = "slab"
+      dummy = 1.0
       pbl_tke(:,:,:) = 1.e-8 
       zmax0(:) = 40.
@@ -745 +725 @@
       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")