Changeset 7484

Timestamp:

2016-12-10T17:27:38+01:00 (7 years ago)

Author:

cetlod

Message:

PISCES gas : minor bug correction

File:

• branches/CNRS/dev_r6526_PISCES_GAS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90 (modified) (2 diffs)

branches/CNRS/dev_r6526_PISCES_GAS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

@@ -221 +221 @@
                zD = ( 10**( -1010. / ( ztc + 273.15 ) - 1.3246) ) * 1.E-4 !  cm2/s to m2/s
                zsch_cos = zkin_vis / zD
+              !  wind speed 
+               zws  = wndm(ji,jj) * wndm(ji,jj)
+               ! Compute the piston velocity for O2 and CO2
+               zkgwan = 0.251 * zws
+               zkgwan = zkgwan * xconv * ( 1.- fr_i(ji,jj) ) * tmask(ji,jj,1)
+# if defined key_degrad
+               zkgwan = zkgwan * facvol(ji,jj,1)
+#endif 
+               ! compute gas exchange for DMS, CO, Isoprene, COS and N2O
                zkgdms(ji,jj) = zkgwan * SQRT( 660./ zsch_dms )
                zkgco (ji,jj) = zkgwan * SQRT( 660./ zsch_co )
@@ -334 +343 @@
          !
          IF (lk_gas) THEN
-         IF ( iom_use( "DMSflx" ) )  THEN
-           zw2d(:,:) = zdmsflx(:,:) * 1000 * tmask(:,:,1)
-           CALL iom_put( "DMSflx"  , zw2d )
-         ENDIF 
-         IF ( iom_use( "COflx" ) )  THEN
-           zw2d(:,:) = zcoflx(:,:) * 1000 * tmask(:,:,1)
-           CALL iom_put( "COflx"  , zw2d )
-         ENDIF
-         IF ( iom_use( "ISPflx" ) )  THEN
-           zw2d(:,:) = zispflx(:,:) * 1000 * tmask(:,:,1)
-           CALL iom_put( "ISPflx"  , zw2d )
-         ENDIF
-         IF ( iom_use( "COSflx" ) )  THEN
-           zw2d(:,:) = zcosflx(:,:) * 1000 * tmask(:,:,1)
-           CALL iom_put( "COSflx"  , zw2d )
-         ENDIF
-         IF ( iom_use( "H_cos" ) )  THEN
-           zw2d(:,:) = zH_cos(:,:) * tmask(:,:,1)
-           CALL iom_put( "H_cos"  , zw2d )
-         ENDIF
-         IF ( iom_use( "N2Oflx" ) )  THEN
-           zw2d(:,:) = znflx(:,:) * 1000 * tmask(:,:,1)
-           CALL iom_put( "N2Oflx"  , zw2d )
-         ENDIF
+            IF ( iom_use( "DMSflx" ) )  THEN
+              zw2d(:,:) = zdmsflx(:,:) * 1000 * tmask(:,:,1)
+              CALL iom_put( "DMSflx"  , zw2d )
+            ENDIF 
+            IF ( iom_use( "COflx" ) )  THEN
+              zw2d(:,:) = zcoflx(:,:) * 1000 * tmask(:,:,1)
+              CALL iom_put( "COflx"  , zw2d )
+            ENDIF
+            IF ( iom_use( "ISPflx" ) )  THEN
+              zw2d(:,:) = zispflx(:,:) * 1000 * tmask(:,:,1)
+              CALL iom_put( "ISPflx"  , zw2d )
+            ENDIF
+            IF ( iom_use( "COSflx" ) )  THEN
+              zw2d(:,:) = zcosflx(:,:) * 1000 * tmask(:,:,1)
+              CALL iom_put( "COSflx"  , zw2d )
+            ENDIF
+            IF ( iom_use( "H_cos" ) )  THEN
+              zw2d(:,:) = zH_cos(:,:) * tmask(:,:,1)
+              CALL iom_put( "H_cos"  , zw2d )
+            ENDIF
+            IF ( iom_use( "N2Oflx" ) )  THEN
+              zw2d(:,:) = znflx(:,:) * 1000 * tmask(:,:,1)
+              CALL iom_put( "N2Oflx"  , zw2d )
+            ENDIF
          ENDIF