Changeset 3585
- Timestamp:
- 2012-11-16T18:14:12+01:00 (11 years ago)
- Location:
- branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM
- Files:
-
- 6 edited
-
CONFIG/AMM12/EXP00/namelist (modified) (1 diff)
-
CONFIG/GYRE/EXP00/namelist (modified) (1 diff)
-
CONFIG/ORCA2_LIM/EXP00/namelist (modified) (3 diffs)
-
CONFIG/ORCA2_OFF_PISCES/EXP00/namelist (modified) (1 diff)
-
NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (13 diffs)
-
NEMO/OPA_SRC/OBS/obs_oper.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/AMM12/EXP00/namelist
r3436 r3585 885 885 &namobs ! observation usage switch ('key_diaobs') 886 886 !----------------------------------------------------------------------- 887 ln_t3d = .false. ! Logical switch for T profile observations 888 ln_s3d = .false. ! Logical switch for S profile observations 889 ln_ena = .false. ! Logical switch for ENACT insitu data set 890 ! ! ln_cor Logical switch for Coriolis insitu data set 891 ln_profb = .false. ! Logical switch for feedback insitu data set 892 ln_sla = .false. ! Logical switch for SLA observations 893 894 ln_sladt = .false. ! Logical switch for AVISO SLA data 895 896 ln_slafb = .false. ! Logical switch for feedback SLA data 897 ! ln_ssh Logical switch for SSH observations 898 899 ln_sst = .true. ! Logical switch for SST observations 900 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 901 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 902 903 ln_sstfb = .false. ! Logical switch for feedback SST data 904 ! ln_sss Logical switch for SSS observations 905 ! ln_seaice Logical switch for Sea Ice observations 906 ! ln_vel3d Logical switch for velocity observations 907 ! ln_velavcur Logical switch for velocity daily av. cur. 908 ! ln_velhrcur Logical switch for velocity high freq. cur. 909 ! ln_velavadcp Logical switch for velocity daily av. ADCP 910 ! ln_velhradcp Logical switch for velocity high freq. ADCP 911 ! ln_velfb Logical switch for feedback velocity data 912 ! ln_grid_global Global distribtion of observations 913 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 914 ! grid_search_file Grid search lookup file header 915 ! enactfiles ENACT input observation file names 916 ! coriofiles Coriolis input observation file name 917 ! ! profbfiles: Profile feedback input observation file name 887 ln_t3d = .false. ! Logical switch for T profile observations 888 ln_s3d = .false. ! Logical switch for S profile observations 889 ln_ena = .false. ! Logical switch for ENACT insitu data set 890 ! ! ln_cor Logical switch for Coriolis insitu data set 891 ln_profb = .false. ! Logical switch for feedback insitu data set 892 ln_sla = .false. ! Logical switch for SLA observations 893 894 ln_sladt = .false. ! Logical switch for AVISO SLA data 895 896 ln_slafb = .false. ! Logical switch for feedback SLA data 897 ! ln_ssh Logical switch for SSH observations 898 899 ln_sst = .false. ! Logical switch for SST observations 900 ln_reysst = .false. ! Logical switch for Reynolds observations 901 ln_reysst_night = .false. ! Reynolds observations is a foundation SST 902 ln_ghrsst = .false. ! Logical switch for GHRSST observations 903 ln_ghrsst_night = .false. ! GHRSST is a foundation SST ln_reysst_night 904 ln_sstfb = .false. ! Logical switch for feedback SST data 905 ln_sstfb_night = .false. ! feedback is a foundation SST ln_reysst_night 906 907 ! ln_sss Logical switch for SSS observations 908 ! ln_seaice Logical switch for Sea Ice observations 909 ! ln_vel3d Logical switch for velocity observations 910 ! ln_velavcur Logical switch for velocity daily av. cur. 911 ! ln_velhrcur Logical switch for velocity high freq. cur. 912 ! ln_velavadcp Logical switch for velocity daily av. ADCP 913 ! ln_velhradcp Logical switch for velocity high freq. ADCP 914 ! ln_velfb Logical switch for feedback velocity data 915 ! ln_grid_global Global distribtion of observations 916 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 917 ! grid_search_file Grid search lookup file header 918 ! enactfiles ENACT input observation file names 919 ! coriofiles Coriolis input observation file name 920 ! ! profbfiles: Profile feedback input observation file name 918 921 profbfiles = 'profiles_01.nc' 919 ! ln_profb_enatim Enact feedback input time setting switch920 ! slafilesact Active SLA input observation file name921 ! slafilespas Passive SLA input observation file name922 ! ! slafbfiles: Feedback SLA input observation file name922 ! ln_profb_enatim Enact feedback input time setting switch 923 ! slafilesact Active SLA input observation file name 924 ! slafilespas Passive SLA input observation file name 925 ! ! slafbfiles: Feedback SLA input observation file name 923 926 slafbfiles = 'sla_01.nc' 924 ! sstfiles GHRSST input observation file name925 ! ! sstfbfiles: Feedback SST input observation file name927 ! sstfiles GHRSST input observation file name 928 ! ! sstfbfiles: Feedback SST input observation file name 926 929 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 927 ! seaicefiles Sea Ice input observation file name928 ! velavcurfiles Vel. cur. daily av. input file name929 ! velhvcurfiles Vel. cur. high freq. input file name930 ! velavadcpfiles Vel. ADCP daily av. input file name931 ! velhvadcpfiles Vel. ADCP high freq. input file name932 ! velfbfiles Vel. feedback input observation file name933 ! dobsini Initial date in window YYYYMMDD.HHMMSS934 ! dobsend Final date in window YYYYMMDD.HHMMSS935 ! n1dint Type of vertical interpolation method936 ! n2dint Type of horizontal interpolation method937 ! ln_nea Rejection of observations near land switch938 nmsshc = 0 ! MSSH correction scheme939 ! mdtcorr MDT correction940 ! mdtcutoff MDT cutoff for computed correction941 ln_altbias = .false. ! Logical switch for alt bias942 ln_ignmis = .true. ! Logical switch for ignoring missing files943 ! endailyavtypes ENACT daily average types930 ! seaicefiles Sea Ice input observation file name 931 ! velavcurfiles Vel. cur. daily av. input file name 932 ! velhvcurfiles Vel. cur. high freq. input file name 933 ! velavadcpfiles Vel. ADCP daily av. input file name 934 ! velhvadcpfiles Vel. ADCP high freq. input file name 935 ! velfbfiles Vel. feedback input observation file name 936 ! dobsini Initial date in window YYYYMMDD.HHMMSS 937 ! dobsend Final date in window YYYYMMDD.HHMMSS 938 ! n1dint Type of vertical interpolation method 939 ! n2dint Type of horizontal interpolation method 940 ! ln_nea Rejection of observations near land switch 941 nmsshc = 0 ! MSSH correction scheme 942 ! mdtcorr MDT correction 943 ! mdtcutoff MDT cutoff for computed correction 944 ln_altbias = .false. ! Logical switch for alt bias 945 ln_ignmis = .true. ! Logical switch for ignoring missing files 946 ! endailyavtypes ENACT daily average types 944 947 ln_grid_global = .true. 945 948 ln_grid_search_lookup = .false. 946 / 949 / 947 950 !----------------------------------------------------------------------- 948 951 &nam_asminc ! assimilation increments ('key_asminc') -
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/GYRE/EXP00/namelist
r3436 r3585 877 877 &namobs ! observation usage switch ('key_diaobs') 878 878 !----------------------------------------------------------------------- 879 ln_t3d = .false. ! Logical switch for T profile observations 880 ln_s3d = .false. ! Logical switch for S profile observations 881 ln_ena = .false. ! Logical switch for ENACT insitu data set 882 ! ! ln_cor Logical switch for Coriolis insitu data set 883 ln_profb = .false. ! Logical switch for feedback insitu data set 884 ln_sla = .false. ! Logical switch for SLA observations 885 886 ln_sladt = .false. ! Logical switch for AVISO SLA data 887 888 ln_slafb = .false. ! Logical switch for feedback SLA data 889 ! ln_ssh Logical switch for SSH observations 890 891 ln_sst = .true. ! Logical switch for SST observations 892 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 893 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 894 895 ln_sstfb = .false. ! Logical switch for feedback SST data 896 ! ln_sss Logical switch for SSS observations 897 ! ln_seaice Logical switch for Sea Ice observations 898 ! ln_vel3d Logical switch for velocity observations 899 ! ln_velavcur Logical switch for velocity daily av. cur. 900 ! ln_velhrcur Logical switch for velocity high freq. cur. 901 ! ln_velavadcp Logical switch for velocity daily av. ADCP 902 ! ln_velhradcp Logical switch for velocity high freq. ADCP 903 ! ln_velfb Logical switch for feedback velocity data 904 ! ln_grid_global Global distribtion of observations 905 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 906 ! grid_search_file Grid search lookup file header 907 ! enactfiles ENACT input observation file names 908 ! coriofiles Coriolis input observation file name 909 ! ! profbfiles: Profile feedback input observation file name 879 ln_t3d = .false. ! Logical switch for T profile observations 880 ln_s3d = .false. ! Logical switch for S profile observations 881 ln_ena = .false. ! Logical switch for ENACT insitu data set 882 ! ! ln_cor Logical switch for Coriolis insitu data set 883 ln_profb = .false. ! Logical switch for feedback insitu data set 884 ln_sla = .false. ! Logical switch for SLA observations 885 886 ln_sladt = .false. ! Logical switch for AVISO SLA data 887 888 ln_slafb = .false. ! Logical switch for feedback SLA data 889 ! ln_ssh Logical switch for SSH observations 890 891 ln_sst = .false. ! Logical switch for SST observations 892 ln_reysst = .false. ! Logical switch for Reynolds observations 893 ln_reysst_night = .false. ! Reynolds observations is a foundation SST 894 ln_ghrsst = .false. ! Logical switch for GHRSST observations 895 ln_ghrsst_night = .false. ! GHRSST is a foundation SST ln_reysst_night 896 ln_sstfb = .false. ! Logical switch for feedback SST data 897 ln_sstfb_night = .false. ! feedback is a foundation SST ln_reysst_night 898 899 ! ln_sss Logical switch for SSS observations 900 ! ln_seaice Logical switch for Sea Ice observations 901 ! ln_vel3d Logical switch for velocity observations 902 ! ln_velavcur Logical switch for velocity daily av. cur. 903 ! ln_velhrcur Logical switch for velocity high freq. cur. 904 ! ln_velavadcp Logical switch for velocity daily av. ADCP 905 ! ln_velhradcp Logical switch for velocity high freq. ADCP 906 ! ln_velfb Logical switch for feedback velocity data 907 ! ln_grid_global Global distribtion of observations 908 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 909 ! grid_search_file Grid search lookup file header 910 ! enactfiles ENACT input observation file names 911 ! coriofiles Coriolis input observation file name 912 ! ! profbfiles: Profile feedback input observation file name 910 913 profbfiles = 'profiles_01.nc' 911 ! ln_profb_enatim Enact feedback input time setting switch912 ! slafilesact Active SLA input observation file name913 ! slafilespas Passive SLA input observation file name914 ! ! slafbfiles: Feedback SLA input observation file name914 ! ln_profb_enatim Enact feedback input time setting switch 915 ! slafilesact Active SLA input observation file name 916 ! slafilespas Passive SLA input observation file name 917 ! ! slafbfiles: Feedback SLA input observation file name 915 918 slafbfiles = 'sla_01.nc' 916 ! sstfiles GHRSST input observation file name917 ! ! sstfbfiles: Feedback SST input observation file name919 ! sstfiles GHRSST input observation file name 920 ! ! sstfbfiles: Feedback SST input observation file name 918 921 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 919 ! seaicefiles Sea Ice input observation file name920 ! velavcurfiles Vel. cur. daily av. input file name921 ! velhvcurfiles Vel. cur. high freq. input file name922 ! velavadcpfiles Vel. ADCP daily av. input file name923 ! velhvadcpfiles Vel. ADCP high freq. input file name924 ! velfbfiles Vel. feedback input observation file name925 ! dobsini Initial date in window YYYYMMDD.HHMMSS926 ! dobsend Final date in window YYYYMMDD.HHMMSS927 ! n1dint Type of vertical interpolation method928 ! n2dint Type of horizontal interpolation method929 ! ln_nea Rejection of observations near land switch930 nmsshc = 0 ! MSSH correction scheme931 ! mdtcorr MDT correction932 ! mdtcutoff MDT cutoff for computed correction933 ln_altbias = .false. ! Logical switch for alt bias934 ln_ignmis = .true. ! Logical switch for ignoring missing files935 ! endailyavtypes ENACT daily average types922 ! seaicefiles Sea Ice input observation file name 923 ! velavcurfiles Vel. cur. daily av. input file name 924 ! velhvcurfiles Vel. cur. high freq. input file name 925 ! velavadcpfiles Vel. ADCP daily av. input file name 926 ! velhvadcpfiles Vel. ADCP high freq. input file name 927 ! velfbfiles Vel. feedback input observation file name 928 ! dobsini Initial date in window YYYYMMDD.HHMMSS 929 ! dobsend Final date in window YYYYMMDD.HHMMSS 930 ! n1dint Type of vertical interpolation method 931 ! n2dint Type of horizontal interpolation method 932 ! ln_nea Rejection of observations near land switch 933 nmsshc = 0 ! MSSH correction scheme 934 ! mdtcorr MDT correction 935 ! mdtcutoff MDT cutoff for computed correction 936 ln_altbias = .false. ! Logical switch for alt bias 937 ln_ignmis = .true. ! Logical switch for ignoring missing files 938 ! endailyavtypes ENACT daily average types 936 939 ln_grid_global = .true. 937 940 ln_grid_search_lookup = .false. 938 / 941 / 939 942 !----------------------------------------------------------------------- 940 943 &nam_asminc ! assimilation increments ('key_asminc') -
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist
r3436 r3585 27 27 cn_exp = "ORCA2" ! experience name 28 28 nn_it000 = 1 ! first time step 29 nn_itend = 5475 ! last time step (std 5475)29 nn_itend = 15 ! last time step (std 5475) 30 30 nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 31 31 nn_leapy = 0 ! Leap year calendar (1) or not (0) … … 769 769 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 770 770 ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold 771 jpni = 0! jpni number of processors following i (set automatically if < 1)772 jpnj = 0! jpnj number of processors following j (set automatically if < 1)773 jpnij = 0! jpnij number of local domains (set automatically if < 1)771 jpni = 2 ! jpni number of processors following i (set automatically if < 1) 772 jpnj = 4 ! jpnj number of processors following j (set automatically if < 1) 773 jpnij = 8 ! jpnij number of local domains (set automatically if < 1) 774 774 / 775 775 !----------------------------------------------------------------------- … … 881 881 &namobs ! observation usage switch ('key_diaobs') 882 882 !----------------------------------------------------------------------- 883 ln_t3d = .false. ! Logical switch for T profile observations 884 ln_s3d = .false. ! Logical switch for S profile observations 885 ln_ena = .false. ! Logical switch for ENACT insitu data set 886 ! ! ln_cor Logical switch for Coriolis insitu data set 887 ln_profb = .false. ! Logical switch for feedback insitu data set 888 ln_sla = .false. ! Logical switch for SLA observations 889 890 ln_sladt = .false. ! Logical switch for AVISO SLA data 891 892 ln_slafb = .false. ! Logical switch for feedback SLA data 893 ! ln_ssh Logical switch for SSH observations 894 895 ln_sst = .true. ! Logical switch for SST observations 896 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 897 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 898 899 ln_sstfb = .false. ! Logical switch for feedback SST data 900 ! ln_sss Logical switch for SSS observations 901 ! ln_seaice Logical switch for Sea Ice observations 902 ! ln_vel3d Logical switch for velocity observations 903 ! ln_velavcur Logical switch for velocity daily av. cur. 904 ! ln_velhrcur Logical switch for velocity high freq. cur. 905 ! ln_velavadcp Logical switch for velocity daily av. ADCP 906 ! ln_velhradcp Logical switch for velocity high freq. ADCP 907 ! ln_velfb Logical switch for feedback velocity data 908 ! ln_grid_global Global distribtion of observations 909 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 910 ! grid_search_file Grid search lookup file header 911 ! enactfiles ENACT input observation file names 912 ! coriofiles Coriolis input observation file name 913 ! ! profbfiles: Profile feedback input observation file name 883 ln_t3d = .false. ! Logical switch for T profile observations 884 ln_s3d = .false. ! Logical switch for S profile observations 885 ln_ena = .false. ! Logical switch for ENACT insitu data set 886 ! ! ln_cor Logical switch for Coriolis insitu data set 887 ln_profb = .false. ! Logical switch for feedback insitu data set 888 ln_sla = .false. ! Logical switch for SLA observations 889 890 ln_sladt = .false. ! Logical switch for AVISO SLA data 891 892 ln_slafb = .false. ! Logical switch for feedback SLA data 893 ! ln_ssh Logical switch for SSH observations 894 895 ln_sst = .true. ! Logical switch for SST observations 896 ln_reysst = .true. ! Logical switch for Reynolds observations 897 ln_reysst_night = .true. ! Reynolds observations is a foundation SST 898 ln_ghrsst = .false. ! Logical switch for GHRSST observations 899 ln_ghrsst_night = .false. ! GHRSST is a foundation SST ln_reysst_night 900 ln_sstfb = .false. ! Logical switch for feedback SST data 901 ln_sstfb_night = .false. ! feedback is a foundation SST ln_reysst_night 902 903 ! ln_sss Logical switch for SSS observations 904 ! ln_seaice Logical switch for Sea Ice observations 905 ! ln_vel3d Logical switch for velocity observations 906 ! ln_velavcur Logical switch for velocity daily av. cur. 907 ! ln_velhrcur Logical switch for velocity high freq. cur. 908 ! ln_velavadcp Logical switch for velocity daily av. ADCP 909 ! ln_velhradcp Logical switch for velocity high freq. ADCP 910 ! ln_velfb Logical switch for feedback velocity data 911 ! ln_grid_global Global distribtion of observations 912 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 913 ! grid_search_file Grid search lookup file header 914 ! enactfiles ENACT input observation file names 915 ! coriofiles Coriolis input observation file name 916 ! ! profbfiles: Profile feedback input observation file name 914 917 profbfiles = 'profiles_01.nc' 915 ! ln_profb_enatim Enact feedback input time setting switch916 ! slafilesact Active SLA input observation file name917 ! slafilespas Passive SLA input observation file name918 ! ! slafbfiles: Feedback SLA input observation file name918 ! ln_profb_enatim Enact feedback input time setting switch 919 ! slafilesact Active SLA input observation file name 920 ! slafilespas Passive SLA input observation file name 921 ! ! slafbfiles: Feedback SLA input observation file name 919 922 slafbfiles = 'sla_01.nc' 920 ! sstfiles GHRSST input observation file name921 ! ! sstfbfiles: Feedback SST input observation file name923 ! sstfiles GHRSST input observation file name 924 ! ! sstfbfiles: Feedback SST input observation file name 922 925 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 923 ! seaicefiles Sea Ice input observation file name924 ! velavcurfiles Vel. cur. daily av. input file name925 ! velhvcurfiles Vel. cur. high freq. input file name926 ! velavadcpfiles Vel. ADCP daily av. input file name927 ! velhvadcpfiles Vel. ADCP high freq. input file name928 ! velfbfiles Vel. feedback input observation file name929 ! dobsini Initial date in window YYYYMMDD.HHMMSS930 ! dobsend Final date in window YYYYMMDD.HHMMSS931 ! n1dint Type of vertical interpolation method932 ! n2dint Type of horizontal interpolation method933 ! ln_nea Rejection of observations near land switch934 nmsshc = 0 ! MSSH correction scheme935 ! mdtcorr MDT correction936 ! mdtcutoff MDT cutoff for computed correction937 ln_altbias = .false. ! Logical switch for alt bias938 ln_ignmis = .true. ! Logical switch for ignoring missing files939 ! endailyavtypes ENACT daily average types926 ! seaicefiles Sea Ice input observation file name 927 ! velavcurfiles Vel. cur. daily av. input file name 928 ! velhvcurfiles Vel. cur. high freq. input file name 929 ! velavadcpfiles Vel. ADCP daily av. input file name 930 ! velhvadcpfiles Vel. ADCP high freq. input file name 931 ! velfbfiles Vel. feedback input observation file name 932 ! dobsini Initial date in window YYYYMMDD.HHMMSS 933 ! dobsend Final date in window YYYYMMDD.HHMMSS 934 ! n1dint Type of vertical interpolation method 935 ! n2dint Type of horizontal interpolation method 936 ! ln_nea Rejection of observations near land switch 937 nmsshc = 0 ! MSSH correction scheme 938 ! mdtcorr MDT correction 939 ! mdtcutoff MDT cutoff for computed correction 940 ln_altbias = .false. ! Logical switch for alt bias 941 ln_ignmis = .true. ! Logical switch for ignoring missing files 942 ! endailyavtypes ENACT daily average types 940 943 ln_grid_global = .true. 941 944 ln_grid_search_lookup = .false. -
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist
r3436 r3585 878 878 &namobs ! observation usage switch ('key_diaobs') 879 879 !----------------------------------------------------------------------- 880 ln_t3d = .false. ! Logical switch for T profile observations 881 ln_s3d = .false. ! Logical switch for S profile observations 882 ln_ena = .false. ! Logical switch for ENACT insitu data set 883 ! ! ln_cor Logical switch for Coriolis insitu data set 884 ln_profb = .false. ! Logical switch for feedback insitu data set 885 ln_sla = .false. ! Logical switch for SLA observations 886 887 ln_sladt = .false. ! Logical switch for AVISO SLA data 888 889 ln_slafb = .false. ! Logical switch for feedback SLA data 890 ! ln_ssh Logical switch for SSH observations 891 892 ln_sst = .true. ! Logical switch for SST observations 893 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 894 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 895 896 ln_sstfb = .false. ! Logical switch for feedback SST data 897 ! ln_sss Logical switch for SSS observations 898 ! ln_seaice Logical switch for Sea Ice observations 899 ! ln_vel3d Logical switch for velocity observations 900 ! ln_velavcur Logical switch for velocity daily av. cur. 901 ! ln_velhrcur Logical switch for velocity high freq. cur. 902 ! ln_velavadcp Logical switch for velocity daily av. ADCP 903 ! ln_velhradcp Logical switch for velocity high freq. ADCP 904 ! ln_velfb Logical switch for feedback velocity data 905 ! ln_grid_global Global distribtion of observations 906 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 907 ! grid_search_file Grid search lookup file header 908 ! enactfiles ENACT input observation file names 909 ! coriofiles Coriolis input observation file name 910 ! ! profbfiles: Profile feedback input observation file name 880 ln_t3d = .false. ! Logical switch for T profile observations 881 ln_s3d = .false. ! Logical switch for S profile observations 882 ln_ena = .false. ! Logical switch for ENACT insitu data set 883 ! ! ln_cor Logical switch for Coriolis insitu data set 884 ln_profb = .false. ! Logical switch for feedback insitu data set 885 ln_sla = .false. ! Logical switch for SLA observations 886 887 ln_sladt = .false. ! Logical switch for AVISO SLA data 888 889 ln_slafb = .false. ! Logical switch for feedback SLA data 890 ! ln_ssh Logical switch for SSH observations 891 892 ln_sst = .false. ! Logical switch for SST observations 893 ln_reysst = .false. ! ln_reysst Logical switch for Reynolds observations 894 ln_reysst_night = .false. ! Reynolds observations is a foundation SST ln_reysst_night = ', ln_reysst_night 895 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 896 ln_ghrsst_night = .false. ! GHRSST is a foundation SST ln_reysst_night = ', ln_reysst_night 897 ln_sstfb = .false. ! Logical switch for feedback SST data 898 ln_sstfb_night = .false. ! feedback is a foundation SST ln_reysst_night = ', ln_reysst_night 899 900 ! ln_sss Logical switch for SSS observations 901 ! ln_seaice Logical switch for Sea Ice observations 902 ! ln_vel3d Logical switch for velocity observations 903 ! ln_velavcur Logical switch for velocity daily av. cur. 904 ! ln_velhrcur Logical switch for velocity high freq. cur. 905 ! ln_velavadcp Logical switch for velocity daily av. ADCP 906 ! ln_velhradcp Logical switch for velocity high freq. ADCP 907 ! ln_velfb Logical switch for feedback velocity data 908 ! ln_grid_global Global distribtion of observations 909 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 910 ! grid_search_file Grid search lookup file header 911 ! enactfiles ENACT input observation file names 912 ! coriofiles Coriolis input observation file name 913 ! ! profbfiles: Profile feedback input observation file name 911 914 profbfiles = 'profiles_01.nc' 912 ! ln_profb_enatim Enact feedback input time setting switch913 ! slafilesact Active SLA input observation file name914 ! slafilespas Passive SLA input observation file name915 ! ! slafbfiles: Feedback SLA input observation file name915 ! ln_profb_enatim Enact feedback input time setting switch 916 ! slafilesact Active SLA input observation file name 917 ! slafilespas Passive SLA input observation file name 918 ! ! slafbfiles: Feedback SLA input observation file name 916 919 slafbfiles = 'sla_01.nc' 917 ! sstfiles GHRSST input observation file name918 ! ! sstfbfiles: Feedback SST input observation file name920 ! sstfiles GHRSST input observation file name 921 ! ! sstfbfiles: Feedback SST input observation file name 919 922 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 920 ! seaicefiles Sea Ice input observation file name921 ! velavcurfiles Vel. cur. daily av. input file name922 ! velhvcurfiles Vel. cur. high freq. input file name923 ! velavadcpfiles Vel. ADCP daily av. input file name924 ! velhvadcpfiles Vel. ADCP high freq. input file name925 ! velfbfiles Vel. feedback input observation file name926 ! dobsini Initial date in window YYYYMMDD.HHMMSS927 ! dobsend Final date in window YYYYMMDD.HHMMSS928 ! n1dint Type of vertical interpolation method929 ! n2dint Type of horizontal interpolation method930 ! ln_nea Rejection of observations near land switch931 nmsshc = 0 ! MSSH correction scheme932 ! mdtcorr MDT correction933 ! mdtcutoff MDT cutoff for computed correction934 ln_altbias = .false. ! Logical switch for alt bias935 ln_ignmis = .true. ! Logical switch for ignoring missing files936 ! endailyavtypes ENACT daily average types923 ! seaicefiles Sea Ice input observation file name 924 ! velavcurfiles Vel. cur. daily av. input file name 925 ! velhvcurfiles Vel. cur. high freq. input file name 926 ! velavadcpfiles Vel. ADCP daily av. input file name 927 ! velhvadcpfiles Vel. ADCP high freq. input file name 928 ! velfbfiles Vel. feedback input observation file name 929 ! dobsini Initial date in window YYYYMMDD.HHMMSS 930 ! dobsend Final date in window YYYYMMDD.HHMMSS 931 ! n1dint Type of vertical interpolation method 932 ! n2dint Type of horizontal interpolation method 933 ! ln_nea Rejection of observations near land switch 934 nmsshc = 0 ! MSSH correction scheme 935 ! mdtcorr MDT correction 936 ! mdtcutoff MDT cutoff for computed correction 937 ln_altbias = .false. ! Logical switch for alt bias 938 ln_ignmis = .true. ! Logical switch for ignoring missing files 939 ! endailyavtypes ENACT daily average types 937 940 ln_grid_global = .true. 938 941 ln_grid_search_lookup = .false. -
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r3411 r3585 59 59 60 60 !! * Module variables 61 LOGICAL, PUBLIC :: ln_t3d !: Logical switch for temperature profiles 62 LOGICAL, PUBLIC :: ln_s3d !: Logical switch for salinity profiles 63 LOGICAL, PUBLIC :: ln_ena !: Logical switch for the ENACT data set 64 LOGICAL, PUBLIC :: ln_cor !: Logical switch for the Coriolis data set 65 LOGICAL, PUBLIC :: ln_profb !: Logical switch for profile feedback datafiles 66 LOGICAL, PUBLIC :: ln_sla !: Logical switch for sea level anomalies 67 LOGICAL, PUBLIC :: ln_sladt !: Logical switch for SLA from AVISO files 68 LOGICAL, PUBLIC :: ln_slafb !: Logical switch for SLA from feedback files 69 LOGICAL, PUBLIC :: ln_sst !: Logical switch for sea surface temperature 70 LOGICAL, PUBLIC :: ln_reysst !: Logical switch for Reynolds sea surface temperature 71 LOGICAL, PUBLIC :: ln_ghrsst !: Logical switch for GHRSST data 72 LOGICAL, PUBLIC :: ln_sstfb !: Logical switch for SST from feedback files 73 LOGICAL, PUBLIC :: ln_seaice !: Logical switch for sea ice concentration 74 LOGICAL, PUBLIC :: ln_vel3d !: Logical switch for velocity component (u,v) observations 75 LOGICAL, PUBLIC :: ln_velavcur !: Logical switch for raw daily averaged netCDF current meter vel. data 76 LOGICAL, PUBLIC :: ln_velhrcur !: Logical switch for raw high freq netCDF current meter vel. data 77 LOGICAL, PUBLIC :: ln_velavadcp !: Logical switch for raw daily averaged netCDF ADCP vel. data 78 LOGICAL, PUBLIC :: ln_velhradcp !: Logical switch for raw high freq netCDF ADCP vel. data 79 LOGICAL, PUBLIC :: ln_velfb !: Logical switch for velocities from feedback files 80 LOGICAL, PUBLIC :: ln_ssh !: Logical switch for sea surface height 81 LOGICAL, PUBLIC :: ln_sss !: Logical switch for sea surface salinity 82 LOGICAL, PUBLIC :: ln_nea !: Remove observations near land 83 LOGICAL, PUBLIC :: ln_altbias !: Logical switch for altimeter bias 84 LOGICAL, PUBLIC :: ln_ignmis !: Logical switch for ignoring missing files 85 LOGICAL, PUBLIC :: ln_s_at_t !: Logical switch to compute model S at T observations 61 LOGICAL, PUBLIC :: ln_t3d !: Logical switch for temperature profiles 62 LOGICAL, PUBLIC :: ln_s3d !: Logical switch for salinity profiles 63 LOGICAL, PUBLIC :: ln_ena !: Logical switch for the ENACT data set 64 LOGICAL, PUBLIC :: ln_cor !: Logical switch for the Coriolis data set 65 LOGICAL, PUBLIC :: ln_profb !: Logical switch for profile feedback datafiles 66 LOGICAL, PUBLIC :: ln_sla !: Logical switch for sea level anomalies 67 LOGICAL, PUBLIC :: ln_sladt !: Logical switch for SLA from AVISO files 68 LOGICAL, PUBLIC :: ln_slafb !: Logical switch for SLA from feedback files 69 LOGICAL, PUBLIC :: ln_sst !: Logical switch for sea surface temperature 70 LOGICAL, PUBLIC :: ln_reysst !: Logical switch for Reynolds sea surface temperature 71 LOGICAL, PUBLIC :: ln_reysst_night !: Reynolds SST is a foundation SST 72 LOGICAL, PUBLIC :: ln_ghrsst !: Logical switch for GHRSST data 73 LOGICAL, PUBLIC :: ln_ghrsst_night !: GHRSST observations is a foundation SST 74 LOGICAL, PUBLIC :: ln_sstfb !: Logical switch for SST from feedback files 75 LOGICAL, PUBLIC :: ln_sstfb_night !: SST from feedback is a foundation SST 76 LOGICAL, PUBLIC :: ln_seaice !: Logical switch for sea ice concentration 77 LOGICAL, PUBLIC :: ln_vel3d !: Logical switch for velocity component (u,v) observations 78 LOGICAL, PUBLIC :: ln_velavcur !: Logical switch for raw daily averaged netCDF current meter vel. data 79 LOGICAL, PUBLIC :: ln_velhrcur !: Logical switch for raw high freq netCDF current meter vel. data 80 LOGICAL, PUBLIC :: ln_velavadcp !: Logical switch for raw daily averaged netCDF ADCP vel. data 81 LOGICAL, PUBLIC :: ln_velhradcp !: Logical switch for raw high freq netCDF ADCP vel. data 82 LOGICAL, PUBLIC :: ln_velfb !: Logical switch for velocities from feedback files 83 LOGICAL, PUBLIC :: ln_ssh !: Logical switch for sea surface height 84 LOGICAL, PUBLIC :: ln_sss !: Logical switch for sea surface salinity 85 LOGICAL, PUBLIC :: ln_nea !: Remove observations near land 86 LOGICAL, PUBLIC :: ln_altbias !: Logical switch for altimeter bias 87 LOGICAL, PUBLIC :: ln_ignmis !: Logical switch for ignoring missing files 88 LOGICAL, PUBLIC :: ln_s_at_t !: Logical switch to compute model S at T observations 86 89 87 90 REAL(KIND=dp), PUBLIC :: dobsini !: Observation window start date YYYYMMDD.HHMMSS … … 107 110 & ld_velav !: Velocity data is daily averaged 108 111 LOGICAL, DIMENSION(:), ALLOCATABLE :: & 109 & l d_sstnight !: SST observation corresponds to night mean112 & l_sstnight !: SST observation corresponds to night mean 110 113 111 114 !!---------------------------------------------------------------------- … … 159 162 NAMELIST/namobs/ln_ena, ln_cor, ln_profb, ln_t3d, ln_s3d, & 160 163 & ln_sla, ln_sladt, ln_slafb, & 161 & ln_ssh, ln_sst, ln_ss tfb, ln_sss, ln_nea,&164 & ln_ssh, ln_sst, ln_sss, ln_nea, & 162 165 & enactfiles, coriofiles, profbfiles, & 163 166 & slafilesact, slafilespas, slafbfiles, & … … 166 169 & dobsini, dobsend, n1dint, n2dint, & 167 170 & nmsshc, mdtcorr, mdtcutoff, & 168 & ln_reysst, ln_ghrsst, reysstname, reysstfmt, & 171 & ln_reysst, ln_reysst_night, & 172 & ln_ghrsst, ln_ghrsst_night, & 173 & ln_sstfb , ln_sstfb_night, & 174 & reysstname, reysstfmt, & 169 175 & ln_grid_search_lookup, & 170 176 & grid_search_file, grid_search_res, & … … 212 218 ln_sst = .FALSE. 213 219 ln_seaice = .FALSE. 214 ln_reysst = .FALSE. 215 ln_ghrsst = .FALSE. 220 ln_reysst = .FALSE. 221 ln_reysst_night = .FALSE. 222 ln_ghrsst = .FALSE. 223 ln_ghrsst_night = .FALSE. 224 ln_sstfb = .FALSE. 225 ln_sstfb_night = .FALSE. 216 226 ln_sss = .FALSE. 217 227 ln_profb = .FALSE. … … 220 230 ln_sladt = .TRUE. 221 231 ln_slafb = .FALSE. 222 ln_sstfb = .FALSE.223 232 ln_velavcur = .FALSE. 224 233 ln_velhrcur = .FALSE. … … 355 364 WRITE(numout,*) ' Logical switch for SST observations ln_sst = ', ln_sst 356 365 WRITE(numout,*) ' Logical switch for Reynolds observations ln_reysst = ', ln_reysst 366 WRITE(numout,*) ' Reynolds observations is a foundation SST ln_reysst_night = ', ln_reysst_night 357 367 WRITE(numout,*) ' Logical switch for GHRSST observations ln_ghrsst = ', ln_ghrsst 368 WRITE(numout,*) ' GHRSST observations is a foundation SST ln_ghrsst_night = ', ln_ghrsst_night 358 369 WRITE(numout,*) ' Logical switch for feedback SST data ln_sstfb = ', ln_sstfb 370 WRITE(numout,*) ' Feedback SST data is a foundation SST ln_ghrsst_night = ', ln_sstfb_night 359 371 WRITE(numout,*) ' Logical switch for SSS observations ln_sss = ', ln_sss 360 372 WRITE(numout,*) ' Logical switch for Sea Ice observations ln_seaice = ', ln_seaice … … 721 733 722 734 ! - Sea surface temperature 723 IF 735 IF( ln_sst ) THEN 724 736 725 737 ! Set the number of variables for sst to 1 … … 739 751 ALLOCATE(sstdata(nsstsets)) 740 752 ALLOCATE(sstdatqc(nsstsets)) 741 ALLOCATE(l d_sstnight(nsstsets))742 sstdata(:)%nsurf =0743 sstdatqc(:)%nsurf =0744 l d_sstnight(:)=.false.753 ALLOCATE(l_sstnight(nsstsets)) 754 sstdata(:)%nsurf = 0 755 sstdatqc(:)%nsurf = 0 756 l_sstnight(:) = .false. 745 757 746 758 nsstsets = 0 … … 750 762 nsstsets = nsstsets + 1 751 763 752 l d_sstnight(nsstsets) = .TRUE.764 l_sstnight(nsstsets) = ln_reysst_night 753 765 754 766 CALL obs_rea_sst_rey( reysstname, reysstfmt, sstdata(nsstsets), & … … 758 770 & ln_nea ) 759 771 760 ENDIF772 ENDIF 761 773 762 IF (ln_ghrsst) THEN774 IF (ln_ghrsst) THEN 763 775 764 776 nsstsets = nsstsets + 1 765 777 766 l d_sstnight(nsstsets) = .TRUE.778 l_sstnight(nsstsets) = ln_ghrsst_night 767 779 768 780 CALL obs_rea_sst( 1, sstdata(nsstsets), jnumsst, & … … 783 795 nsstsets = nsstsets + 1 784 796 785 l d_sstnight(nsstsets) = .TRUE.797 l_sstnight(nsstsets) = ln_sstfb_night 786 798 787 799 CALL obs_rea_sst( 0, sstdata(nsstsets), 1, & … … 1105 1117 & kstp, jpi, jpj, nit000, idaystp, & 1106 1118 & tsn(:,:,1,jp_tem), tmask(:,:,1), & 1107 & n2dint, l d_sstnight(jsstset) )1119 & n2dint, l_sstnight(jsstset) ) 1108 1120 END DO 1109 1121 ENDIF -
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90
r3411 r3585 663 663 & psstn, & ! Model SST field 664 664 & psstmask ! Land-sea mask 665 LOGICAL, INTENT(IN) :: ld_nightav 665 666 666 667 !! * Local declarations … … 694 695 & igrdi, & 695 696 & igrdj 696 LOGICAL, INTENT(IN) :: ld_nightav697 697 698 698 !----------------------------------------------------------------------- … … 705 705 IF ( ld_nightav ) THEN 706 706 707 ! Initialize array for night mean 708 709 IF ( kt .EQ. 0 ) THEN 710 ALLOCATE ( icount_sstnight(kpi,kpj) ) 711 ALLOCATE ( imask_night(kpi,kpj) ) 712 ALLOCATE ( zintmp(kpi,kpj) ) 713 ALLOCATE ( zouttmp(kpi,kpj) ) 714 ALLOCATE ( zmeanday(kpi,kpj) ) 715 nday_qsr = -1 ! initialisation flag for nbc_dcy 716 ENDIF 717 718 ! Initialize daily mean for first timestep 719 idayend = MOD( kt - kit000 + 1, kdaystp ) 720 721 ! Added kt == 0 test to catch restart case 722 IF ( idayend == 1 .OR. kt == 0) THEN 723 IF (lwp) WRITE(numout,*) 'Reset sstdatqc%vdmean on time-step: ',kt 707 ! Initialize array for night mean 708 709 IF ( kt .EQ. 0 ) THEN 710 ALLOCATE ( icount_sstnight(kpi,kpj) ) 711 ALLOCATE ( imask_night(kpi,kpj) ) 712 ALLOCATE ( zintmp(kpi,kpj) ) 713 ALLOCATE ( zouttmp(kpi,kpj) ) 714 ALLOCATE ( zmeanday(kpi,kpj) ) 715 nday_qsr = -1 ! initialisation flag for nbc_dcy 716 ENDIF 717 718 ! Initialize daily mean for first timestep 719 idayend = MOD( kt - kit000 + 1, kdaystp ) 720 721 ! Added kt == 0 test to catch restart case 722 IF ( idayend == 1 .OR. kt == 0) THEN 723 IF (lwp) WRITE(numout,*) 'Reset sstdatqc%vdmean on time-step: ',kt 724 DO jj = 1, jpj 725 DO ji = 1, jpi 726 sstdatqc%vdmean(ji,jj) = 0.0 727 zmeanday(ji,jj) = 0.0 728 icount_sstnight(ji,jj) = 0 729 END DO 730 END DO 731 ENDIF 732 733 zintmp(:,:) = 0.0 734 zouttmp(:,:) = sbc_dcy( zintmp(:,:), .TRUE. ) 735 imask_night(:,:) = INT( zouttmp(:,:) ) 736 724 737 DO jj = 1, jpj 725 738 DO ji = 1, jpi 726 sstdatqc%vdmean(ji,jj) = 0.0 727 zmeanday(ji,jj) = 0.0 728 icount_sstnight(ji,jj) = 0 739 ! Increment the temperature field for computing night mean and counter 740 sstdatqc%vdmean(ji,jj) = sstdatqc%vdmean(ji,jj) & 741 & + psstn(ji,jj)*imask_night(ji,jj) 742 zmeanday(ji,jj) = zmeanday(ji,jj) + psstn(ji,jj) 743 icount_sstnight(ji,jj) = icount_sstnight(ji,jj) + imask_night(ji,jj) 729 744 END DO 730 745 END DO 731 ENDIF732 733 zintmp(:,:) = 0.0734 zouttmp(:,:) = sbc_dcy( zintmp(:,:), .TRUE. )735 imask_night(:,:) = INT( zouttmp(:,:) )736 737 DO jj = 1, jpj738 DO ji = 1, jpi739 ! Increment the temperature field for computing night mean and counter740 sstdatqc%vdmean(ji,jj) = sstdatqc%vdmean(ji,jj) &741 & + psstn(ji,jj)*imask_night(ji,jj)742 zmeanday(ji,jj) = zmeanday(ji,jj) + psstn(ji,jj)743 icount_sstnight(ji,jj) = icount_sstnight(ji,jj) + imask_night(ji,jj)744 END DO745 END DO746 746 747 ! Compute the daily mean at the end of day748 749 zdaystp = 1.0 / REAL( kdaystp )750 751 IF ( idayend == 0 ) THEN752 DO jj = 1, jpj753 DO ji = 1, jpi754 ! Test if "no night" point755 IF ( icount_sstnight(ji,jj) .NE. 0 ) THEN756 sstdatqc%vdmean(ji,jj) = sstdatqc%vdmean(ji,jj) &747 ! Compute the daily mean at the end of day 748 749 zdaystp = 1.0 / REAL( kdaystp ) 750 751 IF ( idayend == 0 ) THEN 752 DO jj = 1, jpj 753 DO ji = 1, jpi 754 ! Test if "no night" point 755 IF ( icount_sstnight(ji,jj) .NE. 0 ) THEN 756 sstdatqc%vdmean(ji,jj) = sstdatqc%vdmean(ji,jj) & 757 757 & / icount_sstnight(ji,jj) 758 ELSE 759 sstdatqc%vdmean(ji,jj) = zmeanday(ji,jj) * zdaystp 760 ENDIF 758 ELSE 759 sstdatqc%vdmean(ji,jj) = zmeanday(ji,jj) * zdaystp 760 ENDIF 761 END DO 761 762 END DO 762 END DO 763 ENDIF 764 765 ENDIF 763 ENDIF 764 765 ENDIF !ld_nightav 766 766 767 767 ! Get the data for interpolation
Note: See TracChangeset
for help on using the changeset viewer.
