Changeset 12281


Ignore:
Timestamp:
2019-12-21T11:47:12+01:00 (5 weeks ago)
Author:
cetlod
Message:

bugfix in dev_r11943_MERGE_2019 : initialise local array in PISCES, see ticket #2357

Location:
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zprod.F90

    r12258 r12281  
    9090      !  Allocate temporary workspace 
    9191      ! 
    92       zprorcan  (:,:,jpk) = 0._wp ; zprorcad  (:,:,jpk) = 0._wp ; zprofed (:,:,jpk) = 0._wp 
    93       zprofen   (:,:,jpk) = 0._wp ; zysopt    (:,:,jpk) = 0._wp 
    94       zpronewn  (:,:,jpk) = 0._wp ; zpronewd  (:,:,jpk) = 0._wp ; zprdia  (:,:,jpk) = 0._wp 
    95       zprbio    (:,:,jpk) = 0._wp ; zprdch    (:,:,jpk) = 0._wp ; zprnch  (:,:,jpk) = 0._wp  
    96       zmxl_fac  (:,:,jpk) = 0._wp ; zmxl_chl  (:,:,jpk) = 0._wp  
    97       zpligprod1(:,:,jpk) = 0._wp ; zpligprod2(:,:,jpk) = 0._wp  
     92      zprorcan  (:,:,:) = 0._wp ; zprorcad  (:,:,:) = 0._wp ; zprofed (:,:,:) = 0._wp 
     93      zprofen   (:,:,:) = 0._wp ; zysopt    (:,:,:) = 0._wp 
     94      zpronewn  (:,:,:) = 0._wp ; zpronewd  (:,:,:) = 0._wp ; zprdia  (:,:,:) = 0._wp 
     95      zprbio    (:,:,:) = 0._wp ; zprdch    (:,:,:) = 0._wp ; zprnch  (:,:,:) = 0._wp  
     96      zmxl_fac  (:,:,:) = 0._wp ; zmxl_chl  (:,:,:) = 0._wp  
     97      zpligprod1(:,:,:) = 0._wp ; zpligprod2(:,:,:) = 0._wp  
    9898 
    9999      ! Computation of the optimal production 
     
    322322     ! 
    323323     IF( ln_ligand ) THEN 
    324          zpligprod1(:,:,jpk) = 0._wp    ;    zpligprod2(:,:,jpk) = 0._wp 
     324         zpligprod1(:,:,:) = 0._wp    ;    zpligprod2(:,:,:) = 0._wp 
    325325         DO jk = 1, jpkm1 
    326326            DO jj = 1, jpj 
     
    353353       CALL iom_put( "PFeN"    , zprofen(:,:,:) * zfact * tmask(:,:,:)  ) ! biogenic iron production by nanophyto 
    354354       CALL iom_put( "PFeD"    , zprofed(:,:,:) * zfact * tmask(:,:,:)  ) ! biogenic iron production by  diatomes 
    355        CALL iom_put( "LPRODP"  , zpligprod1(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
    356        CALL iom_put( "LDETP"   , zpligprod2(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     355       IF( ln_ligand ) THEN 
     356         CALL iom_put( "LPRODP"  , zpligprod1(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     357         CALL iom_put( "LDETP"   , zpligprod2(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     358       ENDIF 
    357359       CALL iom_put( "Mumax"   , zprmaxn(:,:,:) * tmask(:,:,:)  ) ! Maximum growth rate 
    358360       CALL iom_put( "MuN"     , zprbio(:,:,:) * xlimphy(:,:,:) * tmask(:,:,:) ) ! Realized growth rate for nanophyto 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p5zprod.F90

    r12258 r12281  
    9999      IF( ln_timing )   CALL timing_start('p5z_prod') 
    100100      ! 
    101       zprorcan(:,:,jpk) = 0._wp ; zprorcap(:,:,jpk) = 0._wp ; zprorcad(:,:,jpk) = 0._wp 
    102       zcroissn(:,:,jpk) = 0._wp ; zcroissp(:,:,jpk) = 0._wp ; zcroissd(:,:,jpk) = 0._wp 
    103       zprofed (:,:,jpk) = 0._wp ; zprofep (:,:,jpk) = 0._wp ; zprofen (:,:,jpk) = 0._wp 
    104       zpronewn(:,:,jpk) = 0._wp ; zpronewp(:,:,jpk) = 0._wp ; zpronewd(:,:,jpk) = 0._wp 
    105       zproregn(:,:,jpk) = 0._wp ; zproregp(:,:,jpk) = 0._wp ; zproregd(:,:,jpk) = 0._wp  
    106       zpropo4n(:,:,jpk) = 0._wp ; zpropo4p(:,:,jpk) = 0._wp ; zpropo4d(:,:,jpk) = 0._wp 
    107       zprdia  (:,:,jpk) = 0._wp ; zprpic  (:,:,jpk) = 0._wp ; zprbio  (:,:,jpk) = 0._wp 
    108       zprodopn(:,:,jpk) = 0._wp ; zprodopp(:,:,jpk) = 0._wp ; zprodopd(:,:,jpk) = 0._wp 
    109       zysopt  (:,:,jpk) = 0._wp  
    110       zrespn  (:,:,jpk) = 0._wp ; zrespp  (:,:,jpk) = 0._wp ; zrespd  (:,:,jpk) = 0._wp  
     101      zprorcan(:,:,:) = 0._wp ; zprorcap(:,:,:) = 0._wp ; zprorcad(:,:,:) = 0._wp 
     102      zcroissn(:,:,:) = 0._wp ; zcroissp(:,:,:) = 0._wp ; zcroissd(:,:,:) = 0._wp 
     103      zprofed (:,:,:) = 0._wp ; zprofep (:,:,:) = 0._wp ; zprofen (:,:,:) = 0._wp 
     104      zpronewn(:,:,:) = 0._wp ; zpronewp(:,:,:) = 0._wp ; zpronewd(:,:,:) = 0._wp 
     105      zproregn(:,:,:) = 0._wp ; zproregp(:,:,:) = 0._wp ; zproregd(:,:,:) = 0._wp  
     106      zpropo4n(:,:,:) = 0._wp ; zpropo4p(:,:,:) = 0._wp ; zpropo4d(:,:,:) = 0._wp 
     107      zprdia  (:,:,:) = 0._wp ; zprpic  (:,:,:) = 0._wp ; zprbio  (:,:,:) = 0._wp 
     108      zprodopn(:,:,:) = 0._wp ; zprodopp(:,:,:) = 0._wp ; zprodopd(:,:,:) = 0._wp 
     109      zysopt  (:,:,:) = 0._wp  
     110      zrespn  (:,:,:) = 0._wp ; zrespp  (:,:,:) = 0._wp ; zrespd  (:,:,:) = 0._wp  
    111111 
    112112      ! Computation of the optimal production 
     
    444444     ! 
    445445     IF( ln_ligand ) THEN 
    446          zpligprod1(:,:,jpk) = 0._wp    ;    zpligprod2(:,:,jpk) = 0._wp              
     446         zpligprod1(:,:,:) = 0._wp    ;    zpligprod2(:,:,:) = 0._wp              
    447447         DO jk = 1, jpkm1 
    448448            DO jj = 1, jpj 
     
    478478       CALL iom_put( "PFeN"    , zprofen(:,:,:) * zfact * tmask(:,:,:)  ) ! biogenic iron production by nanophyto 
    479479       CALL iom_put( "PFeD"    , zprofed(:,:,:) * zfact * tmask(:,:,:)  ) ! biogenic iron production by  diatomes 
    480        CALL iom_put( "LPRODP"  , zpligprod1(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
    481        CALL iom_put( "LDETP"   , zpligprod2(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     480       IF( ln_ligand ) THEN 
     481         CALL iom_put( "LPRODP"  , zpligprod1(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     482         CALL iom_put( "LDETP"   , zpligprod2(:,:,:) * 1e9 * zfact * tmask(:,:,:) ) 
     483       ENDIF 
    482484       CALL iom_put( "Mumax"   , zprmaxn(:,:,:) * tmask(:,:,:)  ) ! Maximum growth rate 
    483485       CALL iom_put( "MuP"     , zprpic(:,:,:) * xlimpic(:,:,:) * tmask(:,:,:) ) ! Realized growth rate for picophyto 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/trcwri.F90

    r11949 r12281  
    4747      IF( ln_timing )   CALL timing_start('trc_wri') 
    4848      ! 
    49       IF( l_offline .AND. kt == nittrc000 .AND. lwp ) THEN    ! WRITE root name in date.file for use by postpro 
    50          CALL dia_nam( clhstnam, nn_writetrc,' ' ) 
    51          CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    52          WRITE(inum,*) clhstnam 
    53          CLOSE(inum) 
     49      IF( l_offline ) THEN    ! WRITE root name in date.file for use by postpro 
     50         IF(  kt == nittrc000 .AND. lwp ) THEN    ! WRITE root name in date.file for use by postpro 
     51           CALL dia_nam( clhstnam, nn_writetrc,' ' ) 
     52           CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     53           WRITE(inum,*) clhstnam 
     54           CLOSE(inum) 
     55        ENDIF 
     56 
     57       ! Output of initial vertical scale factor 
     58       CALL iom_put( "e3t_0", e3t_0(:,:,:) ) 
     59       CALL iom_put( "e3u_0", e3u_0(:,:,:) ) 
     60       CALL iom_put( "e3v_0", e3v_0(:,:,:) ) 
     61       ! 
     62       CALL iom_put( "e3t" , e3t(:,:,:,Kmm) ) 
     63       CALL iom_put( "e3u" , e3u(:,:,:,Kmm) ) 
     64       CALL iom_put( "e3v" , e3v(:,:,:,Kmm) ) 
     65       ! 
    5466      ENDIF 
    5567      ! write the tracer concentrations in the file 
Note: See TracChangeset for help on using the changeset viewer.