Ignore:
Timestamp:
Jul 23, 2024, 3:29:36 PM (8 weeks ago)
Author:
abarral
Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/1D_nudge_sandu_astex.h

    r5099 r5103  
    1       do l = 1, llm
    2 
    3 ! au dessus de 700hPa, on relaxe vers profil init
    4 !      on fait l'hypothese que dans ce cas, il n'y a plus d'eau liq. au dessus 700hpa
    5 !      donc la relaxation en thetal et qt devient relaxation en tempe et qv
    6        relax_u(l)=(u(l)-u_mod(l))/tau_sandu  ! pour u et v on relaxe sur tte la colonne
    7        relax_v(l)=(v(l)-v_mod(l))/tau_sandu
    8        relax_q(l,1)=0.
    9        relax_q(l,2)=0.
    10        relax_thl(l)=0.
    11 !      print *,'nudge: l tau_sandu u u_mod',l,tau_sandu,u(l),u_mod(l)
    12 
    13        if (l.ge.llm700) then
    14          relax_q(l,1)=(q(l,1)-q_mod(l))/tau_sandu
    15          relax_q(l,2)=0.
    16          relax_thl(l)=(temp(l)-t_mod(l))/tau_sandu
    17        endif
    18 !     print *,'l dq1 relax dqadv',l,dq(l,1),relax_q(l,1),d_q_adv(l,1)
    19 !     print *,'l dq2 relax dqadv',l,dq(l,2),relax_q(l,2),d_q_adv(l,2)
    20 !     print *,'l dt  relax dtadv',l,dt_phys(l),relax_thl(l),d_t_adv(l)
    21       enddo
    22 
    23         u(1:mxcalc)=u(1:mxcalc) + timestep*( du_age(1:mxcalc)+             &         
    24      &              du_phys(1:mxcalc) - relax_u(1:mxcalc))
    25         v(1:mxcalc)=v(1:mxcalc) + timestep*( dv_age(1:mxcalc)+              &
    26      &               dv_phys(1:mxcalc) - relax_v(1:mxcalc))
    27         q(1:mxcalc,1)=q(1:mxcalc,1)+timestep*(                              &
    28      &      dq(1:mxcalc,1) - relax_q(1:mxcalc,1)+d_q_adv(1:mxcalc,1))
    29         q(1:mxcalc,2)=q(1:mxcalc,2)+timestep*(                              &
    30      &      dq(1:mxcalc,2) - relax_q(1:mxcalc,2)+d_q_adv(1:mxcalc,2))
    31         temp(1:mxcalc)=temp(1:mxcalc)+timestep*(                            &
    32      &      dt_phys(1:mxcalc)-relax_thl(1:mxcalc)+d_t_adv(1:mxcalc))
    33 
    34 
    35       ! print*,'OLDLMDZ1D IOPH'
    36       !      CALL iophys_ecrit('relax_thl',klev,'relax_thl','m/s',relax_thl)
    37       !      CALL iophys_ecrit('d_t_adv',klev,'d_t_adv','m/s',d_t_adv)
    38       !      CALL iophys_ecrit('temp',klev,'temp','m/s',temp)
    39       !      CALL iophys_ecrit('q',klev,'q','m/s',q(:,1))
    40       !      CALL iophys_ecrit('relax_q',klev,'relax_q','m/s',relax_q(:,1))
    41       !      CALL iophys_ecrit('d_q_adv',klev,'d_q_adv','m/s',d_q_adv(:,1))
    42 
Note: See TracChangeset for help on using the changeset viewer.