- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zsink.F90
r10425 r11949 49 49 !!---------------------------------------------------------------------- 50 50 51 SUBROUTINE p4z_sink ( kt, knt )51 SUBROUTINE p4z_sink ( kt, knt, Kbb, Kmm, Krhs ) 52 52 !!--------------------------------------------------------------------- 53 53 !! *** ROUTINE p4z_sink *** … … 59 59 !!--------------------------------------------------------------------- 60 60 INTEGER, INTENT(in) :: kt, knt 61 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 61 62 INTEGER :: ji, jj, jk 62 63 CHARACTER (len=25) :: charout … … 83 84 DO ji = 1,jpi 84 85 zmax = MAX( heup_01(ji,jj), hmld(ji,jj) ) 85 zfact = MAX( 0., gdepw _n(ji,jj,jk+1) - zmax ) / wsbio2scale86 zfact = MAX( 0., gdepw(ji,jj,jk+1,Kmm) - zmax ) / wsbio2scale 86 87 wsbio4(ji,jj,jk) = wsbio2 + MAX(0., ( wsbio2max - wsbio2 )) * zfact 87 88 END DO … … 104 105 ! Compute the sedimentation term using p4zsink2 for all the sinking particles 105 106 ! ----------------------------------------------------- 106 CALL trc_sink( kt, wsbio3, sinking , jppoc, rfact2 )107 CALL trc_sink( kt, wsbio3, sinkfer , jpsfe, rfact2 )108 CALL trc_sink( kt, wsbio4, sinking2, jpgoc, rfact2 )109 CALL trc_sink( kt, wsbio4, sinkfer2, jpbfe, rfact2 )110 CALL trc_sink( kt, wsbio4, sinksil , jpgsi, rfact2 )111 CALL trc_sink( kt, wsbio4, sinkcal , jpcal, rfact2 )107 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinking , jppoc, rfact2 ) 108 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkfer , jpsfe, rfact2 ) 109 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2, jpgoc, rfact2 ) 110 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinkfer2, jpbfe, rfact2 ) 111 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinksil , jpgsi, rfact2 ) 112 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinkcal , jpcal, rfact2 ) 112 113 113 114 IF( ln_p5z ) THEN … … 119 120 ! Compute the sedimentation term using p4zsink2 for all the sinking particles 120 121 ! ----------------------------------------------------- 121 CALL trc_sink( kt, wsbio3, sinkingn , jppon, rfact2 )122 CALL trc_sink( kt, wsbio3, sinkingp , jppop, rfact2 )123 CALL trc_sink( kt, wsbio4, sinking2n, jpgon, rfact2 )124 CALL trc_sink( kt, wsbio4, sinking2p, jpgop, rfact2 )122 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkingn , jppon, rfact2 ) 123 CALL trc_sink( kt, Kbb, Kmm, wsbio3, sinkingp , jppop, rfact2 ) 124 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2n, jpgon, rfact2 ) 125 CALL trc_sink( kt, Kbb, Kmm, wsbio4, sinking2p, jpgop, rfact2 ) 125 126 ENDIF 126 127 … … 175 176 WRITE(charout, FMT="('sink')") 176 177 CALL prt_ctl_trc_info(charout) 177 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)178 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 178 179 ENDIF 179 180 !
Note: See TracChangeset
for help on using the changeset viewer.