- Timestamp:
- 2016-12-10T17:27:38+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6526_PISCES_GAS/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90
r7483 r7484 221 221 zD = ( 10**( -1010. / ( ztc + 273.15 ) - 1.3246) ) * 1.E-4 ! cm2/s to m2/s 222 222 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 223 232 zkgdms(ji,jj) = zkgwan * SQRT( 660./ zsch_dms ) 224 233 zkgco (ji,jj) = zkgwan * SQRT( 660./ zsch_co ) … … 334 343 ! 335 344 IF (lk_gas) THEN 336 IF ( iom_use( "DMSflx" ) ) THEN337 zw2d(:,:) = zdmsflx(:,:) * 1000 * tmask(:,:,1)338 CALL iom_put( "DMSflx" , zw2d )339 ENDIF340 IF ( iom_use( "COflx" ) ) THEN341 zw2d(:,:) = zcoflx(:,:) * 1000 * tmask(:,:,1)342 CALL iom_put( "COflx" , zw2d )343 ENDIF344 IF ( iom_use( "ISPflx" ) ) THEN345 zw2d(:,:) = zispflx(:,:) * 1000 * tmask(:,:,1)346 CALL iom_put( "ISPflx" , zw2d )347 ENDIF348 IF ( iom_use( "COSflx" ) ) THEN349 zw2d(:,:) = zcosflx(:,:) * 1000 * tmask(:,:,1)350 CALL iom_put( "COSflx" , zw2d )351 ENDIF352 IF ( iom_use( "H_cos" ) ) THEN353 zw2d(:,:) = zH_cos(:,:) * tmask(:,:,1)354 CALL iom_put( "H_cos" , zw2d )355 ENDIF356 IF ( iom_use( "N2Oflx" ) ) THEN357 zw2d(:,:) = znflx(:,:) * 1000 * tmask(:,:,1)358 CALL iom_put( "N2Oflx" , zw2d )359 ENDIF345 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 360 369 ENDIF 361 370
Note: See TracChangeset
for help on using the changeset viewer.