Ignore:
Timestamp:
08/31/12 15:41:37 (12 years ago)
Author:
cholod
Message:

Load NEMO_TMP into vendor/nemo/current.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • vendor/nemo/current/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zlys.F90

    r4 r44  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    44    !! $Id: p4zlys.F90 3295 2012-01-30 15:49:07Z cetlod $  
     44   !! $Id: p4zlys.F90 3321 2012-03-05 17:10:55Z cetlod $  
    4545   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
    … …  
    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.