Changeset 6225 for branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
- Timestamp:
- 2016-01-08T10:35:19+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
r4521 r6225 27 27 PUBLIC trc_ini_pisces ! called by trcini.F90 module 28 28 29 30 # include "top_substitute.h90"31 29 !!---------------------------------------------------------------------- 32 30 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 71 69 USE p4zmort ! Mortality terms for phytoplankton 72 70 USE p4zlys ! Calcite saturation 71 USE p4zsed ! Sedimentation & burial 73 72 ! 74 73 REAL(wp), SAVE :: sco2 = 2.312e-3_wp 75 REAL(wp), SAVE :: alka0 = 2.42 3e-3_wp74 REAL(wp), SAVE :: alka0 = 2.426e-3_wp 76 75 REAL(wp), SAVE :: oxyg0 = 177.6e-6_wp 77 REAL(wp), SAVE :: po4 = 2.1 74e-6_wp76 REAL(wp), SAVE :: po4 = 2.165e-6_wp 78 77 REAL(wp), SAVE :: bioma0 = 1.000e-8_wp 79 REAL(wp), SAVE :: silic1 = 91. 65e-6_wp80 REAL(wp), SAVE :: no3 = 3 1.04e-6_wp * 7.625_wp78 REAL(wp), SAVE :: silic1 = 91.51e-6_wp 79 REAL(wp), SAVE :: no3 = 30.9e-6_wp * 7.625_wp 81 80 ! 82 81 INTEGER :: ji, jj, jk, ierr … … 97 96 ierr = ierr + p4z_rem_alloc() 98 97 ierr = ierr + p4z_flx_alloc() 98 ierr = ierr + p4z_sed_alloc() 99 99 ! 100 100 IF( lk_mpp ) CALL mpp_sum( ierr ) 101 101 IF( ierr /= 0 ) CALL ctl_stop( 'STOP', 'pisces_alloc: unable to allocate PISCES arrays' ) 102 102 ! 103 ryyss = nyear_len(1) * rday ! number of seconds per year 104 r1_ryyss = 1. / ryyss 105 ! 103 106 104 107 CALL p4z_sms_init ! Maint routine 105 108 ! ! Time-step 106 rfact = rdttrc(1) ! ---------107 rfactr = 1. / rfact108 rfact2 = rfact / FLOAT( nrdttrc )109 rfact2r = 1. / rfact2110 111 IF(lwp) WRITE(numout,*) ' Passive Tracer time step rfact = ', rfact, ' rdt = ', rdttra(1)112 IF(lwp) write(numout,*) ' PISCES Biology time step rfact2 = ', rfact2113 114 115 109 116 110 ! Set biological ratios … … 162 156 END IF 163 157 164 ! Time step duration for biology165 xstep = rfact2 / rday166 158 167 159 CALL p4z_sink_init ! vertical flux of particulate organic matter
Note: See TracChangeset
for help on using the changeset viewer.