Changeset 5469

Timestamp:

10/08/20 14:27:16 (4 years ago)

Author:

oboucher

Message:

change of rho according to Gurvan's

File:

• TOOLS/WATER_BUDGET/waterbudget.sh (modified) (4 diffs)

TOOLS/WATER_BUDGET/waterbudget.sh

@@ -6 +6 @@
 #--30/11/2016
 #--updated 10/05/2019
+#--updated July 2020 => Gurvan says dOCE_vol has to be computed with rho_oce/rho_liq
 
 #--only work for NEMO 1 degree ! mods needed for NEMO025
@@ -94 +95 @@
 rho_oce=1027.
 #--average ratio of salty water to fresh water
-fresh2salt=1.035
+#fresh2salt=1.035
 #--conversion from /day to /s 
 day2sec=86400.
@@ -212 +213 @@
   echo 'OCE vol beg =' ${OCE_vol_daybeg} >> ${fileout}
   echo 'OCE vol end =' ${OCE_vol_dayend} >> ${fileout}
-  dOCE_vol=`python -c "print (${OCE_vol_dayend}-${OCE_vol_daybeg})/${nbyear}/${nbdayinyr}/${day2sec}/1.e6"`
+  dOCE_vol=`python -c "print  ${rho_oce}/${rho_liq}*(${OCE_vol_dayend}-${OCE_vol_daybeg})/${nbyear}/${nbdayinyr}/${day2sec}/1.e6"`
 else
   maxtimestep=`ncdump -h ${file_oce_mthend} | grep currently | grep -o -E '[0-9]+'`
@@ -219 +220 @@
   echo 'OCE vol beg =' ${OCE_vol_mthbeg} >> ${fileout}
   echo 'OCE vol end =' ${OCE_vol_mthend} >> ${fileout}
-  dOCE_vol=`python -c "print (${OCE_vol_mthend}-${OCE_vol_mthbeg})/${nbyear}/${nbdayinyr}/${day2sec}/1.e6"`
+  dOCE_vol=`python -c "print  ${rho_oce}/${rho_liq}*(${OCE_vol_mthend}-${OCE_vol_mthbeg})/${nbyear}/${nbdayinyr}/${day2sec}/1.e6"`
 fi
 echo 'These estimates of d OCE (converted to Sv) are approximate'  >> ${fileout}