Ignore:
Timestamp:
2018-11-30T16:38:17+01:00 (3 years ago)
Author:
aumont
Message:

Various bug fixes and improvements

File:
1 edited

Legend:

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

    r10069 r10362  
    6262      CHARACTER (len=25) :: charout 
    6363      REAL(wp), DIMENSION(jpi,jpj  )   :: totprod, totthick, totcons  
    64       REAL(wp), DIMENSION(jpi,jpj,jpk)   :: zremipoc, zremigoc, zorem3, ztremint 
     64      REAL(wp), DIMENSION(jpi,jpj,jpk)   :: zremipoc, zremigoc, zorem3, ztremint, zfolimi 
    6565      REAL(wp), DIMENSION(jpi,jpj,jpk,jcpoc) :: alphag 
    6666      !!--------------------------------------------------------------------- 
     
    9090      orem  (:,:,:)   = 0. 
    9191      ztremint(:,:,:) = 0. 
     92      zfolimi (:,:,:) = 0. 
    9293 
    9394      DO jn = 1, jcpoc 
     
    209210                  tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zorem2 
    210211                  tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer2 
     212                  zfolimi(ji,jj,jk)   = zofer2 
    211213               END DO 
    212214            END DO 
     
    239241                  tra(ji,jj,jk,jpgop) = tra(ji,jj,jk,jpgop) - zopop2 * (1. + solgoc) 
    240242                  tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zofer2 * (1. + solgoc) 
     243                  zfolimi(ji,jj,jk)   = zofer2 
    241244               END DO 
    242245            END DO 
     
    414417                    tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zorem 
    415418                    tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zofer 
     419                    zfolimi(ji,jj,jk)   = zfolimi(ji,jj,jk) + zofer 
    416420                  ENDIF 
    417421               END DO 
     
    439443                tra(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zopop  
    440444                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zofer  
     445                zfolimi(ji,jj,jk)   = zfolimi(ji,jj,jk) + zofer 
    441446             END DO 
    442447           END DO 
     
    449454          CALL iom_put( "REMINP" , zremipoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate 
    450455          CALL iom_put( "REMING" , zremigoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate 
     456          CALL iom_put( "REMINF" , zfolimi(:,:,:)  * tmask(:,:,:)  * 1.e+9 * zrfact2 )  ! Remineralisation rate 
    451457        ENDIF 
    452458     ENDIF 
Note: See TracChangeset for help on using the changeset viewer.