Changeset 8607 for branches/NERC
- Timestamp:
- 2017-10-10T12:43:01+02:00 (7 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag_slice.F90
r8442 r8607 48 48 jdms, ocal_ccd, xpar, xze, & 49 49 zb_co2_flx, zb_dms_srf, & 50 zn_co2_flx, zn_dms_srf, zn_chl_srf 50 zn_co2_flx, zn_dms_srf, zn_chl_srf, & 51 scl_chl, chl_out 51 52 USE trc, ONLY: med_diag 52 53 … … 70 71 !! no need to output, just pass to cpl var 71 72 IF (lk_oasis) THEN 72 zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * 1.0E-6 !! surf Chl in Kg-chl/m3 as needed for cpl 73 IF (chl_out.eq.1) THEN 74 !! export and scale surface chl 75 zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 76 !! surf Chl in Kg-chl/m3 as needed for cpl 77 ELSEIF (chl_out.eq.2) THEN 78 !! export and scale mld chl 79 zn_chl_srf(:,:) = fchl_ml * scl_chl * 1.0E-6 80 !! mld Chl in Kg-chl/m3 as needed for cpl 81 ENDIF 73 82 chloro_out_cpl(:,:) = zn_chl_srf(:,:) !! Coupling Chl 74 83 END IF -
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
r8521 r8607 51 51 zb_sed_n, zb_sed_si, & 52 52 zn_sed_c, zn_sed_ca, zn_sed_fe, & 53 zn_sed_n, zn_sed_si 53 zn_sed_n, zn_sed_si, zn_chl_srf 54 54 USE trc, ONLY: med_diag, nittrc000 55 55 USE trcnam_trp, ONLY: ln_trcadv_cen2, ln_trcadv_tvd … … 241 241 # endif 242 242 IF ( med_diag%INVTN%dgsave ) THEN 243 CALL iom_put( "INVTN" , ftot_n ) 243 !! HACK - Coupled Chl sent to UM 244 CALL iom_put( "INVTN" , zn_chl_srf ) 245 !! CALL iom_put( "INVTN" , ftot_n ) 244 246 ENDIF 245 247 IF ( med_diag%INVTSI%dgsave ) THEN -
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90
r8132 r8607 177 177 INTEGER :: jdms_model !: choice of DMS model passed to atmosphere 178 178 !! 1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 179 !! 179 !! FOR UKESM 180 REAL(wp) :: scl_chl !: scaling factor for tuned Chl passed to the UM 181 INTEGER :: chl_out !: select Chl field exported and scaled for UM: 182 !: 1- Surface Chl ; 2- MLD Chl 180 183 !! 181 184 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: remdmp !: depth dependent damping coefficient of passive tracers -
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90
r8521 r8607 85 85 & jriver_n,jriver_si,jriver_c,jriver_alk,jriver_dep, & 86 86 & xsdiss, & 87 & sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model 87 & sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model, & 88 & scl_chl, chl_out 88 89 #if defined key_roam 89 90 NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit, & … … 246 247 jdms_input = 0 247 248 jdms_model = 0 249 scl_chl = 0.5 250 chl_out = 1 248 251 249 252 !REWIND(numnatm) … … 400 403 !! 1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM 401 404 !! 405 !! UKESM1 - exported Chl to UM 406 !! scl_chl : scaling factor to tune the chl field sent to the UM 407 !! chl_out : select the chl field to send at the UM: 408 !! 1- Surf Chl ; 2- MLD Chl 409 402 410 IF(lwp) THEN 403 411 !! … … 909 917 !! 910 918 !! UKESM1 - new diagnostics !! Jpalm; AXY (08/07/15) 911 WRITE(numout,*) '=== UKESM1-related parameters' 919 WRITE(numout,*) '=== UKESM1-related parameters ===' 920 WRITE(numout,*) ' ---- --- ---' 921 WRITE(numout,*) '=== UKESM1 -- coupled DMS to the UM' 912 922 WRITE(numout,*) & 913 923 & ' include DMS diagnostic?, jdms = ', jdms … … 935 945 & ' choice of DMS model passed to atmosphere, jdms_model = Anderson et al. (2001; modified)' 936 946 endif 947 WRITE(numout,*) '=== UKESM1 -- coupled Chl to the UM' 948 WRITE(numout,*) & 949 & ' Scaling factor to export tuned Chl to the UM scl_chl = ', scl_chl 950 IF (chl_out .eq. 1) THEN 951 WRITE(numout,*) & 952 & ' Chl field to be scaled and sent to the UM: chl_out = Surface Chl field ' 953 ELSEIF (chl_out .eq. 2) THEN 954 WRITE(numout,*) & 955 & ' Chl field to be scaled and sent to the UM: chl_out = MLD Chl field ' 956 ENDIF 937 957 !! 938 958 ENDIF -
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r8427 r8607 282 282 ELSE 283 283 IF(lwp) WRITE(numout,*) 'Chl surf concentration - setting to zero ...' 284 zn_chl_srf(:,:) = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6284 zn_chl_srf(:,:) = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 285 285 ENDIF 286 286 IF (lk_oasis) THEN
Note: See TracChangeset
for help on using the changeset viewer.