Changeset 3321


Ignore:
Timestamp:
2012-03-05T18:10:55+01:00 (9 years ago)
Author:
cetlod
Message:

trunk:bugfix on pH diagnostic in PISCES, see ticket #931

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zlys.F90

    r3295 r3321  
    150150         IF( lk_iomput ) THEN 
    151151            zrfact2 = 1.e3 * rfact2r 
    152             CALL iom_put( "PH"    , hi      (:,:,:)                  * tmask(:,:,:) ) 
    153             CALL iom_put( "CO3"   , zco3    (:,:,:) * 1e+3           * tmask(:,:,:) ) 
    154             CALL iom_put( "CO3sat", aksp    (:,:,:) * 1e+3 / calcon  * tmask(:,:,:) ) 
    155             CALL iom_put( "DCAL"  , zcaldiss(:,:,:) * zrfact2 * tmask(:,:,:) ) 
     152            CALL iom_put( "PH"    , -1. * LOG10( hi(:,:,:) )                * tmask(:,:,:) ) 
     153            CALL iom_put( "CO3"   ,        zco3    (:,:,:) * 1e+3           * tmask(:,:,:) ) 
     154            CALL iom_put( "CO3sat",        aksp    (:,:,:) * 1e+3 / calcon  * tmask(:,:,:) ) 
     155            CALL iom_put( "DCAL"  ,        zcaldiss(:,:,:) * zrfact2        * tmask(:,:,:) ) 
    156156         ELSE 
    157             trc3d(:,:,:,jp_pcs0_3d    ) = hi  (:,:,:)          * tmask(:,:,:) 
    158             trc3d(:,:,:,jp_pcs0_3d + 1) = zco3(:,:,:)          * tmask(:,:,:) 
    159             trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon * tmask(:,:,:) 
     157            trc3d(:,:,:,jp_pcs0_3d    ) = -1. * LOG10( hi(:,:,:) ) * tmask(:,:,:) 
     158            trc3d(:,:,:,jp_pcs0_3d + 1) = zco3(:,:,:)              * tmask(:,:,:) 
     159            trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon     * tmask(:,:,:) 
    160160         ENDIF 
    161161         ! 
Note: See TracChangeset for help on using the changeset viewer.