Changeset 3585


Ignore:
Timestamp:
2012-11-16T18:14:12+01:00 (8 years ago)
Author:
cbricaud
Message:

commit modifications due to review for SST night operator

Location:
branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/AMM12/EXP00/namelist

    r3436 r3585  
    885885&namobs       !  observation usage switch                               ('key_diaobs') 
    886886!----------------------------------------------------------------------- 
    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 
    918921   profbfiles = 'profiles_01.nc' 
    919                            !     ln_profb_enatim         Enact feedback input time setting switch     
    920                            !     slafilesact             Active SLA input observation file name 
    921                            !     slafilespas             Passive SLA input observation file name  
    922    !                       ! slafbfiles: Feedback SLA input observation file name  
     922                             !     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 
    923926   slafbfiles = 'sla_01.nc' 
    924                            !     sstfiles                GHRSST input observation file name        
    925    !                       ! sstfbfiles: Feedback SST input observation file name  
     927                             !     sstfiles                GHRSST input observation file name 
     928   !                         ! sstfbfiles: Feedback SST input observation file name 
    926929   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    927                            !     seaicefiles             Sea Ice input observation file name  
    928                            !     velavcurfiles           Vel. cur. daily av. input file name   
    929                            !     velhvcurfiles           Vel. cur. high freq. input file name   
    930                            !     velavadcpfiles          Vel. ADCP daily av. input file name     
    931                            !     velhvadcpfiles          Vel. ADCP high freq. input file name  
    932                            !     velfbfiles              Vel. feedback input observation file name  
    933                            !     dobsini                 Initial date in window YYYYMMDD.HHMMSS        
    934                            !     dobsend                 Final date in window YYYYMMDD.HHMMSS          
    935                            !     n1dint                  Type of vertical interpolation method         
    936                            !     n2dint                  Type of horizontal interpolation method        
    937                            !     ln_nea                  Rejection of observations near land switch     
    938    nmsshc     = 0          ! MSSH correction scheme                          
    939                            !     mdtcorr                 MDT  correction                                
    940                            !     mdtcutoff               MDT cutoff for computed correction           
    941    ln_altbias = .false.    ! Logical switch for alt bias                 
    942    ln_ignmis  = .true.     ! Logical switch for ignoring missing files    
    943                            !     endailyavtypes   ENACT daily average types                     
     930                             !     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 
    944947   ln_grid_global = .true. 
    945948   ln_grid_search_lookup = .false. 
    946 /  
     949/ 
    947950!----------------------------------------------------------------------- 
    948951&nam_asminc   !   assimilation increments                               ('key_asminc') 
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/GYRE/EXP00/namelist

    r3436 r3585  
    877877&namobs       !  observation usage switch                               ('key_diaobs') 
    878878!----------------------------------------------------------------------- 
    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  
    910913   profbfiles = 'profiles_01.nc' 
    911                            !     ln_profb_enatim         Enact feedback input time setting switch     
    912                            !     slafilesact             Active SLA input observation file name 
    913                            !     slafilespas             Passive SLA input observation file name  
    914    !                       ! slafbfiles: Feedback SLA input observation file name  
     914                             !     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  
    915918   slafbfiles = 'sla_01.nc' 
    916                            !     sstfiles                GHRSST input observation file name        
    917    !                       ! sstfbfiles: Feedback SST input observation file name  
     919                             !     sstfiles                GHRSST input observation file name        
     920   !                         ! sstfbfiles: Feedback SST input observation file name  
    918921   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    919                            !     seaicefiles             Sea Ice input observation file name  
    920                            !     velavcurfiles           Vel. cur. daily av. input file name   
    921                            !     velhvcurfiles           Vel. cur. high freq. input file name   
    922                            !     velavadcpfiles          Vel. ADCP daily av. input file name     
    923                            !     velhvadcpfiles          Vel. ADCP high freq. input file name  
    924                            !     velfbfiles              Vel. feedback input observation file name  
    925                            !     dobsini                 Initial date in window YYYYMMDD.HHMMSS        
    926                            !     dobsend                 Final date in window YYYYMMDD.HHMMSS          
    927                            !     n1dint                  Type of vertical interpolation method         
    928                            !     n2dint                  Type of horizontal interpolation method        
    929                            !     ln_nea                  Rejection of observations near land switch     
    930    nmsshc     = 0          ! MSSH correction scheme                          
    931                            !     mdtcorr                 MDT  correction                                
    932                            !     mdtcutoff               MDT cutoff for computed correction           
    933    ln_altbias = .false.    ! Logical switch for alt bias                 
    934    ln_ignmis  = .true.     ! Logical switch for ignoring missing files    
    935                            !     endailyavtypes   ENACT daily average types                     
     922                             !     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 
    936939   ln_grid_global = .true. 
    937940   ln_grid_search_lookup = .false. 
    938 /  
     941/ 
    939942!----------------------------------------------------------------------- 
    940943&nam_asminc   !   assimilation increments                               ('key_asminc') 
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist

    r3436 r3585  
    2727   cn_exp      =  "ORCA2"  !  experience name  
    2828   nn_it000    =       1   !  first time step 
    29    nn_itend    =    5475   !  last  time step (std 5475) 
     29   nn_itend    =      15   !  last  time step (std 5475) 
    3030   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    3131   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
     
    769769   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    770770   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) 
    774774/ 
    775775!----------------------------------------------------------------------- 
     
    881881&namobs       !  observation usage switch                               ('key_diaobs') 
    882882!----------------------------------------------------------------------- 
    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  
    914917   profbfiles = 'profiles_01.nc' 
    915                            !     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  
     918                             !     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  
    919922   slafbfiles = 'sla_01.nc' 
    920                            !     sstfiles                GHRSST input observation file name        
    921    !                       ! sstfbfiles: Feedback SST input observation file name  
     923                             !     sstfiles                GHRSST input observation file name        
     924   !                         ! sstfbfiles: Feedback SST input observation file name  
    922925   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    923                            !     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                     
     926                             !     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                     
    940943   ln_grid_global = .true. 
    941944   ln_grid_search_lookup = .false. 
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r3436 r3585  
    878878&namobs       !  observation usage switch                               ('key_diaobs') 
    879879!----------------------------------------------------------------------- 
    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  
    911914   profbfiles = 'profiles_01.nc' 
    912                            !     ln_profb_enatim         Enact feedback input time setting switch     
    913                            !     slafilesact             Active SLA input observation file name 
    914                            !     slafilespas             Passive SLA input observation file name  
    915    !                       ! slafbfiles: Feedback SLA input observation file name  
     915                             !     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  
    916919   slafbfiles = 'sla_01.nc' 
    917                            !     sstfiles                GHRSST input observation file name        
    918    !                       ! sstfbfiles: Feedback SST input observation file name  
     920                             !     sstfiles                GHRSST input observation file name        
     921   !                         ! sstfbfiles: Feedback SST input observation file name  
    919922   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    920                            !     seaicefiles             Sea Ice input observation file name  
    921                            !     velavcurfiles           Vel. cur. daily av. input file name   
    922                            !     velhvcurfiles           Vel. cur. high freq. input file name   
    923                            !     velavadcpfiles          Vel. ADCP daily av. input file name     
    924                            !     velhvadcpfiles          Vel. ADCP high freq. input file name  
    925                            !     velfbfiles              Vel. feedback input observation file name  
    926                            !     dobsini                 Initial date in window YYYYMMDD.HHMMSS        
    927                            !     dobsend                 Final date in window YYYYMMDD.HHMMSS          
    928                            !     n1dint                  Type of vertical interpolation method         
    929                            !     n2dint                  Type of horizontal interpolation method        
    930                            !     ln_nea                  Rejection of observations near land switch     
    931    nmsshc     = 0          ! MSSH correction scheme                          
    932                            !     mdtcorr                 MDT  correction                                
    933                            !     mdtcutoff               MDT cutoff for computed correction           
    934    ln_altbias = .false.    ! Logical switch for alt bias                 
    935    ln_ignmis  = .true.     ! Logical switch for ignoring missing files    
    936                            !     endailyavtypes   ENACT daily average types                     
     923                             !     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                     
    937940   ln_grid_global = .true. 
    938941   ln_grid_search_lookup = .false. 
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r3411 r3585  
    5959 
    6060   !! * 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 
    8689 
    8790   REAL(KIND=dp), PUBLIC :: dobsini   !: Observation window start date YYYYMMDD.HHMMSS 
     
    107110      & ld_velav     !: Velocity data is daily averaged 
    108111   LOGICAL, DIMENSION(:), ALLOCATABLE :: & 
    109       & ld_sstnight  !: SST observation corresponds to night mean 
     112      & l_sstnight  !: SST observation corresponds to night mean 
    110113 
    111114   !!---------------------------------------------------------------------- 
     
    159162      NAMELIST/namobs/ln_ena, ln_cor, ln_profb, ln_t3d, ln_s3d,       & 
    160163         &            ln_sla, ln_sladt, ln_slafb,                     & 
    161          &            ln_ssh, ln_sst, ln_sstfb, ln_sss, ln_nea,       & 
     164         &            ln_ssh, ln_sst, ln_sss, ln_nea,                 & 
    162165         &            enactfiles, coriofiles, profbfiles,             & 
    163166         &            slafilesact, slafilespas, slafbfiles,           & 
     
    166169         &            dobsini, dobsend, n1dint, n2dint,               & 
    167170         &            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,                          & 
    169175         &            ln_grid_search_lookup,                          & 
    170176         &            grid_search_file, grid_search_res,              & 
     
    212218      ln_sst = .FALSE. 
    213219      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. 
    216226      ln_sss = .FALSE. 
    217227      ln_profb = .FALSE. 
     
    220230      ln_sladt = .TRUE. 
    221231      ln_slafb = .FALSE. 
    222       ln_sstfb = .FALSE. 
    223232      ln_velavcur = .FALSE. 
    224233      ln_velhrcur = .FALSE. 
     
    355364         WRITE(numout,*) '             Logical switch for SST observations                ln_sst = ', ln_sst 
    356365         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     
    357367         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     
    358369         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     
    359371         WRITE(numout,*) '             Logical switch for SSS observations                ln_sss = ', ln_sss 
    360372         WRITE(numout,*) '             Logical switch for Sea Ice observations         ln_seaice = ', ln_seaice 
     
    721733 
    722734      !  - Sea surface temperature 
    723       IF ( ln_sst ) THEN 
     735      IF( ln_sst ) THEN 
    724736 
    725737         ! Set the number of variables for sst to 1 
     
    739751         ALLOCATE(sstdata(nsstsets)) 
    740752         ALLOCATE(sstdatqc(nsstsets)) 
    741          ALLOCATE(ld_sstnight(nsstsets)) 
    742          sstdata(:)%nsurf=0 
    743          sstdatqc(:)%nsurf=0     
    744          ld_sstnight(:)=.false. 
     753         ALLOCATE(l_sstnight(nsstsets)) 
     754         sstdata(:)%nsurf  = 0 
     755         sstdatqc(:)%nsurf = 0     
     756         l_sstnight(:)    = .false. 
    745757 
    746758         nsstsets = 0 
     
    750762            nsstsets = nsstsets + 1 
    751763 
    752             ld_sstnight(nsstsets) = .TRUE. 
     764            l_sstnight(nsstsets) = ln_reysst_night 
    753765 
    754766            CALL obs_rea_sst_rey( reysstname, reysstfmt, sstdata(nsstsets), & 
     
    758770               &              ln_nea ) 
    759771 
    760         ENDIF 
     772         ENDIF 
    761773         
    762         IF (ln_ghrsst) THEN 
     774         IF (ln_ghrsst) THEN 
    763775         
    764776            nsstsets = nsstsets + 1 
    765777 
    766             ld_sstnight(nsstsets) = .TRUE. 
     778            l_sstnight(nsstsets) = ln_ghrsst_night 
    767779           
    768780            CALL obs_rea_sst( 1, sstdata(nsstsets), jnumsst, & 
     
    783795               nsstsets = nsstsets + 1 
    784796 
    785                ld_sstnight(nsstsets) = .TRUE. 
     797               l_sstnight(nsstsets) = ln_sstfb_night 
    786798             
    787799               CALL obs_rea_sst( 0, sstdata(nsstsets), 1, & 
     
    11051117               &              kstp, jpi, jpj, nit000, idaystp,  & 
    11061118               &              tsn(:,:,1,jp_tem), tmask(:,:,1),  & 
    1107                &              n2dint, ld_sstnight(jsstset) ) 
     1119               &              n2dint, l_sstnight(jsstset) ) 
    11081120         END DO 
    11091121      ENDIF 
  • branches/2012/dev_r3342_MERCATOR7_SST/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r3411 r3585  
    663663         & psstn,  &    ! Model SST field 
    664664         & psstmask     ! Land-sea mask 
     665      LOGICAL, INTENT(IN) :: ld_nightav 
    665666 
    666667      !! * Local declarations 
     
    694695         & igrdi, & 
    695696         & igrdj 
    696       LOGICAL, INTENT(IN) :: ld_nightav 
    697697 
    698698      !----------------------------------------------------------------------- 
     
    705705      IF ( ld_nightav ) THEN 
    706706 
    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 
    724737         DO jj = 1, jpj 
    725738            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) 
    729744            END DO 
    730745         END DO 
    731       ENDIF 
    732  
    733       zintmp(:,:) = 0.0 
    734       zouttmp(:,:) = sbc_dcy( zintmp(:,:), .TRUE. ) 
    735       imask_night(:,:) = INT( zouttmp(:,:) ) 
    736  
    737       DO jj = 1, jpj 
    738          DO ji = 1, jpi 
    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) 
    744          END DO 
    745       END DO 
    746746    
    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) & 
     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) & 
    757757                    &                        / 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 
    761762            END DO 
    762          END DO 
    763       ENDIF 
    764  
    765       ENDIF 
     763         ENDIF 
     764 
     765      ENDIF !ld_nightav 
    766766 
    767767      ! Get the data for interpolation 
Note: See TracChangeset for help on using the changeset viewer.