- Timestamp:
- 2019-11-22T15:29:17+01:00 (5 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/p4zlys.F90
r11536 r11949 43 43 CONTAINS 44 44 45 SUBROUTINE p4z_lys( kt, knt )45 SUBROUTINE p4z_lys( kt, knt, Kbb, Krhs ) 46 46 !!--------------------------------------------------------------------- 47 47 !! *** ROUTINE p4z_lys *** … … 54 54 !!--------------------------------------------------------------------- 55 55 INTEGER, INTENT(in) :: kt, knt ! ocean time step and ??? 56 INTEGER, INTENT(in) :: Kbb, Krhs ! time level indices 56 57 ! 57 58 INTEGER :: ji, jj, jk, jn … … 72 73 ! ------------------------------------------- 73 74 74 CALL solve_at_general( zhinit, zhi )75 CALL solve_at_general( zhinit, zhi, Kbb ) 75 76 76 77 DO jk = 1, jpkm1 77 78 DO jj = 1, jpj 78 79 DO ji = 1, jpi 79 zco3(ji,jj,jk) = tr b(ji,jj,jk,jpdic) * ak13(ji,jj,jk) * ak23(ji,jj,jk) / (zhi(ji,jj,jk)**2 &80 zco3(ji,jj,jk) = tr(ji,jj,jk,jpdic,Kbb) * ak13(ji,jj,jk) * ak23(ji,jj,jk) / (zhi(ji,jj,jk)**2 & 80 81 & + ak13(ji,jj,jk) * zhi(ji,jj,jk) + ak13(ji,jj,jk) * ak23(ji,jj,jk) + rtrn ) 81 82 hi (ji,jj,jk) = zhi(ji,jj,jk) * rhop(ji,jj,jk) / 1000. … … 109 110 ! (ACCORDING TO THIS FORMULATION ALSO SOME PARTICULATE 110 111 ! CACO3 GETS DISSOLVED EVEN IN THE CASE OF OVERSATURATION) 111 zdispot = kdca * zexcess * tr b(ji,jj,jk,jpcal)112 zdispot = kdca * zexcess * tr(ji,jj,jk,jpcal,Kbb) 112 113 ! CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3], 113 114 ! AND [SUM(CO2)] DUE TO CACO3 DISSOLUTION/PRECIPITATION 114 115 zcaldiss(ji,jj,jk) = zdispot * rfact2 / rmtss ! calcite dissolution 115 116 ! 116 tr a(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2. * zcaldiss(ji,jj,jk)117 tr a(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zcaldiss(ji,jj,jk)118 tr a(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zcaldiss(ji,jj,jk)117 tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + 2. * zcaldiss(ji,jj,jk) 118 tr(ji,jj,jk,jpcal,Krhs) = tr(ji,jj,jk,jpcal,Krhs) - zcaldiss(ji,jj,jk) 119 tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) + zcaldiss(ji,jj,jk) 119 120 END DO 120 121 END DO … … 132 133 WRITE(charout, FMT="('lys ')") 133 134 CALL prt_ctl_trc_info(charout) 134 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)135 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 135 136 ENDIF 136 137 !
Note: See TracChangeset
for help on using the changeset viewer.