Changeset 9374
- Timestamp:
- 2018-03-06T13:35:36+01:00 (7 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM
- Files:
-
- 4 edited
-
CONFIG/SHARED/field_def_bgc.xml (modified) (1 diff)
-
NEMO/TOP_SRC/CFC/trcini_cfc.F90 (modified) (4 diffs)
-
NEMO/TOP_SRC/MEDUSA/air_sea.F90 (modified) (3 diffs)
-
NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/CONFIG/SHARED/field_def_bgc.xml
r9311 r9374 86 86 <field id="PHN_E3T" long_name="PHN * E3T" unit="mmol-N/m2" > PHN * e3t </field > 87 87 <field id="PHN_INV" field_ref="PHN_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 88 <field id="SURF_PHN" field_ref="PHN" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />88 <field id="SURF_PHN" field_ref="PHN" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 89 89 <field id="PHD" long_name="diatom phytoplankton" unit="mmol-N/m3" /> 90 90 <field id="PHD_E3T" long_name="PHD * E3T" unit="mmol-N/m2" > PHD * e3t </field > 91 91 <field id="PHD_INV" field_ref="PHD_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 92 <field id="SURF_PHD" field_ref="PHD" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />92 <field id="SURF_PHD" field_ref="PHD" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 93 93 <field id="ZMI" long_name="micro zooplankton" unit="mmol-N/m3" /> 94 94 <field id="ZMI_E3T" long_name="ZMI * E3T" unit="mmol-N/m2" > ZMI * e3t </field > 95 95 <field id="ZMI_INV" field_ref="ZMI_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 96 <field id="SURF_ZMI" field_ref="ZMI" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />96 <field id="SURF_ZMI" field_ref="ZMI" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 97 97 <field id="ZME" long_name="meso zooplankton" unit="mmol-N/m3" /> 98 98 <field id="ZME_E3T" long_name="ZME * E3T" unit="mmol-N/m2" > ZME * e3t </field > 99 99 <field id="ZME_INV" field_ref="ZME_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 100 <field id="SURF_ZME" field_ref="ZME" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />100 <field id="SURF_ZME" field_ref="ZME" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 101 101 <field id="DIN" long_name="dissolved inorganic nitrogen" unit="mmol-N/m3" /> 102 102 <field id="DIN_E3T" long_name="DIN * E3T" unit="mmol-N/m2" > DIN * e3t </field > 103 103 <field id="DIN_INV" field_ref="DIN_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 104 <field id="SURF_DIN" field_ref="DIN" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />104 <field id="SURF_DIN" field_ref="DIN" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 105 105 <field id="SIL" long_name="dissolved silicic acid" unit="mmolSi/m3" /> 106 106 <field id="SIL_E3T" long_name="SIL * E3T" unit="mmolSi/m2" > SIL * e3t </field > 107 107 <field id="SIL_INV" field_ref="SIL_E3T" unit="mmolSi/m2" grid_ref="vert_sum" /> 108 <field id="SURF_SIL" field_ref="SIL" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />108 <field id="SURF_SIL" field_ref="SIL" unit="mmolSi/m3" grid_ref="grid_T_surface_extract" /> 109 109 <field id="FER" long_name="dissolved iron" unit="mmolFe/m3" /> 110 110 <field id="FER_E3T" long_name="FER * E3T" unit="mmolFe/m2" > FER * e3t </field > 111 111 <field id="FER_INV" field_ref="FER_E3T" unit="mmolFe/m2" grid_ref="vert_sum" /> 112 <field id="SURF_FER" field_ref="FER" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />112 <field id="SURF_FER" field_ref="FER" unit="mmolFe/m3" grid_ref="grid_T_surface_extract" /> 113 113 <field id="DET" long_name="detrital nitrogen" unit="mmol-N/m3" /> 114 114 <field id="DET_E3T" long_name="DET * E3T" unit="mmol-N/m2" > DET * e3t </field > 115 115 <field id="DET_INV" field_ref="DET_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 116 <field id="SURF_DET" field_ref="DET" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />116 <field id="SURF_DET" field_ref="DET" unit="mmol-N/m3" grid_ref="grid_T_surface_extract" /> 117 117 <field id="PDS" long_name="biogenic silicon in diatom phytoplankton" unit="mmolSi/m3" /> 118 118 <field id="PDS_E3T" long_name="PDS * E3T" unit="mmolSi/m2" > PDS * e3t </field > 119 119 <field id="PDS_INV" field_ref="PDS_E3T" unit="mmolSi/m2" grid_ref="vert_sum" /> 120 <field id="SURF_PDS" field_ref="PDS" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />120 <field id="SURF_PDS" field_ref="PDS" unit="mmolSi/m3" grid_ref="grid_T_surface_extract" /> 121 121 <field id="DTC" long_name="detrital carbon" unit="mmol-C/m3" /> 122 122 <field id="DTC_E3T" long_name="DTC * E3T" unit="mmol-C/m2" > DTC * e3t </field > 123 123 <field id="DTC_INV" field_ref="DTC_E3T" unit="mmol-C/m2" grid_ref="vert_sum" /> 124 <field id="SURF_DTC" field_ref="DTC" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />124 <field id="SURF_DTC" field_ref="DTC" unit="mmol-C/m3" grid_ref="grid_T_surface_extract" /> 125 125 <field id="DiC" long_name="dissolved inorganic carbon" unit="mmol-C/m3" /> 126 126 <field id="DiC_E3T" long_name="DiC * E3T" unit="mmol-C/m2" > DiC * e3t </field > 127 127 <field id="DiC_INV" field_ref="DiC_E3T" unit="mmol-C/m2" grid_ref="vert_sum" /> 128 <field id="SURF_DiC" field_ref="DiC" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />128 <field id="SURF_DiC" field_ref="DiC" unit="mmol-C/m3" grid_ref="grid_T_surface_extract" /> 129 129 <field id="ALK" long_name="alkalinity" unit="meq/m3" /> 130 130 <field id="ALK_E3T" long_name="ALK * E3T" unit="meq/m2" > ALK * e3t </field > 131 131 <field id="ALK_INV" field_ref="ALK_E3T" unit="meq/m2" grid_ref="vert_sum" /> 132 <field id="SURF_ALK" field_ref="ALK" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />132 <field id="SURF_ALK" field_ref="ALK" unit="meq/m3" grid_ref="grid_T_surface_extract" /> 133 133 <field id="OXY" long_name="dissolved oxygen" unit="mmolO2/m3" /> 134 134 <field id="OXY_E3T" long_name="OXY * E3T" unit="mmolO2/m2" > OXY * e3t </field > 135 135 <field id="OXY_INV" field_ref="OXY_E3T" unit="mmolO2/m2" grid_ref="vert_sum" /> 136 <field id="SURF_OXY" field_ref="OXY" unit="m g Chl/m3" grid_ref="grid_T_surface_extract" />136 <field id="SURF_OXY" field_ref="OXY" unit="mmolO2/m3" grid_ref="grid_T_surface_extract" /> 137 137 138 138 <!-- Ideal traceurs : variables available with key_idtra --> -
branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r8280 r9374 44 44 !! ** Method : - Read the namcfc namelist and check the parameter values 45 45 !!---------------------------------------------------------------------- 46 INTEGER :: ji, jj, jn, jl, jm, js, io, i err46 INTEGER :: ji, jj, jn, jl, jm, js, io, iostatus, ierr 47 47 INTEGER :: iskip = 7 ! number of 1st descriptor lines 48 48 REAL(wp) :: zyy, zyd … … 57 57 58 58 CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 59 REWIND(inum)60 59 61 60 ! compute the number of year in the file 62 61 ! file starts in 1931 do jn represent the year in the century 62 iostatus = 0 63 63 jn = 31 64 DO 65 READ(inum,'(1x)', END=100)64 DO WHILE ( iostatus == 0 ) 65 READ(inum,'(1x)', IOSTAT=iostatus, END=100) 66 66 jn = jn + 1 67 END DO 67 ENDDO 68 IF( iostatus .NE. 0 ) THEN 69 !! Error while reading CFC input file 70 CALL ctl_stop('trc_ini_cfc: & 71 & Error on the 1st reading of cfc1112sf6.atm') 72 RETURN 73 ENDIF 68 74 100 jpyear = jn - 1 - iskip 69 75 IF ( lwp) WRITE(numout,*) ' ', jpyear ,' years read' … … 72 78 ALLOCATE( p_cfc(jpyear,jphem,3), STAT=ierr ) 73 79 IF( ierr > 0 ) THEN 74 CALL ctl_stop( 'trc_ini_cfc: unable to allocate p_cfc array' ) ; RETURN 80 CALL ctl_stop( 'trc_ini_cfc: unable to allocate p_cfc array' ) 81 RETURN 75 82 ENDIF 76 83 IF( trc_sms_cfc_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'trc_ini_cfc: unable to allocate CFC arrays' ) … … 101 108 ! file starts in 1931 do jn represent the year in the century.jhh 102 109 ! Read file till the end 103 jn = 31 104 DO 105 READ(inum,*, IOSTAT=io) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), & 110 DO jn = 31, jpyear 111 READ(inum, '( F6.1, 6F7.2)', IOSTAT=io) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), & 106 112 & p_cfc(jn,1,3), p_cfc(jn,2,1), & 107 113 & p_cfc(jn,2,2), p_cfc(jn,2,3) 108 IF( io < 0 ) exit 109 jn = jn + 1 114 IF( io .NE.0 ) THEN 115 !! Error while reading CFC input file 116 CALL ctl_stop('trc_ini_cfc: & 117 & Error on the 2nd reading of cfc1112sf6.atm') 118 RETURN 119 ENDIF 110 120 END DO 111 121 -
branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
r9344 r9374 67 67 # endif 68 68 USE in_out_manager, ONLY: lwp, numout 69 USE oce, ONLY: PCO2a_in_cpl70 69 USE par_kind, ONLY: wp 71 70 USE par_oce, ONLY: jpi, jpim1, jpj, jpjm1 72 USE sbc_oce, ONLY: fr_i, lk_oasis,qsr, wndm71 USE sbc_oce, ONLY: fr_i, qsr, wndm 73 72 USE sms_medusa, ONLY: jdms, jdms_input, jdms_model, & 74 73 jriver_alk, jriver_c, & … … 152 151 !! OPEN wet point IF..THEN loop 153 152 IF (tmask(ji,jj,1) == 1) then 154 IF (lk_oasis) THEN155 !!! Jpalm test on atm xCO2156 IF ( (f_xco2a(ji,jj) .GT. 10000.0 ).OR.(f_xco2a(ji,jj) .LE. 0.0 ) ) THEN157 IF(lwp) THEN158 WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj), &159 ' -- ji =', mig(ji),' jj = ', mjg(jj)160 ENDIF161 CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic coupled atm xCO2 ' )162 ENDIF163 ENDIF164 153 !! 165 154 !! AXY (23/06/15): as part of an effort to update the … … 264 253 !! Enormous dpco2 - even if all carb chem caract are OK. 265 254 !! and this check stops the model. 266 !! ==> let's run the model without it and see how it goes. 255 !! --Input checks are already more than enough to stop the 256 !! model if carb chem goes crazy. 257 !! we remove the mocsy output checks 267 258 !!=================== 268 259 !!IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR. & -
branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r9327 r9374 303 303 WRITE(numout,*) '** MEDUSA Atm xCO2 given by the UM **' 304 304 f_xco2a(:,:) = PCO2a_in_cpl(:,:) 305 !! Check the xCO2 from the UM is OK 306 !! piece of code moved from air-sea.F90 307 !!--- 308 DO jj = 2,jpjm1 309 DO ji = 2,jpim1 310 !! OPEN wet point IF..THEN loop 311 IF (tmask(ji,jj,1) == 1) then 312 !!! Jpalm test on atm xCO2 313 IF ( (f_xco2a(ji,jj) .GT. 10000.0 ).OR. & 314 (f_xco2a(ji,jj) .LE. 0.0 ) ) THEN 315 IF(lwp) THEN 316 WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj), & 317 ' -- ji =', mig(ji),' jj = ', mjg(jj) 318 ENDIF 319 CALL ctl_stop( 'MEDUSA - Air-Sea :', & 320 'unrealistic coupled atm xCO2 ' ) 321 ENDIF 322 ENDIF 323 ENDDO 324 ENDDO 325 !!--- 305 326 ELSEIF (lk_pi_co2) THEN 306 327 !! OCMIP pre-industrial xCO2
Note: See TracChangeset
for help on using the changeset viewer.
