diff --git a/gcm_forecast.tmpl b/gcm_forecast.tmpl index 4b043b28..9c9f3d68 100755 --- a/gcm_forecast.tmpl +++ b/gcm_forecast.tmpl @@ -555,6 +555,18 @@ if( $AGCM_LM != 72 ) then cp $file dummy cat dummy | sed -e "s|/L72/|/L${AGCM_LM}/|g" | sed -e "s|z72|z${AGCM_LM}|g" > $file end + + # We have to do something special for L186. By default all of the GOCART emissions have: + # pressure_lid_in_hPa: 0.01 + # which is the pressure lid in hPa for all other GEOS vertical levels. But L186 has a lid + # at 0.00025 Pa. So we need to change this value in the GOCART emissions files for L186. + if ( $AGCM_LM == 186 ) then + set files = `/bin/ls *.rc` + foreach file ($files) + cp $file dummy + cat dummy | sed -e '/pressure_lid_in_hPa:/c\pressure_lid_in_hPa: 2.5e-6' > $file + end + endif endif # Rename big ExtData files that are not needed diff --git a/gcm_regress.j b/gcm_regress.j index 1ed03381..b02894dc 100755 --- a/gcm_regress.j +++ b/gcm_regress.j @@ -382,6 +382,26 @@ if( $EXTDATA2G_TRUE == 1 ) then endif +if( $AGCM_LM != 72 ) then + set files = `/bin/ls *.yaml` + foreach file ($files) + cp $file dummy + cat dummy | sed -e "s|/L72/|/L${AGCM_LM}/|g" | sed -e "s|z72|z${AGCM_LM}|g" > $file + end + + # We have to do something special for L186. By default all of the GOCART emissions have: + # pressure_lid_in_hPa: 0.01 + # which is the pressure lid in hPa for all other GEOS vertical levels. But L186 has a lid + # at 0.00025 Pa. So we need to change this value in the GOCART emissions files for L186. + if ( $AGCM_LM == 186 ) then + set files = `/bin/ls *.rc` + foreach file ($files) + cp $file dummy + cat dummy | sed -e '/pressure_lid_in_hPa:/c\pressure_lid_in_hPa: 2.5e-6' > $file + end + endif +endif + # Move GOCART to use RRTMGP Bands # ------------------------------- # UNCOMMENT THE LINES BELOW IF RUNNING RRTMGP diff --git a/gcm_run.j b/gcm_run.j index 4daa7804..efa592a1 100755 --- a/gcm_run.j +++ b/gcm_run.j @@ -776,6 +776,18 @@ if( $AGCM_LM != 72 ) then cp $file dummy cat dummy | sed -e "s|/L72/|/L${AGCM_LM}/|g" | sed -e "s|z72|z${AGCM_LM}|g" > $file end + + # We have to do something special for L186. By default all of the GOCART emissions have: + # pressure_lid_in_hPa: 0.01 + # which is the pressure lid in hPa for all other GEOS vertical levels. But L186 has a lid + # at 0.00025 Pa. So we need to change this value in the GOCART emissions files for L186. + if ( $AGCM_LM == 186 ) then + set files = `/bin/ls *.rc` + foreach file ($files) + cp $file dummy + cat dummy | sed -e '/pressure_lid_in_hPa:/c\pressure_lid_in_hPa: 2.5e-6' > $file + end + endif endif # Rename big ExtData files that are not needed diff --git a/gcm_run_benchmark.j b/gcm_run_benchmark.j index d3f1c871..22c97ee0 100755 --- a/gcm_run_benchmark.j +++ b/gcm_run_benchmark.j @@ -649,6 +649,26 @@ if( ${EMISSIONS} == AMIP_EMISSIONS ) then endif +if( $AGCM_LM != 72 ) then + set files = `/bin/ls *.yaml` + foreach file ($files) + cp $file dummy + cat dummy | sed -e "s|/L72/|/L${AGCM_LM}/|g" | sed -e "s|z72|z${AGCM_LM}|g" > $file + end + + # We have to do something special for L186. By default all of the GOCART emissions have: + # pressure_lid_in_hPa: 0.01 + # which is the pressure lid in hPa for all other GEOS vertical levels. But L186 has a lid + # at 0.00025 Pa. So we need to change this value in the GOCART emissions files for L186. + if ( $AGCM_LM == 186 ) then + set files = `/bin/ls *.rc` + foreach file ($files) + cp $file dummy + cat dummy | sed -e '/pressure_lid_in_hPa:/c\pressure_lid_in_hPa: 2.5e-6' > $file + end + endif +endif + # Rename big ExtData files that are not needed # -------------------------------------------- set SC_TRUE = `grep -i '^\s*ENABLE_STRATCHEM:\s*\.TRUE\.' GEOS_ChemGridComp.rc | wc -l`