Changeset 12377 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zsink.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/TOP/PISCES/P4Z/p4zsink.F90
r12276 r12377 38 38 INTEGER :: ik100 39 39 40 !! * Substitutions 41 # include "do_loop_substitute.h90" 40 42 !!---------------------------------------------------------------------- 41 43 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 49 51 !!---------------------------------------------------------------------- 50 52 51 SUBROUTINE p4z_sink ( kt, knt )53 SUBROUTINE p4z_sink ( kt, knt, Kbb, Kmm, Krhs ) 52 54 !!--------------------------------------------------------------------- 53 55 !! *** ROUTINE p4z_sink *** … … 59 61 !!--------------------------------------------------------------------- 60 62 INTEGER, INTENT(in) :: kt, knt 63 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 61 64 INTEGER :: ji, jj, jk 62 65 CHARACTER (len=25) :: charout … … 77 80 ! by data and from the coagulation theory 78 81 ! ----------------------------------------------------------- 79 DO jk = 1, jpkm1 80 DO jj = 1, jpj 81 DO ji = 1,jpi 82 zmax = MAX( heup_01(ji,jj), hmld(ji,jj) ) 83 zfact = MAX( 0., gdepw_n(ji,jj,jk+1) - zmax ) / wsbio2scale 84 wsbio4(ji,jj,jk) = wsbio2 + MAX(0., ( wsbio2max - wsbio2 )) * zfact 85 END DO 86 END DO 87 END DO 82 DO_3D_11_11( 1, jpkm1 ) 83 zmax = MAX( heup_01(ji,jj), hmld(ji,jj) ) 84 zfact = MAX( 0., gdepw(ji,jj,jk+1,Kmm) - zmax ) / wsbio2scale 85 wsbio4(ji,jj,jk) = wsbio2 + MAX(0., ( wsbio2max - wsbio2 )) * zfact 86 END_3D 88 87 89 88 ! limit the values of the sinking speeds to avoid numerical instabilities … … 102 101 ! Compute the sedimentation term using p4zsink2 for all the sinking particles 103 102 ! ----------------------------------------------------- 104 CALL trc_sink( kt, wsbio3, sinking , jppoc, rfact2 )105 CALL trc_sink( kt, wsbio3, sinkfer , jpsfe, rfact2 )106 CALL trc_sink( kt, wsbio4, sinking2, jpgoc, rfact2 )107 CALL trc_sink( kt, wsbio4, sinkfer2, jpbfe, rfact2 )108 CALL trc_sink( kt, wsbio4, sinksil , jpgsi, rfact2 )109 CALL trc_sink( kt, wsbio4, sinkcal , jpcal, rfact2 )103 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinking , jppoc, rfact2 ) 104 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkfer , jpsfe, rfact2 ) 105 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2, jpgoc, rfact2 ) 106 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinkfer2, jpbfe, rfact2 ) 107 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinksil , jpgsi, rfact2 ) 108 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinkcal , jpcal, rfact2 ) 110 109 111 110 IF( ln_p5z ) THEN … … 117 116 ! Compute the sedimentation term using p4zsink2 for all the sinking particles 118 117 ! ----------------------------------------------------- 119 CALL trc_sink( kt, wsbio3, sinkingn , jppon, rfact2 )120 CALL trc_sink( kt, wsbio3, sinkingp , jppop, rfact2 )121 CALL trc_sink( kt, wsbio4, sinking2n, jpgon, rfact2 )122 CALL trc_sink( kt, wsbio4, sinking2p, jpgop, rfact2 )118 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkingn , jppon, rfact2 ) 119 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkingp , jppop, rfact2 ) 120 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2n, jpgon, rfact2 ) 121 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2p, jpgop, rfact2 ) 123 122 ENDIF 124 123 … … 142 141 ENDIF 143 142 ! 144 IF( ln_ctl) THEN ! print mean trends (used for debugging)143 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging) 145 144 WRITE(charout, FMT="('sink')") 146 145 CALL prt_ctl_trc_info(charout) 147 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)146 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 148 147 ENDIF 149 148 !
Note: See TracChangeset
for help on using the changeset viewer.