Changeset 12281 for NEMO/branches/2019
- Timestamp:
- 2019-12-21T11:47:12+01:00 (4 years ago)
- 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 90 90 ! Allocate temporary workspace 91 91 ! 92 zprorcan (:,:, jpk) = 0._wp ; zprorcad (:,:,jpk) = 0._wp ; zprofed (:,:,jpk) = 0._wp93 zprofen (:,:, jpk) = 0._wp ; zysopt (:,:,jpk) = 0._wp94 zpronewn (:,:, jpk) = 0._wp ; zpronewd (:,:,jpk) = 0._wp ; zprdia (:,:,jpk) = 0._wp95 zprbio (:,:, jpk) = 0._wp ; zprdch (:,:,jpk) = 0._wp ; zprnch (:,:,jpk) = 0._wp96 zmxl_fac (:,:, jpk) = 0._wp ; zmxl_chl (:,:,jpk) = 0._wp97 zpligprod1(:,:, jpk) = 0._wp ; zpligprod2(:,:,jpk) = 0._wp92 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 98 98 99 99 ! Computation of the optimal production … … 322 322 ! 323 323 IF( ln_ligand ) THEN 324 zpligprod1(:,:, jpk) = 0._wp ; zpligprod2(:,:,jpk) = 0._wp324 zpligprod1(:,:,:) = 0._wp ; zpligprod2(:,:,:) = 0._wp 325 325 DO jk = 1, jpkm1 326 326 DO jj = 1, jpj … … 353 353 CALL iom_put( "PFeN" , zprofen(:,:,:) * zfact * tmask(:,:,:) ) ! biogenic iron production by nanophyto 354 354 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 357 359 CALL iom_put( "Mumax" , zprmaxn(:,:,:) * tmask(:,:,:) ) ! Maximum growth rate 358 360 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 99 99 IF( ln_timing ) CALL timing_start('p5z_prod') 100 100 ! 101 zprorcan(:,:, jpk) = 0._wp ; zprorcap(:,:,jpk) = 0._wp ; zprorcad(:,:,jpk) = 0._wp102 zcroissn(:,:, jpk) = 0._wp ; zcroissp(:,:,jpk) = 0._wp ; zcroissd(:,:,jpk) = 0._wp103 zprofed (:,:, jpk) = 0._wp ; zprofep (:,:,jpk) = 0._wp ; zprofen (:,:,jpk) = 0._wp104 zpronewn(:,:, jpk) = 0._wp ; zpronewp(:,:,jpk) = 0._wp ; zpronewd(:,:,jpk) = 0._wp105 zproregn(:,:, jpk) = 0._wp ; zproregp(:,:,jpk) = 0._wp ; zproregd(:,:,jpk) = 0._wp106 zpropo4n(:,:, jpk) = 0._wp ; zpropo4p(:,:,jpk) = 0._wp ; zpropo4d(:,:,jpk) = 0._wp107 zprdia (:,:, jpk) = 0._wp ; zprpic (:,:,jpk) = 0._wp ; zprbio (:,:,jpk) = 0._wp108 zprodopn(:,:, jpk) = 0._wp ; zprodopp(:,:,jpk) = 0._wp ; zprodopd(:,:,jpk) = 0._wp109 zysopt (:,:, jpk) = 0._wp110 zrespn (:,:, jpk) = 0._wp ; zrespp (:,:,jpk) = 0._wp ; zrespd (:,:,jpk) = 0._wp101 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 111 111 112 112 ! Computation of the optimal production … … 444 444 ! 445 445 IF( ln_ligand ) THEN 446 zpligprod1(:,:, jpk) = 0._wp ; zpligprod2(:,:,jpk) = 0._wp446 zpligprod1(:,:,:) = 0._wp ; zpligprod2(:,:,:) = 0._wp 447 447 DO jk = 1, jpkm1 448 448 DO jj = 1, jpj … … 478 478 CALL iom_put( "PFeN" , zprofen(:,:,:) * zfact * tmask(:,:,:) ) ! biogenic iron production by nanophyto 479 479 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 482 484 CALL iom_put( "Mumax" , zprmaxn(:,:,:) * tmask(:,:,:) ) ! Maximum growth rate 483 485 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 47 47 IF( ln_timing ) CALL timing_start('trc_wri') 48 48 ! 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 ! 54 66 ENDIF 55 67 ! write the tracer concentrations in the file
Note: See TracChangeset
for help on using the changeset viewer.