Hi guys, I’m currently calculating the density differences across the open ocean and have run into a little issue... Slightly confused myself by finding two methods for calculating density on the expressions page, am I correct in thinking that: yes? Use “./ECP85_sal.nc” yes? Use “./ECP85_temp.nc” yes? SET REGION/I=1:36/J=1:36 yes? LET salt = ocn_sal[d=1] yes? LET y250strat = dens[k=1,l=250] – dens[k=10,l=250] yes? List y250strat Would be less accurate than calculating with the following script, because it’s generating the density based on surface water salinity and temperature values?: yes? use "./ECP85_sal.nc" yes? use "./ECP85_temp.nc" yes? let temp=ocn_temp[d=2]*1.00024 !convert to different temp scale for accuracy yes? let salt=ocn_sal[d=1] yes? let potemp=THETA_FO(SALT, TEMP,Z[G=TEMP],0) !Define the potential temperature for surface ocean yes? let potemp2=THETA_FO(SALT,TEMP,Z[G=TEMP],2303.5) !Define the potential temperature for benthic ocean yes? let surfdens=rho_un(salt,potemp,0) !calculate density at surface 0 = dBar pressure yes? let benthdens=rho_un(salt,potemp2,2303.5) !calculate density at 2303.5m deep, with 2303.5dBar pressure yes? let y250strat= surfdens[k=1,l=250] - benthdens[k=12,l=250] yes? list y250strat Also is there anything obviously wrong with the script I’ve made? First time calculating ocean densities! Cheers, Josh Sent from Mail for Windows 10 |