- Timestamp:
- 2013-10-18T11:14:35+02:00 (11 years ago)
- Location:
- branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
r3295 r4075 61 61 !! ** Purpose : Initialisation of the PISCES biochemical model 62 62 !!---------------------------------------------------------------------- 63 ! 64 INTEGER :: ji, jj, jk 65 REAL(wp) :: zcaralk, zbicarb, zco3 66 REAL(wp) :: ztmas, ztmas1 67 !!---------------------------------------------------------------------- 63 68 64 IF(lwp) WRITE(numout,*) 69 65 IF(lwp) WRITE(numout,*) ' trc_ini_pisces : PISCES biochemical model initialisation' … … 134 130 ENDIF 135 131 136 IF( .NOT. ln_rsttr ) THEN137 ! Initialization of chemical variables of the carbon cycle138 ! --------------------------------------------------------139 DO jk = 1, jpk140 DO jj = 1, jpj141 DO ji = 1, jpi142 ztmas = tmask(ji,jj,jk)143 ztmas1 = 1. - tmask(ji,jj,jk)144 zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / ( 1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) ) )145 zco3 = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1146 zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk )147 hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1148 END DO149 END DO150 END DO151 !152 END IF153 154 132 ! Time step duration for biology 155 133 xstep = rfact2 / rday -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90
r3320 r4075 60 60 INTEGER, INTENT( in ) :: kt ! ocean time-step index 61 61 !! 62 INTEGER :: jnt, jn, jl 62 INTEGER :: ji, jj, jk, jn, jl, jnt 63 REAL(wp) :: zcaralk, zbicarb, zco3 64 REAL(wp) :: ztmas, ztmas1 63 65 CHARACTER (len=25) :: charout 64 66 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrdpis … … 66 68 ! 67 69 IF( nn_timing == 1 ) CALL timing_start('trc_sms_pisces') 70 ! 71 IF( kt == nittrc000 .AND. .NOT. ln_rsttr ) THEN 72 ! Initialization of chemical variables of the carbon cycle 73 ! -------------------------------------------------------- 74 DO jk = 1, jpk 75 DO jj = 1, jpj 76 DO ji = 1, jpi 77 ztmas = tmask(ji,jj,jk) 78 ztmas1 = 1. - tmask(ji,jj,jk) 79 zcaralk = trn(ji,jj,jk,jptal) - borat(ji,jj,jk) / ( 1. + 1.E-8 / ( rtrn + akb3(ji,jj,jk) ) ) 80 zco3 = ( zcaralk - trn(ji,jj,jk,jpdic) ) * ztmas + 0.5e-3 * ztmas1 81 zbicarb = ( 2. * trn(ji,jj,jk,jpdic) - zcaralk ) 82 hi(ji,jj,jk) = ( ak23(ji,jj,jk) * zbicarb / zco3 ) * ztmas + 1.e-9 * ztmas1 83 END DO 84 END DO 85 END DO 86 ! 87 END IF 88 68 89 ! 69 90 IF( ln_pisdmp .AND. MOD( kt - nn_dttrc, nn_pisdmp ) == 0 ) CALL trc_sms_pisces_dmp( kt ) ! Relaxation of some tracers
Note: See TracChangeset
for help on using the changeset viewer.