Changeset 1858 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/TRA/trabbc.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/trabbc.F90 (modified) (2 diffs)

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

@@ -71 +71 @@
       !!----------------------------------------------------------------------
       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
@@ -98 +98 @@
             DO ji = 2, jpim1
 #endif
-               zqgh_trd = ro0cpr * qgh_trd0(ji,jj) / fse3t(ji,jj,nbotlevt(ji,jj))
+               zqgh_trd = r1_rau0_rcp * qgh_trd0(ji,jj) / fse3t(ji,jj,nbotlevt(ji,jj))
                ta(ji,jj,nbotlevt(ji,jj)) = ta(ji,jj,nbotlevt(ji,jj)) + zqgh_trd
             END DO