- Timestamp:
- 2021-03-04T16:29:51+01:00 (3 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90
r13725 r14586 318 318 REAL(wp) :: rpig !: chla/chla+phea ratio (NAMELIST) 319 319 LOGICAL :: ln_rgb !: use RGB light scheme rather than 2-band (NAMELIST) 320 LOGICAL :: ln_chlbio !: use input chl for bio rgb 320 321 321 322 INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:) :: neln !: number of levels in the euphotic layer -
branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90
r13725 r14586 94 94 & xthetarem,xo2min 95 95 #endif 96 NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb 96 NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig,ln_rgb,ln_chlbio 97 97 INTEGER :: jl, jn 98 98 INTEGER :: ios ! Local integer output status for namelist read … … 1098 1098 rpig = 0. 1099 1099 ln_rgb = .false. 1100 ln_chlbio = .true. 1100 1101 1101 1102 !READ(numnatm,natopt) … … 1122 1123 WRITE(numout,*) ' chla/chla+phea ratio rpig = ',rpig 1123 1124 WRITE(numout,*) ' use RGB scheme ln_rgb = ',ln_rgb 1125 WRITE(numout,*) ' ln_chlbio = ',ln_chlbio 1124 1126 WRITE(numout,*) ' ' 1125 1127 -
branches/UKMO/dev_r5518_GO6_package_FOAMv14_readchl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcopt_medusa.F90
r14388 r14586 126 126 DO jj = 1, jpj 127 127 DO ji = 1, jpi 128 !totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd) 129 totchl = chl_for_qsr(ji,jj,jk) 128 IF ( ln_chlbio ) THEN 129 totchl = chl_for_qsr(ji,jj,jk) 130 ELSE 131 totchl = trn(ji,jj,jk,jpchn) + trn(ji,jj,jk,jpchd) 132 ENDIF 130 133 totchl = MIN( 10. , MAX( 0.05, totchl ) ) 131 134 irgb = NINT( 41 + 20.* LOG10( totchl ) + 1.e-15 )
Note: See TracChangeset
for help on using the changeset viewer.