Changeset 10127 for NEMO/trunk/src


Ignore:
Timestamp:
2018-09-13T17:27:54+02:00 (2 years ago)
Author:
cetlod
Message:

Bugfix to allow the use of on-the-fly interpolation for PISCES ancillary data, see ticket #2003

Location:
NEMO/trunk/src/TOP/PISCES/P4Z
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zsbc.F90

    r10111 r10127  
    7575   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   hydrofe          !: Hydrothermal vent supply of iron 
    7676 
    77    REAL(wp), PUBLIC ::   rivalkinput, rivdicinput, nitdepinput, sumdepsi 
    78    REAL(wp), PUBLIC ::   rivdininput, rivdipinput, rivdsiinput 
     77   REAL(wp), PUBLIC :: sedsilfrac, sedcalfrac 
     78   REAL(wp), PUBLIC :: rivalkinput, rivdicinput 
     79   REAL(wp), PUBLIC :: rivdininput, rivdipinput, rivdsiinput 
    7980 
    8081   !! * Substitutions 
     
    110111         IF( kt == nit000 .OR. ( kt /= nit000 .AND. ntimes_dust > 1 ) ) THEN 
    111112            CALL fld_read( kt, 1, sf_dust ) 
    112             IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN   ;   dust(:,:) = sf_dust(1)%fnow(:,:,1) 
    113             ELSE                                               ;   dust(:,:) = sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) 
     113            IF( nn_ice_tr == -1 .AND. .NOT.ln_ironice ) THEN   ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) ) 
     114            ELSE                                               ;   dust(:,:) = MAX( rtrn, sf_dust(1)%fnow(:,:,1) * ( 1.-fr_i(:,:) ) ) 
    114115            ENDIF 
    115116         ENDIF 
     
    174175             zcoef = rno3 * 14E6 * ryyss 
    175176             CALL fld_read( kt, 1, sf_ndepo ) 
    176              nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1)  
     177             nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) ) 
    177178         ENDIF 
    178179         IF( .NOT.ln_linssh ) THEN 
    179180           zcoef = rno3 * 14E6 * ryyss 
    180            nitdep(:,:) = sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1)  
     181           nitdep(:,:) = MAX( rtrn, sf_ndepo(1)%fnow(:,:,1) / zcoef / e3t_n(:,:,1) ) 
    181182         ENDIF 
    182183      ENDIF 
     
    210211      REAL(wp), DIMENSION(nbtimes) :: zsteps                 ! times records 
    211212      REAL(wp), DIMENSION(:), ALLOCATABLE :: rivinput 
    212       REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zdust, zndepo, zriver, zcmask 
     213      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zriver, zcmask 
    213214      ! 
    214215      CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files 
     
    309310            CALL iom_open (  TRIM( sn_dust%clname ) , numdust ) 
    310311            CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust)  ! get number of record in file 
    311             ALLOCATE( zdust(jpi,jpj,ntimes_dust) ) 
    312             DO jm = 1, ntimes_dust 
    313                CALL iom_get( numdust, jpdom_data, TRIM( sn_dust%clvar ), zdust(:,:,jm), jm ) 
    314             END DO 
    315             CALL iom_close( numdust ) 
    316             ztimes_dust = 1._wp / REAL(ntimes_dust, wp)  
    317             sumdepsi = 0.e0 
    318             DO jm = 1, ntimes_dust 
    319                sumdepsi = sumdepsi + glob_sum( zdust(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_dust ) 
    320             END DO 
    321             sumdepsi = sumdepsi / ( nyear_len(1) * rday ) * 12. * 8.8 * 0.075 * mfrac / 28.1  
    322             DEALLOCATE( zdust) 
    323          ENDIF 
    324       ELSE 
    325          sumdepsi  = 0._wp 
     312         END IF 
    326313      END IF 
    327314 
     
    421408            CALL iom_open ( TRIM( sn_ndepo%clname ), numdepo ) 
    422409            CALL iom_gettime( numdepo, zsteps, kntime=ntimes_ndep) 
    423             ALLOCATE( zndepo(jpi,jpj,ntimes_ndep) ) 
    424             DO jm = 1, ntimes_ndep 
    425                CALL iom_get( numdepo, jpdom_data, TRIM( sn_ndepo%clvar ), zndepo(:,:,jm), jm ) 
    426             END DO 
    427             CALL iom_close( numdepo ) 
    428             ztimes_ndep = 1._wp / REAL(ntimes_ndep, wp)  
    429             nitdepinput = 0._wp 
    430             DO jm = 1, ntimes_ndep 
    431               nitdepinput = nitdepinput + glob_sum( zndepo(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_ndep ) 
    432             ENDDO 
    433             nitdepinput = nitdepinput / rno3 / 14E6  
    434             DEALLOCATE( zndepo) 
    435          ENDIF 
    436       ELSE 
    437          nitdepinput = 0._wp 
     410         ENDIF 
    438411      ENDIF 
    439412 
     
    527500         WRITE(numout,*) '    DIC Supply : ', rivdicinput*1E3*12./1E12     ,' TgC/yr' 
    528501         WRITE(numout,*)  
    529          WRITE(numout,*) '    Total input of elements from atmospheric supply' 
    530          WRITE(numout,*) '    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    531          WRITE(numout,*) '    N Supply   : ', nitdepinput*rno3*1E3/1E12*14.,' TgN/yr' 
    532          WRITE(numout,*)  
    533       ENDIF 
     502      ENDIF 
     503      ! 
     504      sedsilfrac = 0.03     ! percentage of silica loss in the sediments 
     505      sedcalfrac = 0.6      ! percentage of calcite loss in the sediments 
    534506      ! 
    535507   END SUBROUTINE p4z_sbc_init 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90

    r10068 r10127  
    4949      ! 
    5050      INTEGER, INTENT(in) ::   kt, knt ! ocean time step 
    51       INTEGER  ::   ji, jj, jk, ikt 
    52       REAL(wp) ::   zsumsedsi, zsumsedpo4, zsumsedcal 
    53       REAL(wp) ::   zrivalk, zrivsil, zrivno3 
     51      INTEGER  ::  ji, jj, jk, ikt 
     52      REAL(wp) ::  zrivalk, zrivsil, zrivno3 
    5453      REAL(wp) ::  zwflux, zfminus, zfplus 
    5554      REAL(wp) ::  zlim, zfact, zfactcal 
     
    6261      ! 
    6362      CHARACTER (len=25) :: charout 
    64       REAL(wp), DIMENSION(jpi,jpj    ) :: zwork1, zwork2, zwork3 
    65       REAL(wp), DIMENSION(jpi,jpj    ) :: zdenit2d, zbureff 
     63      REAL(wp), DIMENSION(jpi,jpj    ) :: zdenit2d, zbureff, zwork 
    6664      REAL(wp), DIMENSION(jpi,jpj    ) :: zwsbio3, zwsbio4, zwscal 
    6765      REAL(wp), DIMENSION(jpi,jpj    ) :: zsedcal, zsedsi, zsedc 
     
    8280      zdenit2d(:,:) = 0.e0 
    8381      zbureff (:,:) = 0.e0 
    84       zwork1  (:,:) = 0.e0 
    85       zwork2  (:,:) = 0.e0 
    86       zwork3  (:,:) = 0.e0 
     82      zwork   (:,:) = 0.e0 
    8783      zsedsi  (:,:) = 0.e0 
    8884      zsedcal (:,:) = 0.e0 
     
    255251                   &     + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) ) * 1E6 
    256252                 zbureff(ji,jj) = 0.013 + 0.53 * zflx**2 / ( 7.0 + zflx )**2 
    257                 ENDIF 
    258               END DO 
    259            END DO  
    260  
    261            ! Loss of biogenic silicon, Caco3 organic carbon in the sediments.  
    262            ! First, the total loss is computed. 
    263            ! The factor for calcite comes from the alkalinity effect 
    264            ! ------------------------------------------------------------- 
    265            DO jj = 1, jpj 
    266               DO ji = 1, jpi 
    267                  IF( tmask(ji,jj,1) == 1 ) THEN 
    268                     ikt = mbkt(ji,jj)  
    269                     zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwsbio4(ji,jj) 
    270                     zwork2(ji,jj) = trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj)  
    271                     ! For calcite, burial efficiency is made a function of saturation 
    272                     zfactcal      = MIN( excess(ji,jj,ikt), 0.2 ) 
    273                     zfactcal      = MIN( 1., 1.3 * ( 0.2 - zfactcal ) / ( 0.4 - zfactcal ) ) 
    274                     zwork3(ji,jj) = trb(ji,jj,ikt,jpcal) * zwscal(ji,jj) * 2.e0 * zfactcal 
    275                 ENDIF 
    276             END DO 
    277          END DO 
    278          zsumsedsi  = glob_sum( zwork1(:,:) * e1e2t(:,:) ) * r1_rday 
    279          zsumsedpo4 = glob_sum( zwork2(:,:) * e1e2t(:,:) ) * r1_rday 
    280          zsumsedcal = glob_sum( zwork3(:,:) * e1e2t(:,:) ) * r1_rday 
     253              ENDIF 
     254            END DO 
     255         END DO  
    281256         ! 
    282257      ENDIF 
     
    285260      ! Thus, the amount of silica lost in the sediments equal the supply at the surface (dust+rivers) 
    286261      ! ------------------------------------------------------ 
    287       IF( .NOT.lk_sed )  zrivsil =  1._wp - ( sumdepsi + rivdsiinput * r1_ryyss ) / ( zsumsedsi + rtrn ) 
     262      IF( .NOT.lk_sed )  zrivsil = 1._wp - sedsilfrac 
    288263 
    289264      DO jj = 1, jpj 
     
    312287               zfactcal = MIN( excess(ji,jj,ikt), 0.2 ) 
    313288               zfactcal = MIN( 1., 1.3 * ( 0.2 - zfactcal ) / ( 0.4 - zfactcal ) ) 
    314                zrivalk  =  1._wp - ( rivalkinput * r1_ryyss ) * zfactcal / ( zsumsedcal + rtrn ) 
     289               zrivalk  = sedcalfrac * zfactcal 
    315290               tra(ji,jj,ikt,jptal) =  tra(ji,jj,ikt,jptal) + zcaloss * zrivalk * 2.0 
    316291               tra(ji,jj,ikt,jpdic) =  tra(ji,jj,ikt,jpdic) + zcaloss * zrivalk 
     
    492467            IF( iom_use("Nfix"   ) ) CALL iom_put( "Nfix", nitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) )  ! nitrogen fixation  
    493468            IF( iom_use("INTNFIX") ) THEN   ! nitrogen fixation rate in ocean ( vertically integrated ) 
    494                zwork1(:,:) = 0. 
     469               zwork(:,:) = 0. 
    495470               DO jk = 1, jpkm1 
    496                  zwork1(:,:) = zwork1(:,:) + nitrpot(:,:,jk) * nitrfix * rno3 * zfact * e3t_n(:,:,jk) * tmask(:,:,jk) 
     471                 zwork(:,:) = zwork(:,:) + nitrpot(:,:,jk) * nitrfix * rno3 * zfact * e3t_n(:,:,jk) * tmask(:,:,jk) 
    497472               ENDDO 
    498                CALL iom_put( "INTNFIX" , zwork1 )  
     473               CALL iom_put( "INTNFIX" , zwork )  
    499474            ENDIF 
    500475            IF( iom_use("SedCal" ) ) CALL iom_put( "SedCal", zsedcal(:,:) * zfact ) 
Note: See TracChangeset for help on using the changeset viewer.