Changeset 10261


Ignore:
Timestamp:
2018-10-31T16:30:10+01:00 (19 months ago)
Author:
dford
Message:

Ensure BGC diagnostics are properly land-masked.

Location:
branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/TOP_SRC/FABM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90

    r10241 r10261  
    151151            CALL write_25hourm_xml(xml_unit,model%horizontal_diagnostic_variables(jn)) 
    152152         END DO 
    153          WRITE (xml_unit,'(A)') '  <field id="visib" long_name="visibility" unit="m" default_value="-.200E+21" grid_ref="grid_T_3D" />' 
    154          WRITE (xml_unit,'(A)') '  <field id="visib25h" long_name="visibility 25-hour mean" unit="m" default_value="-.200E+21" grid_ref="grid_T_3D" />' 
    155          WRITE (xml_unit,'(A)') '  <field id="top_visib" long_name="Top-level visibility" unit="m" default_value="-.200E+21" grid_ref="grid_T_2D" />' 
    156          WRITE (xml_unit,'(A)') '  <field id="mid_visib" long_name="Middle-level visibility" unit="m" default_value="-.200E+21" grid_ref="grid_T_2D" />' 
    157          WRITE (xml_unit,'(A)') '  <field id="bot_visib" long_name="Bottom-level visibility" unit="m" default_value="-.200E+21" grid_ref="grid_T_2D" />' 
     153         WRITE (xml_unit,'(A)') '  <field id="visib" long_name="visibility" unit="m" default_value="1.e+20" grid_ref="grid_T_3D" />' 
     154         WRITE (xml_unit,'(A)') '  <field id="visib25h" long_name="visibility 25-hour mean" unit="m" default_value="1.e+20" grid_ref="grid_T_3D" />' 
     155         WRITE (xml_unit,'(A)') '  <field id="top_visib" long_name="Top-level visibility" unit="m" default_value="1.e+20" grid_ref="grid_T_2D" />' 
     156         WRITE (xml_unit,'(A)') '  <field id="mid_visib" long_name="Middle-level visibility" unit="m" default_value="1.e+20" grid_ref="grid_T_2D" />' 
     157         WRITE (xml_unit,'(A)') '  <field id="bot_visib" long_name="Bottom-level visibility" unit="m" default_value="1.e+20" grid_ref="grid_T_2D" />' 
    158158         WRITE (xml_unit,1000) ' </field_group>' 
    159159 
     
    238238      ENDIF 
    239239 
    240       WRITE (missing_value,'(E9.3)') variable%missing_value 
     240      WRITE (missing_value,'(E9.3)') 1.e+20 
    241241      WRITE (string_dimensions,'(I1)') number_dimensions 
    242242      SELECT CASE (number_dimensions) 
     
    261261      CHARACTER(LEN=20) :: missing_value 
    262262 
    263       WRITE (missing_value,'(E9.3)') variable%missing_value 
     263      WRITE (missing_value,'(E9.3)') 1.e+20 
    264264      WRITE (xml_unit,'(A)') '  <field id="top_'//TRIM(variable%name)//'" long_name="Top-level '//TRIM(variable%long_name)//'" unit="'//TRIM(variable%units)//'" default_value="'//TRIM(ADJUSTL(missing_value))//'" grid_ref="grid_T_2D"/>' 
    265265      WRITE (xml_unit,'(A)') '  <field id="mid_'//TRIM(variable%name)//'" long_name="Middle-level '//TRIM(variable%long_name)//'" unit="'//TRIM(variable%units)//'" default_value="'//TRIM(ADJUSTL(missing_value))//'" grid_ref="grid_T_2D"/>' 
  • branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/TOP_SRC/FABM/trcwri_fabm.F90

    r10241 r10261  
    8181      INTEGER, INTENT( in )               :: kt 
    8282      INTEGER              :: jn 
     83       
     84      REAL(wp), DIMENSION(jpi,jpj,jpk)    :: zw3d 
    8385 
    8486#if defined key_tracer_budget 
     
    105107             CALL iom_put( model%diagnostic_variables(jn)%name, fabm_get_bulk_diagnostic_data(model,jn)) 
    106108      END DO 
    107       CALL iom_put( 'visib', visib(:,:,:) ) 
     109      zw3d(:,:,:) = visib(:,:,:)*tmask(:,:,:) + 1.e+20*(1.0-tmask(:,:,:)) 
     110      CALL iom_put( 'visib', zw3d ) 
    108111 
    109112      ! write 2D diagnostics in the file 
Note: See TracChangeset for help on using the changeset viewer.