Ignore:
Timestamp:
May 21, 2025, 11:54:30 AM (6 weeks ago)
Author:
Laurent Fairhead
Message:

Correction of bug for revision 5662 for the 1D

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified LMDZ6/trunk/libf/phylmd/dyn1d/scm.f90

    r5627 r5668  
    55   USE clesphys_mod_h
    66   USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin
     7   USE phyetat0_get_mod, ONLY : phyetat0_get, phyetat0_srf
    78   USE phys_state_var_mod, ONLY : phys_state_var_init, phys_state_var_end, &
    89       clwcon, detr_therm, &
     
    2425       ql_ancien, qs_ancien, qbs_ancien, cf_ancien, rvc_ancien, &
    2526       prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
    26        u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_
     27       u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_, &
     28       hice, tice, bilg_cumul
    2729
    2830
     
    231233      CHARACTER(len=80) :: abort_message
    232234      CHARACTER(len=20) :: modname = 'scm'
     235      LOGICAL :: found
    233236
    234237      print*,'VOUS ENTREZ DANS LE 1D FORMAT STANDARD'
     
    613616
    614617!-----------------------------------------------------------------------
    615         call pbl_surface_init(fder, snsrf, qsurfsrf, tsoil)
     618  !GG
     619  ! Sea ice
     620  !IF (iflag_seaice == 2) THEN
     621
     622  found=phyetat0_get(hice,"hice","Ice thickness",0.)
     623  IF (.NOT. found) THEN
     624       PRINT*, "phyetat0: Le champ <hice> est absent"
     625       PRINT*, "Initialisation a hice=1m "
     626       hice(:)=1.0
     627  END IF
     628  found=phyetat0_get(tice,"tice","Sea Ice temperature",0.)
     629  IF (.NOT. found) THEN
     630       PRINT*, "phyetat0: Le champ <tice> est absent"
     631       PRINT*, "Initialisation a tsol_sic"
     632       tice(:)=ftsol(:,is_sic)
     633  END IF
     634  found=phyetat0_get(bilg_cumul,"bilg_cumul","Flux conductivite + transmit sea-ice",0.)
     635  IF (.NOT. found) THEN
     636       PRINT*, "phyetat0: Le champ <bilg_cumul> est absent"
     637       PRINT*, "Initialisation a zero"
     638       bilg_cumul(:)=0.0
     639  END IF
     640
     641  !END IF
     642  !GG
     643        call pbl_surface_init(fder, snsrf, qsurfsrf, tsoil, hice, tice, bilg_cumul)
    616644
    617645!------------------ prepare limit conditions for limit.nc -----------------
Note: See TracChangeset for help on using the changeset viewer.