Changeset 995 for LMDZ4/trunk/libf/dyn3d


Ignore:
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

Location:
LMDZ4/trunk/libf/dyn3d
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/dyn3d/abort_gcm.F

    r524 r995  
    1919C         ierr    = severity of situation ( = 0 normal )
    2020
    21       character*20 modname
     21      character*(*) modname
    2222      integer ierr
    23       character*80 message
     23      character*(*) message
    2424
    2525      write(lunout,*) 'in abort_gcm'
  • LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

    r988 r995  
    1111      USE fonte_neige_mod
    1212      USE pbl_surface_mod
    13       USE ocean_slab_mod
    14       USE surface_data,     ONLY : ocean, ok_veget
    1513      USE phys_state_var_mod
    1614      !
     
    4644      REAL :: alblw(klon,nbsrf)
    4745      REAL :: tsoil(klon,nsoilmx,nbsrf)
    48 cIM "slab" ocean
    49       REAL :: tslab(klon), seaice(klon)
    5046      REAL :: frugs(klon,nbsrf), agesno(klon,nbsrf)
    5147      REAL :: rugmer(klon)
     
    113109      real zx_tmp_2d(iim,jjp1)
    114110
    115       REAL :: dummy(klon, nbsrf)
     111      REAL :: dummy
    116112
    117113      logical              :: ok_newmicro
     
    152148!      CALL defrun_new(99,.TRUE.,clesphy0)
    153149      CALL conf_gcm( 99, .TRUE. , clesphy0 )
    154       call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
     150      call conf_phys(ok_journe, ok_mensuel, ok_instan,                  &
    155151     &                 ok_hf,                                           &
    156152     &                 solarlong0,seuil_inversion,                      &
     
    418414     .    jjm, rlonu, rlatv , interbar )
    419415      !
    420 cIM "slab" ocean
    421       varname = 'tslab'
    422       tslab(:) = 0.0
    423       CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,tslab,0.0,
    424      .     jjm, rlonu, rlatv , interbar )
    425 c
    426       varname = 'seaice'
    427       seaice(:) = 0.0
    428       CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,seaice,0.0,
    429      .     jjm, rlonu, rlatv , interbar )
    430       !
    431416      varname = 'rugmer'
    432417      rugmer(:) = 0.0
     
    710695      q_ancien = 0.
    711696      agesno = 0.
    712 cIM "slab" ocean
    713       tslab(1:klon) = ftsol(1:klon,is_oce)
    714       seaice = 0. !kg m-2
    715 c
    716697      frugs(1:klon,is_oce) = rugmer(1:klon)
    717698      frugs(1:klon,is_ter) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
     
    729710c avec les valeurs qui vont etre ecrit dans startphy.nc
    730711c
    731       dummy(:,:) = 1.0
    732       ocean = "slab"
     712      dummy = 1.0
    733713      pbl_tke(:,:,:) = 1.e-8
    734714      zmax0(:) = 40.
     
    745725      call pbl_surface_init(qsol, fder, snsrf, qsolsrf,
    746726     $     evap, frugs, agesno, tsoil)
    747       call ocean_slab_init(dummy(1,1), tslab, seaice, dummy(:,:))
    748727
    749728      call phyredem("startphy.nc")
  • LMDZ4/trunk/libf/dyn3d/limit_netcdf.F

    r989 r995  
    11771177      ENDDO
    11781178c
     1179c Ecriture de fichier limit toujours en simple precision, NF_FLOAT
    11791180      PRINT*, 'Ecriture du fichier limit'
    11801181c
     
    12371238      epais(2) = 1
    12381239c
    1239 #ifdef NC_DOUBLE
    1240       ierr = NF_PUT_VAR1_DOUBLE (nid,id_tim,k,DBLE(k))
    1241 c
    1242       IF (newlmt ) THEN
    1243           ierr = NF_PUT_VARA_DOUBLE (nid,id_FOCE,debut,epais
    1244      $        ,pctsrf_t(1,is_oce,k))
    1245           ierr = NF_PUT_VARA_DOUBLE (nid,id_FSIC,debut,epais
    1246      $        ,pctsrf_t(1,is_sic,k))
    1247           ierr = NF_PUT_VARA_DOUBLE (nid,id_FTER,debut,epais
    1248      $        ,pctsrf_t(1,is_ter,k))
    1249           ierr = NF_PUT_VARA_DOUBLE (nid,id_FLIC,debut,epais
    1250      $        ,pctsrf_t(1,is_lic,k))
    1251       ELSE
    1252           ierr = NF_PUT_VARA_DOUBLE (nid,id_NAT,debut,epais
    1253      $        ,phy_nat(1,k))
    1254       ENDIF
    1255 c
    1256       ierr = NF_PUT_VARA_DOUBLE (nid,id_SST,debut,epais,phy_sst(1,k))
    1257       ierr = NF_PUT_VARA_DOUBLE (nid,id_BILS,debut,epais,phy_bil(1,k))
    1258       ierr = NF_PUT_VARA_DOUBLE (nid,id_ALB,debut,epais,phy_alb(1,k))
    1259       ierr = NF_PUT_VARA_DOUBLE (nid,id_RUG,debut,epais,phy_rug(1,k))
    1260 #else
    12611240      ierr = NF_PUT_VAR1_REAL (nid,id_tim,k,FLOAT(k))
    12621241      IF (newlmt ) THEN
     
    12771256      ierr = NF_PUT_VARA_REAL (nid,id_ALB,debut,epais,phy_alb(1,k))
    12781257      ierr = NF_PUT_VARA_REAL (nid,id_RUG,debut,epais,phy_rug(1,k))
    1279 #endif
    12801258c
    12811259      ENDDO
Note: See TracChangeset for help on using the changeset viewer.