Changeset 1858 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/TRA/trasbc.F90

Timestamp:

2010-05-04T10:39:48+02:00 (14 years ago)

Author:

gm

Message:

ticket:#665 : step 1 - heat content of freezing-melting ice

File:

• branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/TRA/trasbc.F90 (modified) (2 diffs)

branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -99 +99 @@
       !!----------------------------------------------------------------------
       USE oce, ONLY :   ztrdt => ua   ! use ua as 3D workspace   
-      USE oce, ONLY :   ztrds => va   ! use va as 3D workspace   
+      USE oce, ONLY :   ztrds => va   !  -  va     -      -   
       !!
       INTEGER, INTENT(in) ::   kt     ! ocean time-step index
@@ -132 +132 @@
 #endif
             IF( lk_vvl) THEN
-               zta = ro0cpr * qns(ji,jj) * zse3t &                   ! temperature : heat flux
+               zta = r1_rau0_rcp * qns(ji,jj) * zse3t &              ! temperature : heat flux
                 &    - emp(ji,jj) * zsrau * tn(ji,jj,1)  * zse3t     ! & cooling/heating effet of EMP flux
                zsa = 0.e0                                            ! No salinity concent./dilut. effect
             ELSE
-               zta = ro0cpr * qns(ji,jj) * zse3t     ! temperature : heat flux
+               zta = r1_rau0_rcp * qns(ji,jj) * zse3t                ! temperature : heat flux
                zsa = emps(ji,jj) * zsrau * sn(ji,jj,1)   * zse3t     ! salinity :  concent./dilut. effect
             ENDIF