[Dart-dev] [5918] DART/branches/development/models/POP/model_mod.f90: The conversion routine wants salinity in psu (g/kg) and

[nancy at ucar.edu]

Revision: 5918
Author:   nancy
Date:     2012-11-05 09:52:56 -0700 (Mon, 05 Nov 2012)
Log Message:
-----------
The conversion routine wants salinity in psu (g/kg) and
the model has it in msu (kg/kg).  multiply by 1000 before
doing the potential temp -> in-situ temp conversion.

Modified Paths:
--------------
    DART/branches/development/models/POP/model_mod.f90

-------------- next part --------------
Modified: DART/branches/development/models/POP/model_mod.f90
===================================================================
--- DART/branches/development/models/POP/model_mod.f90	2012-11-05 16:51:48 UTC (rev 5917)
+++ DART/branches/development/models/POP/model_mod.f90	2012-11-05 16:52:56 UTC (rev 5918)
@@ -3370,13 +3370,13 @@
 endif
 
 
-! salinity
+! salinity - in msu (kg/kg).  converter will want psu (g/kg).
 call do_interp(x, start_index(S_index), hgt_bot, hgt_top, hgt_fract, llon, llat, &
                KIND_SALINITY, salinity_val, istatus)
 if(istatus /= 0) return
 if (debug > 8) print *, 'salinity: ', salinity_val
 
-! potential temperature
+! potential temperature - degrees C.
 call do_interp(x, start_index(T_index), hgt_bot, hgt_top, hgt_fract, llon, llat, &
                KIND_POTENTIAL_TEMPERATURE, potential_temp, istatus)
 if(istatus /= 0) return
@@ -3391,8 +3391,9 @@
                                  hgt_top, pressure(hgt_top), pres_val
 endif
 
-! and finally, convert to sensible (in-situ) temperature
-call insitu_temp(potential_temp, salinity_val, pres_val*10.0_r8, interp_val)
+! and finally, convert to sensible (in-situ) temperature.
+! potential temp in degrees C, pressure in decibars, salinity in psu or pss (g/kg).
+call insitu_temp(potential_temp, salinity_val*1000.0_r8, pres_val*10.0_r8, interp_val)
 if (debug > 2) print *, 's,pt,pres,t: ', salinity_val, potential_temp, pres_val, interp_val
 
 end subroutine compute_temperature