Changeset 5449 for BOL/Replay/replay_equip.sh

Timestamp:

Dec 23, 2024, 3:55:47 AM (10 hours ago)

Author:

fhourdin

Message:

For alternate number of vertical layers

File:

• BOL/Replay/replay_equip.sh (modified) (2 diffs)

BOL/Replay/replay_equip.sh

@@ -117 +117 @@
 case $param in
    vdif_k|vdif_cd|vdif|my_25|vdif_dq) param_ini=vdif_ini ; inimod=simple_vdif_ini ; klon=ngrid ; klev=nlay ;;
-   thermcell_main|thermcell_plume_6A|thermcell_env|thermcell_height|thermcell_dry|\
+   thermcell_plume_6A|thermcell_env|thermcell_height|thermcell_dry|\
       thermcell_closure|thermcell_height|thermcell_dq|thermcell_flux2|thermcell_down| \
       thermcell_updown_dq|thermcell_dtke) \
       param_ini=thermcell_ini ; inimod=lmdz_thermcell_ini ; klon=ngrid ; klev=nlay ;;
    wake|pkupper|wake_popdyn_1|wake_popdyn_2|vdif_kcay|ustarhb) param_ini=wake_ini ; inimod=lmdz_wake_ini ; klon=klon ; klev=klev ;;
+   surf_wind) param_ini=surf_wind_ini ; inimod=lmdz_surf_wind_ini ; klon=klon ; klev=nsrfwnd ;;
    ratqs_main|ratqs_inter|ratqs_oro|ratqs_hetero|ratqs_tke) param_ini=ratqs_ini ; inimod=lmdz_ratqs_ini ; klon=klon ; klev=klev ;;
    lscp|fisrtilp) param_ini=lscp_ini ; inimod=lmdz_lscp_ini ; klon=klon ; klev=klev ;;
@@ -351 +352 @@
    line=`sed -n -e "/CALL wake_ini/=" physiq_mod.F90 | head -1`
    if [[ "$lonlat" = "" ]] ; then
-        pattern="CALL iophys_ini(pdtphys)"
+        if [[ $klev == klev || $klev == nlay ]] ; then
+            klev_=klev
+        else
+            klev_=$klev
+        fi
+        pattern="CALL iophys_ini(pdtphys,$klev_)"
    else
         pattern='call iotd_ini("phys.nc",1,1,klev,'$lonlat',presnivs,1,1,1,0.,pdtphys,calend)'
    fi
+
    sed -i"" -e "${line}s/^/   $pattern ! $replay_comment  ! $replay_comment\n/" physiq_mod.F90
 fi