Ignore:
Timestamp:
Oct 24, 2016, 6:46:30 PM (8 years ago)
Author:
fhourdin
Message:

Corrections pour le cas Dice (Marie-Pierre Lefebvre)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

    r2672 r2683  
    367367      fich_dice='dice_driver.nc'
    368368      call read_dice(fich_dice,nlev_dice,nt_dice                    &
    369      &     ,zz_dice,plev_dice,th_dice,qv_dice,u_dice,v_dice,o3_dice &
     369     &     ,zz_dice,plev_dice,t_dice,qv_dice,u_dice,v_dice,o3_dice &
    370370     &     ,shf_dice,lhf_dice,lwup_dice,swup_dice,tg_dice,ustar_dice&
    371371     &     ,psurf_dice,ug_dice,vg_dice,ht_dice,hq_dice              &
     
    376376!champs initiaux:
    377377      do k=1,nlev_dice
    378          th_dicei(k)=th_dice(k)
     378         t_dicei(k)=t_dice(k)
    379379         qv_dicei(k)=qv_dice(k)
    380380         u_dicei(k)=u_dice(k)
     
    405405
    406406      CALL interp_dice_vertical(play,nlev_dice,nt_dice,plev_dice       &
    407      &         ,th_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei             &
     407     &         ,t_dicei,qv_dicei,u_dicei,v_dicei,o3_dicei             &
    408408     &         ,ht_dicei,hq_dicei,hu_dicei,hv_dicei,w_dicei,omega_dicei&
    409      &         ,th_mod,qv_mod,u_mod,v_mod,o3_mod                       &
     409     &         ,t_mod,qv_mod,u_mod,v_mod,o3_mod                       &
    410410     &         ,ht_mod,hq_mod,hu_mod,hv_mod,w_mod,omega_mod,mxcalc)
    411411
     
    425425      do l = 1, llm
    426426! Ligne du dessous ?? decommenter si on lit theta au lieu de temp
    427        temp(l) = th_mod(l)*(play(l)/pzero)**rkappa
    428 !      temp(l) = t_mod(l)
     427!      temp(l) = th_mod(l)*(play(l)/pzero)**rkappa
     428       temp(l) = t_mod(l)
    429429       q(l,1) = qv_mod(l)
    430430       q(l,2) = 0.0
     
    476476!---------------------------------------------------------------------
    477477
    478 !!!! Si la temperature de surface n'est pas imposée:
     478!!!! Si la temperature de surface n'est pas impos??e:
    479479 
    480480      if (forcing_gabls4) then
     
    531531      write(*,*) 'SST initiale mxcalc: ',tsurf,mxcalc
    532532      do l = 1, llm
    533 ! Ligne du dessous à decommenter si on lit theta au lieu de temp
     533! Ligne du dessous ?? decommenter si on lit theta au lieu de temp
    534534!      temp(l) = th_mod(l)*(play(l)/pzero)**rkappa
    535535       temp(l) = t_mod(l)
     
    909909      endif !forcing_case
    910910!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     911
Note: See TracChangeset for help on using the changeset viewer.