New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11892 for NEMO – NEMO

Changeset 11892 for NEMO


Ignore:
Timestamp:
2019-11-13T09:44:18+01:00 (4 years ago)
Author:
laurent
Message:

Improving test-case STATION_ASF, use of relevant heights and depth for PAPA station.

Location:
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/BDY/bdydta.F90

    r11831 r11892  
    615615            ENDIF 
    616616 
    617             IF( llneed ) THEN                                              ! dta_bdy(jbdy)%xxx will be needed 
     617            IF( llneed .AND. iszdim > 0 ) THEN                             ! dta_bdy(jbdy)%xxx will be needed 
    618618               !                                                           !   -> must be associated with an allocated target 
    619619               ALLOCATE( bf_alias(1)%fnow( iszdim, 1, ipk ) )              ! allocate the target 

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcdcy.F90

    r11672 r11892  
    9595         ENDIF 
    9696      END IF 
    97        
     97 
    9898      ! Setting parameters for each new day: 
    9999      CALL sbc_dcy_param() 
    … …  
    176176         rtmd(:,:) = MOD( (rtmd(:,:) + 1._wp) , 1._wp) 
    177177      ENDIF 
    178        
     178 
    179179      ! If this is a new day, we have to update the dawn, dusk and scaling function 
    180180      !---------------------- 

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

    r11861 r11892  
    2020&namusr_def    !   STATION_ASF user defined namelist 
    2121!----------------------------------------------------------------------- 
    22    rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     22   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     23   !                    ! [PAPA station = 1 m] 
    2324/ 
    2425!----------------------------------------------------------------------- 
    … …  
    6869&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station) 
    6970!----------------------------------------------------------------------- 
    70    rn_lat1d    =      50.1 !  Column latitude 
    71    rn_lon1d    =    215.1  !  Column longitude 
    72    ln_c1d_locpt =  .true.  ! Localization of 1D config in a grid (T) or independant point (F) 
     71   rn_lat1d    =      50.1 !  Column latitude    [PAPA station] 
     72   rn_lon1d    =    215.1  !  Column longitude   [PAPA station] 
    7373/ 
    7474!----------------------------------------------------------------------- 
    … …  
    124124   ln_COARE_3p6 = .true.    ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125125   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    126       rn_zqt      =  2.     !  Air temperature & humidity reference height (m) 
    127       rn_zu       = 10.     !  Wind vector reference height (m) 
     126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
     127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
    128128      ! 
    129129      ! Skin is ONLY available in ECMWF and COARE algorithms: 

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

    r11861 r11892  
    2020&namusr_def    !   STATION_ASF user defined namelist 
    2121!----------------------------------------------------------------------- 
    22    rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     22   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     23   !                    ! [PAPA station = 1 m] 
    2324/ 
    2425!----------------------------------------------------------------------- 
    … …  
    6869&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station) 
    6970!----------------------------------------------------------------------- 
    70    rn_lat1d    =      50.1 !  Column latitude 
    71    rn_lon1d    =    215.1  !  Column longitude 
    72    ln_c1d_locpt =  .true.  ! Localization of 1D config in a grid (T) or independant point (F) 
     71   rn_lat1d    =      50.1 !  Column latitude    [PAPA station] 
     72   rn_lon1d    =    215.1  !  Column longitude   [PAPA station] 
    7373/ 
    7474!----------------------------------------------------------------------- 
    … …  
    124124   ln_COARE_3p6 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125125   ln_ECMWF     = .true.    ! "ECMWF"     algorithm   (IFS cycle 31) 
    126       rn_zqt      =  2.     !  Air temperature & humidity reference height (m) 
    127       rn_zu       = 10.     !  Wind vector reference height (m) 
     126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
     127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
    128128      ! 
    129129      ! Skin is ONLY available in ECMWF and COARE algorithms: 

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

    r11861 r11892  
    2020&namusr_def    !   STATION_ASF user defined namelist 
    2121!----------------------------------------------------------------------- 
    22    rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     22   rn_dept1    =  1.0   ! depth (m) at which the prescribed SST is taken  == depth of first T point! (IMPORTANT for warm-layer impact, the deeper, the more pronounced!) 
     23   !                    ! [PAPA station = 1 m] 
    2324/ 
    2425!----------------------------------------------------------------------- 
    … …  
    6869&namc1d        !   1D configuration options / STATION_ASF               ("key_c1d" default: PAPA station) 
    6970!----------------------------------------------------------------------- 
    70    rn_lat1d    =      50.1 !  Column latitude 
    71    rn_lon1d    =    215.1  !  Column longitude 
    72    ln_c1d_locpt =  .true.  ! Localization of 1D config in a grid (T) or independant point (F) 
     71   rn_lat1d    =      50.1 !  Column latitude    [PAPA station] 
     72   rn_lon1d    =    215.1  !  Column longitude   [PAPA station] 
    7373/ 
    7474!----------------------------------------------------------------------- 
    … …  
    124124   ln_COARE_3p6 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    125125   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    126       rn_zqt      =  2.     !  Air temperature & humidity reference height (m) 
    127       rn_zu       = 10.     !  Wind vector reference height (m) 
     126      rn_zqt    =  2.5      !  Air temperature & humidity reference height (m) [PAPA station] 
     127      rn_zu     =  4.       !  Wind vector reference height (m)                [PAPA station] 
    128128      ! 
    129129      ! Skin is ONLY available in ECMWF and COARE algorithms: 

  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/README.md

    r11861 r11892  
    88```STATION_ASF``` is a demonstration case that mimics an in-situ station (buoy, platform) dedicated to the estimation of surface air-sea fluxes by means of the measurement of traditional meteorological surface parameters. 
    99 
    10 ```STATION_ASF``` is based on the merging of the "single column" and the "standalone surface module" configurations of NEMO (C1D meets SAS). As such, the spatial domain of ```STATION_ASF``` is punctual (1D). 
     10```STATION_ASF``` is based on the merging of the "single column" and the "standalone surface module" configurations of NEMO. In short, it coulb defined as "SAS meets C1D". As such, the spatial domain of ```STATION_ASF``` is punctual (1D, well actually 3 x 3 as in C1D). 
    1111 
    12 ```STATION_ASF``` is therefore a very light and versatile tool to test the different bulk algorithms and cool-skin/warm-layer parameterizations included in NEMO. 
     12```STATION_ASF``` is therefore a versatile tool, and extremely light in terms of computing requirements, to test the different bulk algorithms and cool-skin/warm-layer parameterization options included in NEMO. 
    1313 
    14 As input ```STATION_ASF``` will require the sea surface parameters: 
     14As input ```STATION_ASF``` will require the traditional *bulk* sea surface parameters: 
    1515 
    16 - bulk sea surface temperature (SST) at $z_{SST}$ meters below the surface 
     16- sea surface temperature (SST) at $z_{SST}$ meters below the surface 
    1717- Surface current vector 
    1818- Sea surface salinity 
    … …  
    3131## Physical description 
    3232 
    33  
    34  
    3533### Important namelist parameters speficic to STATION_ASF 
    3634 
    … …  
    4745## Input files to test STATION ASF 
    4846 
    49 Three full years of hourly data from the PAPA buoy can be downloaded here: 
     47Three full years of processed hourly data from the PAPA station (buoy) can be downloaded here: 
    5048https://drive.google.com/file/d/1MxNvjhRHmMrL54y6RX7WIaM9-LGl--ZP/ 
    5149 
    52 These two files are everything you need to play with the set of namelists provided for this test-case. 
     50These three files are everything you need to play with the set of namelists provided for this test-case. 
    5351 
    5452- ```Station_PAPA_50N-145W_atm_hourly.nc```  → contains hourly surface atmospheric state 
    … …  
    5654- ```Station_PAPA_50N-145W_oce_hourly.nc``` → contains hourly sea surface state 
    5755 
     56For station PAPA (50.1 N, 144.9 W), air temperature and humidity are measured at 2.5 m, the wind speed at 4 m, and the SST at 1 m below the surface, hence the following namelist parameters are given: 
     57 
     58- ```rn_dept1 =    1.  ``` (&namusr_def) 
     59- ```rn_lat1d =  50.1 ``` (&namc1d) 
     60- ```rn_lon1d = 215.1``` (&namc1d) 
     61- ```rn_zqt   =   2.5``` (&namsbc_blk) 
     62- ```rn_zu    =    4.``` (&namsbc_blk) 
     63 
     64 
     65 
     66 
     67 
     68 
     69 
     70 
     71 
Note: See TracChangeset for help on using the changeset viewer.