Changeset 3905 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

Timestamp:

2013-05-24T15:40:39+02:00 (11 years ago)

Author:

cetlod

Message:

bugfix : store the flux due to freezing/melting, needed at least for PISCES, see ticket #1100

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

@@ -105 +105 @@
             DO ji = 1, jpi
                zdep    = rfact2 / fse3t(ji,jj,1)
-    !           zwflux  = ( emps(ji,jj) - emp(ji,jj) ) &
-    !           &        * tsn(ji,jj,1,jp_sal) / ( tsn(ji,jj,1,jp_sal) - 6.0 ) / 1000.
-               zwflux = 0.
-               zfminus = MIN( 0., -zwflux ) * trn(ji,jj,1,jpfer) * zdep
-               zfplus  = MAX( 0., -zwflux ) * 10E-9 * zdep
+               zwflux  = fmmflx(ji,jj) / 1000._wp
+               zfminus = MIN( 0._wp, -zwflux ) * trn(ji,jj,1,jpfer) * zdep
+               zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep
                zironice(ji,jj) =  zfplus + zfminus
             END DO