New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7484 for branches/CNRS/dev_r6526_PISCES_GAS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90 – NEMO

Ignore:
Timestamp:
2016-12-10T17:27:38+01:00 (8 years ago)
Author:
cetlod
Message:

PISCES gas : minor bug correction

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/CNRS/dev_r6526_PISCES_GAS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

    r7483 r7484  
    221221               zD = ( 10**( -1010. / ( ztc + 273.15 ) - 1.3246) ) * 1.E-4 !  cm2/s to m2/s 
    222222               zsch_cos = zkin_vis / zD 
     223              !  wind speed  
     224               zws  = wndm(ji,jj) * wndm(ji,jj) 
     225               ! Compute the piston velocity for O2 and CO2 
     226               zkgwan = 0.251 * zws 
     227               zkgwan = zkgwan * xconv * ( 1.- fr_i(ji,jj) ) * tmask(ji,jj,1) 
     228# if defined key_degrad 
     229               zkgwan = zkgwan * facvol(ji,jj,1) 
     230#endif  
     231               ! compute gas exchange for DMS, CO, Isoprene, COS and N2O 
    223232               zkgdms(ji,jj) = zkgwan * SQRT( 660./ zsch_dms ) 
    224233               zkgco (ji,jj) = zkgwan * SQRT( 660./ zsch_co ) 
     
    334343         ! 
    335344         IF (lk_gas) THEN 
    336          IF ( iom_use( "DMSflx" ) )  THEN 
    337            zw2d(:,:) = zdmsflx(:,:) * 1000 * tmask(:,:,1) 
    338            CALL iom_put( "DMSflx"  , zw2d ) 
    339          ENDIF  
    340          IF ( iom_use( "COflx" ) )  THEN 
    341            zw2d(:,:) = zcoflx(:,:) * 1000 * tmask(:,:,1) 
    342            CALL iom_put( "COflx"  , zw2d ) 
    343          ENDIF 
    344          IF ( iom_use( "ISPflx" ) )  THEN 
    345            zw2d(:,:) = zispflx(:,:) * 1000 * tmask(:,:,1) 
    346            CALL iom_put( "ISPflx"  , zw2d ) 
    347          ENDIF 
    348          IF ( iom_use( "COSflx" ) )  THEN 
    349            zw2d(:,:) = zcosflx(:,:) * 1000 * tmask(:,:,1) 
    350            CALL iom_put( "COSflx"  , zw2d ) 
    351          ENDIF 
    352          IF ( iom_use( "H_cos" ) )  THEN 
    353            zw2d(:,:) = zH_cos(:,:) * tmask(:,:,1) 
    354            CALL iom_put( "H_cos"  , zw2d ) 
    355          ENDIF 
    356          IF ( iom_use( "N2Oflx" ) )  THEN 
    357            zw2d(:,:) = znflx(:,:) * 1000 * tmask(:,:,1) 
    358            CALL iom_put( "N2Oflx"  , zw2d ) 
    359          ENDIF 
     345            IF ( iom_use( "DMSflx" ) )  THEN 
     346              zw2d(:,:) = zdmsflx(:,:) * 1000 * tmask(:,:,1) 
     347              CALL iom_put( "DMSflx"  , zw2d ) 
     348            ENDIF  
     349            IF ( iom_use( "COflx" ) )  THEN 
     350              zw2d(:,:) = zcoflx(:,:) * 1000 * tmask(:,:,1) 
     351              CALL iom_put( "COflx"  , zw2d ) 
     352            ENDIF 
     353            IF ( iom_use( "ISPflx" ) )  THEN 
     354              zw2d(:,:) = zispflx(:,:) * 1000 * tmask(:,:,1) 
     355              CALL iom_put( "ISPflx"  , zw2d ) 
     356            ENDIF 
     357            IF ( iom_use( "COSflx" ) )  THEN 
     358              zw2d(:,:) = zcosflx(:,:) * 1000 * tmask(:,:,1) 
     359              CALL iom_put( "COSflx"  , zw2d ) 
     360            ENDIF 
     361            IF ( iom_use( "H_cos" ) )  THEN 
     362              zw2d(:,:) = zH_cos(:,:) * tmask(:,:,1) 
     363              CALL iom_put( "H_cos"  , zw2d ) 
     364            ENDIF 
     365            IF ( iom_use( "N2Oflx" ) )  THEN 
     366              zw2d(:,:) = znflx(:,:) * 1000 * tmask(:,:,1) 
     367              CALL iom_put( "N2Oflx"  , zw2d ) 
     368            ENDIF 
    360369         ENDIF 
    361370 
Note: See TracChangeset for help on using the changeset viewer.