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 2963 for branches/2011/dev_r2787_PISCES_improvment/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90 – NEMO

Ignore:
Timestamp:
2011-10-19T14:16:52+02:00 (13 years ago)
Author:
cetlod
Message:

dev_r2787_PISCES_improvment:include temporal evolution of ocean volume ( vvl )

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/dev_r2787_PISCES_improvment/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

    r2960 r2963  
    144144      REAL(wp) ::  silmean = 91.51     ! mean value of silicate 
    145145      ! 
    146       INTEGER  :: jk 
    147146      REAL(wp) :: zarea, zalksum, zpo4sum, zno3sum, zsilsum 
    148       REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zvolume     ! volume of ocean 
    149147      !!--------------------------------------------------------------------- 
    150148 
    … …  
    157155         !                                                    ! --------------------------- ! 
    158156         ! set total alkalinity, phosphate, nitrate & silicate 
    159          ALLOCATE( zvolume(jpi,jpj,jpk) ) 
    160          !                                                              ! masked grid volume 
    161          DO jk = 1, jpk 
    162             zvolume(:,:,jk) = e1e2t(:,:)  * fse3t(:,:,jk) * tmask(:,:,jk)    
    163          ENDDO  
    164 #if defined key_degrad 
    165          zvolume(:,:,:) = zvolume(:,:,:) * facvol(:,:,:) 
    166 #endif 
    167          zarea          = 1._wp / glob_sum( zvolume(:,:,:) ) * 1e6               
    168  
    169          zalksum = glob_sum( trn(:,:,:,jptal) * zvolume(:,:,:)  ) * zarea 
    170          zpo4sum = glob_sum( trn(:,:,:,jppo4) * zvolume(:,:,:)  ) * zarea / 122. 
    171          zno3sum = glob_sum( trn(:,:,:,jpno3) * zvolume(:,:,:)  ) * zarea / 7.6 
    172          zsilsum = glob_sum( trn(:,:,:,jpsil) * zvolume(:,:,:)  ) * zarea 
     157         zarea          = 1._wp / glob_sum( cvol(:,:,:) ) * 1e6               
     158 
     159         zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea 
     160         zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea / 122. 
     161         zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea / 7.6 
     162         zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea 
    173163  
    174164         IF(lwp) WRITE(numout,*) '       TALK mean : ', zalksum 
Note: See TracChangeset for help on using the changeset viewer.