Ignore:
Timestamp:
Oct 22, 2020, 2:50:18 PM (4 years ago)
Author:
evignon
Message:

Premiere comission Etienne: changements pour le 1D (forcage en Ts au dessus des continents) et inclusion drag arbres dans yamada4_num=6

Location:
LMDZ6/trunk/libf/phylmd/dyn1d
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

    r3686 r3780  
    233233       CALL getin('ok_flux_surf',ok_flux_surf)
    234234
     235!Config  Key  = ok_forc_tsurf
     236!Config  Desc = forcage ou non par la Ts
     237!Config  Def  = false
     238!Config  Help = forcage ou non par la Ts
     239       ok_forc_tsurf=.false.
     240       CALL getin('ok_forc_tsurf',ok_forc_tsurf)
     241
    235242!Config  Key  = ok_prescr_ust
    236243!Config  Desc = ustar impose ou non
    … …  
    239246       ok_prescr_ust = .false.
    240247       CALL getin('ok_prescr_ust',ok_prescr_ust)
     248
     249
     250!Config  Key  = ok_prescr_beta
     251!Config  Desc = betaevap impose ou non
     252!Config  Def  = false
     253!Config  Help = betaevap impose ou non
     254       ok_prescr_beta = .false.
     255       CALL getin('ok_prescr_beta',ok_prescr_beta)
    241256
    242257!Config  Key  = ok_old_disvert
    … …  
    280295!Config  Desc = surface temperature
    281296!Config  Def  = 290.
    282 !Config  Help = not used if type_ts_forcing=1 in lmdz1d.F
     297!Config  Help = surface temperature
    283298       tsurf = 290.
    284299       CALL getin('tsurf',tsurf)
    … …  
    297312       zsurf = 0.
    298313       CALL getin('zsurf',zsurf)
     314! EV pour accord avec format standard       
     315       CALL getin('zorog',zsurf)
     316
    299317
    300318!Config  Key  = rugos
    … …  
    359377       qsolinp = 1.
    360378       CALL getin('qsolinp',qsolinp)
     379
     380
     381
     382!Config  Key  = betaevap
     383!Config  Desc = beta for actual evaporation when prescribed
     384!Config  Def  = 1.0
     385!Config  Help =
     386       betaevap = 1.
     387       CALL getin('betaevap',betaevap)     
    361388
    362389!Config  Key  = zpicinp
    … …  
    520547       CALL getin('forc_ustar',forc_ustar)
    521548       IF (forc_ustar .EQ. 1) ok_prescr_ust=.true.
     549
    522550
    523551!Config  Key  = nudging_u
    … …  
    12481276      END
    12491277
    1250 !======================================================================
    1251        SUBROUTINE read_tsurf1d(knon,sst_out)
    1252 
    1253 ! This subroutine specifies the surface temperature to be used in 1D simulations
    1254 
    1255       USE dimphy, ONLY : klon
    1256 
    1257       INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
    1258       REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
    1259 
    1260        INTEGER :: i
    1261 ! COMMON defined in lmdz1d.F:
    1262        real ts_cur
    1263        common /sst_forcing/ts_cur
    1264 
    1265        DO i = 1, knon
    1266         sst_out(i) = ts_cur
    1267        ENDDO
    1268 
    1269       END SUBROUTINE read_tsurf1d
    1270 
     1278!!======================================================================
     1279!       SUBROUTINE read_tsurf1d(knon,sst_out)
     1280!
     1281!! This subroutine specifies the surface temperature to be used in 1D simulations
     1282!
     1283!      USE dimphy, ONLY : klon
     1284!
     1285!      INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
     1286!      REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
     1287!
     1288!       INTEGER :: i
     1289!! COMMON defined in lmdz1d.F:
     1290!       real ts_cur
     1291!       common /sst_forcing/ts_cur
     1292
     1293!       DO i = 1, knon
     1294!        sst_out(i) = ts_cur
     1295!       ENDDO
     1296!
     1297!      END SUBROUTINE read_tsurf1d
     1298!
    12711299!===============================================================
    12721300      subroutine advect_vert(llm,w,dt,q,plev)

  • LMDZ6/trunk/libf/phylmd/dyn1d/1D_decl_cases.h

    r3686 r3780  
    3434        real w_mod(llm), t_mod(llm),q_mod(llm)
    3535        real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm),ug_mod(llm),vg_mod(llm)
    36         real temp_nudg_mod(llm),qv_nudg_mod(llm),u_nudg_mod(llm),v_nudg_mod(llm)
     36            real temp_nudg_mod(llm),qv_nudg_mod(llm),u_nudg_mod(llm),v_nudg_mod(llm)
    3737        real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm)
    3838        real th_mod(llm)
    3939
    40         real ts_cur
    41         common /sst_forcing/ts_cur ! also in read_tsurf1d.F
     40! EV comment these lines
     41!        real ts_cur
     42!        common /sst_forcing/ts_cur ! also in read_tsurf1d.F
    4243!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    4344! Declarations specifiques au cas RICO

  • LMDZ6/trunk/libf/phylmd/dyn1d/1D_interp_cases.h

    r3686 r3780  
    11
    2          print*,'FORCING CASE forcing_case2'
     2        print*,'FORCING CASE forcing_case2'
    33!       print*,                                                             &
    44!    & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/pdt_cas=',     &
    … …  
    2828     &       ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas                           &
    2929     &       ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas)
    30 
    31              ts_cur = ts_prof_cas
     30! EV tg instead of ts_cur
     31             tg = ts_prof_cas
    3232!            psurf=plev_prof_cas(1)
    3333             psurf=ps_prof_cas

  • LMDZ6/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

    r3686 r3780  
    7070
    7171! initial and boundary conditions :
    72 !      tsurf = ts_prof_cas
     72!     tsurf = ts_prof_cas
    7373      psurf = ps_prof_cas
    74       ts_cur = ts_prof_cas
     74      !EV tg instead of ts_cur
     75      tg = ts_prof_cas
     76      print*, 'tg=', tg
     77
    7578      do l = 1, llm
    7679       temp(l) = t_mod_cas(l)
    … …  
    108111       IF (ok_prescr_ust) THEN
    109112       ust=ustar_prof_cas
    110        print *,'ust=',ust
    111113       ENDIF
    112114
     115
    113116      endif !forcing_SCM

  • LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read_std.F90

    r3688 r3780  
    9292      REAL, ALLOCATABLE :: time_val(:)
    9393
    94       print*,'ON EST VRAIMENT LA'
     94      print*,'ON EST VRAIMENT DASN MOD_1D_CASES_READ_STD'
    9595      fich_cas='cas.nc'
    9696      print*,'fich_cas ',fich_cas
    … …  
    924924!      enddo
    925925
     926!       print*, 'plev_prof_cas', plev_prof_cas
     927!       print*, 'play', play
    926928       do l = 1, llm
    927929
    … …  
    951953
    952954         frac = (plev_prof_cas(k2)-play(l))/(plev_prof_cas(k2)-plev_prof_cas(k1))
     955         
    953956         t_mod_cas(l)= t_prof_cas(k2) - frac*(t_prof_cas(k2)-t_prof_cas(k1))
    954957         theta_mod_cas(l)= th_prof_cas(k2) - frac*(th_prof_cas(k2)-th_prof_cas(k1))
    … …  
    10751078       enddo ! l
    10761079
     1080
    10771081          return
    10781082          end SUBROUTINE interp2_case_vertical_std

  • LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_decl_cases.h

    r3593 r3780  
    3737        real th_mod(llm)
    3838
    39         real ts_cur
    40         common /sst_forcing/ts_cur ! also in read_tsurf1d.F
     39        !real ts_cur
     40        !common /sst_forcing/ts_cur ! also in read_tsurf1d.F
    4141!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    4242! Declarations specifiques au cas RICO

  • LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_interp_cases.h

    r3593 r3780  
    6262     &             ,ts_prof,plev_prof,t_prof,q_prof,u_prof,v_prof,w_prof    &
    6363     &             ,ht_prof,vt_prof,hq_prof,vq_prof)
    64 
    65         if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
     64! EV: tg instead of ts_cur
     65        if (type_ts_forcing.eq.1) tg = ts_prof !
    6666
    6767! vertical interpolation:
    … …  
    113113!     print *,'llm l omega_profd',llm,l,omega_profd(l)
    114114!     enddo
    115 
    116         if (type_ts_forcing.eq.1) ts_cur = tg_prof ! SST used in read_tsurf1d
     115! EV tg instead of ts_cur
     116        if (type_ts_forcing.eq.1) tg = tg_prof ! SST used
    117117
    118118! vertical interpolation:
    … …  
    206206     &             ,ug_gabls4,vg_gabls4,ht_gabls4,hq_gabls4,tg_gabls4                            &
    207207     &             ,ug_profg,vg_profg,ht_profg,hq_profg,tg_profg)
    208 
    209         if (type_ts_forcing.eq.1) ts_cur = tg_prof ! SST used in read_tsurf1d
     208!EV tg instead of ts_cur
     209        if (type_ts_forcing.eq.1) tg = tg_prof ! SST used
    210210
    211211! vertical interpolation:
    … …  
    499499     &             ,nlev_sandu                                              &
    500500     &             ,ts_sandu,ts_prof)
    501 
    502         if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
     501! EV tg instead of ts_cur
     502        if (type_ts_forcing.eq.1) tg = ts_prof ! SST used in read_tsurf1d
    503503
    504504! vertical interpolation:
    … …  
    582582     &             ,ufa_astex,vfa_astex,div_prof,ts_prof,ug_prof,vg_prof    &
    583583     &             ,ufa_prof,vfa_prof)
    584 
    585         if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
    586 
     584! EV tg instead of ts_cur
     585        if (type_ts_forcing.eq.1) tg = ts_prof ! SST used
    587586! vertical interpolation:
    588587      CALL interp_astex_vertical(play,nlev_astex,plev_profa                 &
    … …  
    675674     &       ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas,lat_prof_cas               &
    676675     &       ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas)
    677 
    678              ts_cur = ts_prof_cas
     676! EV tg instead of ts_cur
     677
     678             tg = ts_prof_cas
    679679             psurf=plev_prof_cas(1)
    680680
    … …  
    850850     &       ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas                           &
    851851     &       ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas)
    852 
    853              ts_cur = ts_prof_cas
     852! EV tg instead of ts_cur
     853
     854             tg = ts_prof_cas
    854855!            psurf=plev_prof_cas(1)
    855856             psurf=ps_prof_cas

  • LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_read_forc_cases.h

    r3679 r3780  
    875875
    876876! initial and boundary conditions :
    877 !      tsurf = ts_prof_cas
    878       ts_cur = ts_prof_cas
     877!     tsurf = ts_prof_cas
     878! EV tg instead of ts_cur
     879      tg= ts_prof_cas
    879880      psurf=plev_prof_cas(1)
    880881      write(*,*) 'SST initiale: ',tsurf
    … …  
    965966! initial and boundary conditions :
    966967!      tsurf = ts_prof_cas
    967       ts_cur = ts_prof_cas
     968! EV tg instead of ts_cur
     969      tg = ts_prof_cas
    968970      psurf=plev_prof_cas(1)
    969971      write(*,*) 'SST initiale: ',tsurf
    … …  
    10631065! initial and boundary conditions :
    10641066!      tsurf = ts_prof_cas
    1065       ts_cur = ts_prof_cas
     1067! EV tg instead of ts_cur
     1068
     1069      tg = ts_prof_cas
    10661070      psurf=plev_prof_cas(1)
    10671071      write(*,*) 'SST initiale: ',tsurf

  • LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

    r3594 r3780  
    728728
    729729!Al1 pour SST forced, appell?? depuis ocean_forced_noice
    730       ts_cur = tsurf ! SST used in read_tsurf1d
     730! EV tg instead of ts_cur
     731
     732      tg = tsurf ! SST used in read_tsurf1d
    731733!=====================================================================
    732734! Initialisation de la physique :
    … …  
    791793
    792794        fder=0.
     795        print *, 'snsrf', snsrf
    793796        snsrf(1,:)=snowmass ! masse de neige des sous surface
    794797        qsurfsrf(1,:)=qsurf ! humidite de l'air des sous surface
    … …  
    841844     end if
    842845
    843 
    844846        print*,'nat_surf,pctsrf(1,is_oce),pctsrf(1,is_ter)',nat_surf         &
    845847     &        ,pctsrf(1,is_oce),pctsrf(1,is_ter)
    … …  
    851853              ! 6 albedo, mais on peut quand meme tourner avec
    852854              ! moins. Seules les 2 ou 4 premiers seront lus
     855              print*, 'les albedos sont sont', albedo, falb_dir
    853856        falb_dir=albedo
    854857        falb_dif=albedo
     858             print*, falb_dir
    855859        rugoro=rugos
    856860        t_ancien(1,:)=temp(:)
    … …  
    913917        v_ancien(1,:)=v(:)
    914918
    915 u10m=0.
    916 v10m=0.
    917 ale_wake=0.
    918 ale_bl_stat=0.
     919        u10m=0.
     920        v10m=0.
     921        ale_wake=0.
     922        ale_bl_stat=0.
    919923 
    920924!------------------------------------------------------------------------

  • LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

    r3693 r3780  
    7575      real :: zcufi    = 1.
    7676      real :: zcvfi    = 1.
    77 
    78 !-      real :: nat_surf
    79 !-      logical :: ok_flux_surf
    80 !-      real :: fsens
    81 !-      real :: flat
    82 !-      real :: tsurf
    83 !-      real :: rugos
    84 !-      real :: qsol(1:2)
    85 !-      real :: qsurf
    86 !-      real :: psurf
    87 !-      real :: zsurf
    88 !-      real :: albedo
    89 !-
    90 !-      real :: time     = 0.
    91 !-      real :: time_ini
    92 !-      real :: xlat
    93 !-      real :: xlon
    94 !-      real :: wtsurf
    95 !-      real :: wqsurf
    96 !-      real :: restart_runoff
    97 !-      real :: xagesno
    98 !-      real :: qsolinp
    99 !-      real :: zpicinp
    100 !-
    10177      real :: fnday
    10278      real :: day, daytime
    … …  
    141117        logical :: forcing_case2   = .false.
    142118        logical :: forcing_SCM   = .false.
    143         integer :: type_ts_forcing ! 0 = SST constant; 1 = SST read from a file
    144 !                                                            (cf read_tsurf1d.F)
    145119
    146120!flag forcings
    … …  
    148122        logical :: nudge_thermo=.false.
    149123        logical :: cptadvw=.true.
     124
     125
    150126!=====================================================================
    151127! DECLARATIONS FOR EACH CASE
    … …  
    248224!
    249225      integer :: it_end ! iteration number of the last call
    250 !Al1
     226!Al1,plev,play,phi,phis,presnivs,
    251227      integer ecrit_slab_oc !1=ecrit,-1=lit,0=no file
    252228      data ecrit_slab_oc/-1/
    … …  
    278254      d_v_age(:)=0.
    279255     
     256     
    280257! Initialization of Common turb_forcing
    281258       dtime_frcg = 0.
    … …  
    290267! OPTIONS OF THE 1D SIMULATION (lmdz1d.def => unicol.def)
    291268!---------------------------------------------------------------------
    292 !Al1
    293269        call conf_unicol
    294270!Al1 moves this gcssold var from common fcg_gcssold to
    … …  
    296272! --------------------------------------------------------------------
    297273        close(1)
    298 !Al1
    299274        write(*,*) 'lmdz1d.def lu => unicol.def'
    300275
    … …  
    302277       year_ini_cas=1997
    303278       ! It is possible that those parameters are run twice.
    304 
    305279       ! A REVOIR : LIRE PEUT ETRE AN MOIS JOUR DIRECETEMENT
     280
     281
    306282       call getin('anneeref',year_ini_cas)
    307283       call getin('dayref',day_deb)
    … …  
    309285       call getin('time_ini',heure_ini_cas)
    310286
    311         type_ts_forcing = 0
    312         IF (nat_surf==0) type_ts_forcing=1 ! SST forcee sur OCEAN
    313         print*,'NATURE DE LA SURFACE ',nat_surf
     287      print*,'NATURE DE LA SURFACE ',nat_surf
    314288!
    315289! Initialization of the logical switch for nudging
     290
    316291     jcode = iflag_nudge
    317292     do i = 1,nudge_max
    … …  
    319294       jcode = jcode/10
    320295     enddo
    321 !---------------------------------------------------------------------
     296!-----------------------------------------------------------------------
    322297!  Definition of the run
    323 !---------------------------------------------------------------------
     298!-----------------------------------------------------------------------
    324299
    325300      call conf_gcm( 99, .TRUE. )
    … …  
    343318      allocate( phy_flic(year_len)) ! Fraction de glace
    344319      phy_flic(:)=0.0
     320
     321
    345322!-----------------------------------------------------------------------
    346323!   Choix du calendrier
    … …  
    373350!      Le numero du jour est dans "day". L heure est traitee separement.
    374351!      La date complete est dans "daytime" (l'unite est le jour).
     352
     353
    375354      if (nday>0) then
    376355         fnday=nday
    … …  
    409388! Initialization of dimensions, geometry and initial state
    410389!---------------------------------------------------------------------
    411 !      call init_phys_lmdz(1,1,llm,1,(/1/)) ! job now done via iniphysiq
     390!     call init_phys_lmdz(1,1,llm,1,(/1/)) ! job now done via iniphysiq
    412391!     but we still need to initialize dimphy module (klon,klev,etc.)  here.
    413392      call init_dimphy1D(1,llm)
    … …  
    446425!!! Feedback forcing values for Gateaux differentiation (al1)
    447426!!!=====================================================================
    448 !!! Surface Planck forcing bracketing call radiation
    449 !!      surf_Planck = 0.
    450 !!      surf_Conv   = 0.
    451 !!      write(*,*) 'Gateaux-dif Planck,Conv:',surf_Planck,surf_Conv
    452 !!! a mettre dans le lmdz1d.def ou autre
    453 !!
    454427!!
    455428      qsol = qsolinp
    … …  
    469442      ENDIF
    470443      print*,'Flux sol ',fsens,flat
    471 !!      ok_flux_surf=.false.
    472 !!      fsens=-wtsurf*rcpd*rho(1)
    473 !!      flat=-wqsurf*rlvtt*rho(1)
    474 !!!!
    475444
    476445! Vertical discretization and pressure levels at half and mid levels:
    … …  
    496465      plev =ap+bp*psurf
    497466      play = 0.5*(plev(1:llm)+plev(2:llm+1))
    498       zlay=-rd*300.*log(play/psurf)/rg ! moved after reading profiles
     467      zlay=-rd*300.*log(play/psurf)/rg ! moved after reading profiles.
    499468
    500469      IF (forcing_type .eq. 59) THEN
    … …  
    527496      print*,'mxcalc=',mxcalc
    528497!     print*,'zlay=',zlay(mxcalc)
    529       print*,'play=',play(mxcalc)
    530 
    531 !Al1 pour SST forced, appell?? depuis ocean_forced_noice
    532       ts_cur = tsurf ! SST used in read_tsurf1d
     498!      print*,'play=',play(mxcalc)
     499
     500!! When surface temperature is forced
     501      tg= tsurf ! surface T used in read_tsurf1d
     502
     503
    533504!=====================================================================
    534505! Initialisation de la physique :
    … …  
    546517! airefi,zcufi,zcvfi initialises au debut de ce programme
    547518! rday,ra,rg,rd,rcpd declares dans YOMCST.h et calcules dans suphel.F
     519
     520
    548521      day_step = float(nsplit_phys)*day_step/float(iphysiq)
    549522      write (*,*) 'Time step divided by nsplit_phys (=',nsplit_phys,')'
    … …  
    563536     ! e.g. for cell boundaries, which are meaningless in 1D; so pad these
    564537     ! with '0.' when necessary
     538
    565539      call iniphysiq(iim,jjm,llm, &
    566540           1,comm_lmdz, &
    … …  
    653627              ! 6 albedo, mais on peut quand meme tourner avec
    654628              ! moins. Seules les 2 ou 4 premiers seront lus
     629              print*, 'les albedos sont', albedo
    655630        falb_dir=albedo
    656631        falb_dif=albedo
    … …  
    664639        prw_ancien = 0.
    665640!jyg<
    666 !!        pbl_tke(:,:,:)=1.e-8
     641!!      pbl_tke(:,:,:)=1.e-8
    667642        pbl_tke(:,:,:)=0.
    668         pbl_tke(:,2,:)=1.e-2
     643! EV: pourquoi????
     644!        pbl_tke(:,2,:)=1.e-2
    669645        PRINT *, ' pbl_tke dans lmdz1d '
    670646        if (prt_level .ge. 5) then
    … …  
    675651
    676652!>jyg
    677 
    678653        rain_fall=0.
    679654        snow_fall=0.
    … …  
    715690        v_ancien(1,:)=v(:)
    716691 
    717 u10m=0.
    718 v10m=0.
    719 ale_wake=0.
    720 ale_bl_stat=0.
     692        u10m=0.
     693        v10m=0.
     694        ale_wake=0.
     695        ale_bl_stat=0.
    721696
    722697!------------------------------------------------------------------------
    … …  
    738713! to be set at some arbitratry convenient values.
    739714!------------------------------------------------------------------------
    740 !Al1 =============== restart option ==========================
     715!Al1 =============== restart option ======================================
    741716        if (.not.restart) then
    742717          iflag_pbl = 5
    … …  
    803778       print*,'plev,play,phi,phis,presnivs,u,v,temp,q,omega2'
    804779       print*,'temp(1),q(1,1),u(1),v(1),plev(1),phis :'
    805        print*,temp(1),q(1,1),u(1),v(1),plev(1),phis
     780       print*,temp(1),q(1,1),u(1),v(1),plev(1),phis(1)
    806781! raz for safety
    807782       do l=1,llm
    … …  
    809784       enddo
    810785      endif
    811 !Al1 ================  end restart =================================
     786!======================  end restart =================================
    812787      IF (ecrit_slab_oc.eq.1) then
    813788         open(97,file='div_slab.dat',STATUS='UNKNOWN')
    … …  
    820795       CALL iophys_ini
    821796#endif
     797
     798!=====================================================================
    822799! START OF THE TEMPORAL LOOP :
    823800!=====================================================================
    824801           
    825802      it_end = nint(fnday*day_step)
    826 !test JLD     it_end = 10
    827803      do while(it.le.it_end)
    828804
    … …  
    832808         print*,'PAS DE TEMPS ',timestep
    833809       endif
    834 !Al1 demande de restartphy.nc
    835810       if (it.eq.it_end) lastcall=.True.
    836811
    … …  
    844819!  Geopotential :
    845820!---------------------------------------------------------------------
    846 
     821!        phis(1)=zsurf*RG
     822!        phi(1)=phis(1)+RD*temp(1)*(plev(1)-play(1))/(.5*(plev(1)+play(1)))
    847823        phi(1)=RD*temp(1)*(plev(1)-play(1))/(.5*(plev(1)+play(1)))
     824
    848825        do l = 1, llm-1
    849826          phi(l+1)=phi(l)+RD*(temp(l)+temp(l+1))*                           &
    850827     &    (play(l)-play(l+1))/(play(l)+play(l+1))
    851828        enddo
     829
    852830
    853831!---------------------------------------------------------------------
    … …  
    950928       sfdt = sin(0.5*fcoriolis*timestep)
    951929       cfdt = cos(0.5*fcoriolis*timestep)
    952 !       print *,'fcoriolis,sfdt,cfdt,timestep',fcoriolis,sfdt,cfdt,timestep
    953 !
     930
    954931        d_u_age(1:mxcalc)= -2.*sfdt/timestep*                                &
    955932     &          (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) -                          &
    … …  
    10301007        temp(1:mxcalc)=temp(1:mxcalc)+timestep*(                            &
    10311008     &              dt_phys(1:mxcalc)                                       &
    1032      &             +d_t_adv(1:mxcalc)                                      &
    1033      &             +d_t_nudge(1:mxcalc)                                      &
     1009     &             +d_t_adv(1:mxcalc)                                       &
     1010     &             +d_t_nudge(1:mxcalc)                                     &
    10341011     &             +dt_cooling(1:mxcalc))  ! Taux de chauffage ou refroid.
    10351012
    10361013
    1037 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     1014!=======================================================================
    10381015!! CONSERVE EN ATTENDANT QUE LE CAS EN QUESTION FONCTIONNE EN STD !!
    1039 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    1040 !     endif  ! forcing_sandu or forcing_astex
    1041 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     1016!=======================================================================
    10421017
    10431018        teta=temp*(pzero/play)**rkappa
    1044 !
     1019
    10451020!---------------------------------------------------------------------
    10461021!   Nudge soil temperature if requested
    … …  
    10801055
    10811056!  incremente day time
    1082 !        print*,'daytime bef',daytime,1./day_step
    10831057        daytime = daytime+1./day_step
    1084 !Al1dbg
    10851058        day = int(daytime+0.1/day_step)
    10861059!        time = max(daytime-day,0.0)
    … …  
    10881061!cc        time = real(mod(it,day_step))/day_step
    10891062        time = time_ini/24.+real(mod(it,day_step))/day_step
    1090 !        print*,'daytime nxt time',daytime,time
    10911063        it=it+1
    10921064
    10931065      enddo
    10941066
    1095 !Al1
    10961067      if (ecrit_slab_oc.ne.-1) close(97)
    10971068
    10981069!Al1 Call to 1D equivalent of dynredem (an,mois,jour,heure ?)
    1099 ! -------------------------------------
     1070! ---------------------------------------------------------------------------
    11001071       call dyn1dredem("restart1dyn.nc",                                    &
    11011072     &              plev,play,phi,phis,presnivs,                            &
Note: See TracChangeset for help on using the changeset viewer.