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 3751 for branches/2012/dev_MERGE_2012 – NEMO

Ignore:
Timestamp:
2013-01-14T17:25:50+01:00 (11 years ago)
Author:
cetlod
Message:

minor bug correction in PISCES

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

    r3680 r3751  
    148148         IF( ln_diatrc ) THEN 
    149149            zfact = 1.e+3 * rfact2r 
    150             IF( lk_iomput .AND. jnt == nrdttrc ) THEN 
    151                CALL iom_put( "Irondep", zirondep(:,:,1) * zfact * fse3t(:,:,1) * tmask(:,:,1) ) ! surface downward dust depo of iron 
    152                CALL iom_put( "pdust"  , dust(:,:) / ( wdust * 30.42 * 0.035 )  * tmask(:,:,1) ) ! dust concentration at surface 
     150            IF( lk_iomput ) THEN 
     151               IF( jnt == nrdttrc ) THEN 
     152                  CALL iom_put( "Irondep", zirondep(:,:,1) * zfact * fse3t(:,:,1) * tmask(:,:,1) ) ! surface downward dust depo of iron 
     153                  CALL iom_put( "pdust"  , dust(:,:) / ( wdust * 30.42 * 0.035 )  * tmask(:,:,1) ) ! dust concentration at surface 
     154               ENDIF 
    153155            ELSE 
    154156               trc2d(:,:,jp_pcs0_2d + 11) = zirondep(:,:,1) * zfact * fse3t(:,:,1) * tmask(:,:,1) 
     
    363365      IF( ln_diatrc ) THEN 
    364366         zfact = 1.e+3 * rfact2r 
    365          IF( lk_iomput .AND. jnt == nrdttrc ) THEN 
    366             CALL iom_put( "Nfix"  , znitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) )  ! nitrogen fixation  
    367             CALL iom_put( "Sdenit", zwork4(:,:)               * rno3 * zfact * tmask(:,:,1) )  ! Nitrate reduction in the sediments 
     367         IF( lk_iomput ) THEN 
     368            IF( jnt == nrdttrc ) THEN 
     369               CALL iom_put( "Nfix"  , znitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) )  ! nitrogen fixation  
     370               CALL iom_put( "Sdenit", zwork4(:,:)               * rno3 * zfact * tmask(:,:,1) )  ! Nitrate reduction in the sediments 
     371            ENDIF 
    368372         ELSE 
    369373            trc2d(:,:,jp_pcs0_2d + 12) = znitrpot(:,:,1) * nitrfix * zfact * fse3t(:,:,1) * tmask(:,:,1) 
Note: See TracChangeset for help on using the changeset viewer.