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 3475 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90 – NEMO

Ignore:
Timestamp:
2012-09-13T19:00:23+02:00 (12 years ago)
Author:
cetlod
Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB last updates from PISCES, hopefully... see ticket #972

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r3451 r3475  
    146146                     &              * 1.E3 / po4r / ( 31. * zcoef + rtrn ) 
    147147                  rivdsi(ji,jj) =   sf_river(jr_dsi)%fnow(ji,jj,1)                                    & 
    148                      &              * 1.E3        / ( 28. * zcoef + rtrn ) 
     148                     &              * 1.E3        / ( 28.1 * zcoef + rtrn ) 
    149149               END DO 
    150150            END DO 
    … …  
    325325         ALLOCATE( rivdic(jpi,jpj), rivalk(jpi,jpj), rivdin(jpi,jpj), rivdip(jpi,jpj), rivdsi(jpi,jpj) )  
    326326         ! 
    327          ALLOCATE( sf_river(jpriv), rivinput(jpriv), STAT=ierr1 )           !* allocate and fill sf_sst (forcing structure) with sn_sst 
     327         ALLOCATE( sf_river(jpriv), rivinput(jpriv), STAT=ierr1 )           !* allocate and fill sf_river (forcing structure) with sn_river_ 
    328328         rivinput(:) = 0.0 
    329329 
    … …  
    355355         rivdininput = (rivinput(jr_din) + rivinput(jr_don) ) * 1E3 / rno3 / 14._wp 
    356356         rivdipinput = (rivinput(jr_dip) + rivinput(jr_dop) ) * 1E3 / po4r / 31._wp 
    357          rivdsiinput = rivinput(jr_dsi) * 1E3 / 28._wp 
     357         rivdsiinput = rivinput(jr_dsi) * 1E3 / 28.1_wp 
    358358         rivalkinput = rivinput(jr_dic) * 1E3 / 12._wp 
    359359         ! 
    … …  
    494494         WRITE(numout,*) '    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    495495         WRITE(numout,*) '    N Supply   : ', rivdininput*rno3*1E3/1E12*14.,' TgN/yr' 
    496          WRITE(numout,*) '    Si Supply  : ', rivdsiinput*1E3/1E12*28.,' TgSi/yr' 
     496         WRITE(numout,*) '    Si Supply  : ', rivdsiinput*1E3/1E12*28.1,' TgSi/yr' 
    497497         WRITE(numout,*) '    P Supply   : ', rivdipinput*1E3*po4r/1E12*31.,' TgP/yr' 
    498498         WRITE(numout,*) '    Alk Supply : ', rivalkinput*1E3/1E12,' Teq/yr' 
Note: See TracChangeset for help on using the changeset viewer.