Changeset 5449


Ignore:
Timestamp:
Dec 23, 2024, 3:55:47 AM (5 hours ago)
Author:
fhourdin
Message:

For alternate number of vertical layers

Location:
BOL/Replay
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • BOL/Replay/replay_clean.sh

    r5444 r5449  
    77   \rm -f dump_replay*.h dump_param.h dump_ini.h  input tmp_get_subroutine_arg  tmp_include_line  tmp_include_line*  tmp_include_line*_cpp  tmp_include_line*
    88
    9    for file in `grep "$replay_comment" *${f90} | cut -d: -f1` ; do
     9   for file in `grep "$replay_comment" *f90 *F90 | cut -d: -f1` ; do
    1010      sed -i"" -e "/$replay_comment/d" $file
    1111   done
  • BOL/Replay/replay_equip.sh

    r5442 r5449  
    117117case $param in
    118118   vdif_k|vdif_cd|vdif|my_25|vdif_dq) param_ini=vdif_ini ; inimod=simple_vdif_ini ; klon=ngrid ; klev=nlay ;;
    119    thermcell_main|thermcell_plume_6A|thermcell_env|thermcell_height|thermcell_dry|\
     119   thermcell_plume_6A|thermcell_env|thermcell_height|thermcell_dry|\
    120120      thermcell_closure|thermcell_height|thermcell_dq|thermcell_flux2|thermcell_down| \
    121121      thermcell_updown_dq|thermcell_dtke) \
    122122      param_ini=thermcell_ini ; inimod=lmdz_thermcell_ini ; klon=ngrid ; klev=nlay ;;
    123123   wake|pkupper|wake_popdyn_1|wake_popdyn_2|vdif_kcay|ustarhb) param_ini=wake_ini ; inimod=lmdz_wake_ini ; klon=klon ; klev=klev ;;
     124   surf_wind) param_ini=surf_wind_ini ; inimod=lmdz_surf_wind_ini ; klon=klon ; klev=nsrfwnd ;;
    124125   ratqs_main|ratqs_inter|ratqs_oro|ratqs_hetero|ratqs_tke) param_ini=ratqs_ini ; inimod=lmdz_ratqs_ini ; klon=klon ; klev=klev ;;
    125126   lscp|fisrtilp) param_ini=lscp_ini ; inimod=lmdz_lscp_ini ; klon=klon ; klev=klev ;;
     
    351352   line=`sed -n -e "/CALL wake_ini/=" physiq_mod.F90 | head -1`
    352353   if [[ "$lonlat" = "" ]] ; then
    353         pattern="CALL iophys_ini(pdtphys)"
     354        if [[ $klev == klev || $klev == nlay ]] ; then
     355            klev_=klev
     356        else
     357            klev_=$klev
     358        fi
     359        pattern="CALL iophys_ini(pdtphys,$klev_)"
    354360   else
    355361        pattern='call iotd_ini("phys.nc",1,1,klev,'$lonlat',presnivs,1,1,1,0.,pdtphys,calend)'
    356362   fi
     363
    357364   sed -i"" -e "${line}s/^/   $pattern ! $replay_comment  ! $replay_comment\n/" physiq_mod.F90
    358365fi
Note: See TracChangeset for help on using the changeset viewer.