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

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

Changeset 12609 for branches


Ignore:
Timestamp:
2020-03-26T11:53:23+01:00 (4 years ago)
Author:
dcarneir
Message:

Making SIT changes in global configs to be compatible with build in shelf configs

Location:
branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r12594 r12609  
    6060#endif 
    6161   LOGICAL, PUBLIC :: ln_bkgwri = .FALSE.      !: No output of the background state fields 
     62   LOGICAL, PUBLIC :: ln_avgbkg = .FALSE.      !: No output of the mean background state fields 
    6263   LOGICAL, PUBLIC :: ln_asmiau = .FALSE.      !: No applying forcing with an assimilation increment 
    6364   LOGICAL, PUBLIC :: ln_asmdin = .FALSE.      !: No direct initialization 
     
    133134      REAL(wp), POINTER, DIMENSION(:,:) ::   hdiv   ! 2D workspace 
    134135      !! 
    135       NAMELIST/nam_asminc/ ln_bkgwri, ln_balwri,                           & 
     136      NAMELIST/nam_asminc/ ln_bkgwri, ln_avgbkg, ln_balwri,                & 
    136137         &                 ln_trainc, ln_dyninc, ln_sshinc,                & 
    137          &                 ln_phytobal, ln_slchltotinc, ln_slchldiainc,    & 
    138          &                 ln_slchlnoninc, ln_schltotinc, ln_slphytotinc,  & 
    139          &                 ln_slphydiainc, ln_slphynoninc, ln_spco2inc,    & 
    140          &                 ln_sfco2inc, ln_plchltotinc, ln_pchltotinc,     & 
    141          &                 ln_pno3inc, ln_psi4inc, ln_pdicinc, ln_palkinc, & 
    142          &                 ln_pphinc, ln_po2inc,                           & 
    143138         &                 ln_asmdin, ln_asmiau,                           & 
    144139         &                 nitbkg, nitdin, nitiaustr, nitiaufin, niaufn,   & 
    145          &                 ln_salfix, salfixmin, nn_divdmp,                & 
    146          &                 ln_seaiceinc, ln_sitinc, ln_temnofreeze,        & 
    147          &                 mld_choice_bgc, rn_maxchlinc 
     140         &                 ln_salfix, salfixmin, nn_divdmp 
    148141      !!---------------------------------------------------------------------- 
    149142 
     
    171164         WRITE(numout,*) '   Namelist nam_asminc : set assimilation increment parameters' 
    172165         WRITE(numout,*) '      Logical switch for writing out background state          ln_bkgwri = ', ln_bkgwri 
     166         WRITE(numout,*) '      Logical switch for writing mean background state         ln_avgbkg = ', ln_avgbkg 
    173167         WRITE(numout,*) '      Logical switch for writing out balancing increments      ln_balwri = ', ln_balwri 
    174168         WRITE(numout,*) '      Logical switch for applying tracer increments            ln_trainc = ', ln_trainc 
    175169         WRITE(numout,*) '      Logical switch for applying velocity increments          ln_dyninc = ', ln_dyninc 
    176170         WRITE(numout,*) '      Logical switch for applying SSH increments               ln_sshinc = ', ln_sshinc 
    177          WRITE(numout,*) '      Logical switch for Direct Initialization (DI)            ln_asmdin = ', ln_asmdin 
    178171         WRITE(numout,*) '      Logical switch for applying SIC increments               ln_seaiceinc = ', ln_seaiceinc 
    179172         WRITE(numout,*) '      Logical switch for applying SIT increments               ln_sitinc = ', ln_sitinc 
    180          WRITE(numout,*) '      Logical switch for phytoplankton balancing             ln_phytobal = ', ln_phytobal 
    181          WRITE(numout,*) '      Logical switch for applying slchltot increments     ln_slchltotinc = ', ln_slchltotinc 
    182          WRITE(numout,*) '      Logical switch for applying slchldia increments     ln_slchldiainc = ', ln_slchldiainc 
    183          WRITE(numout,*) '      Logical switch for applying slchlnon increments     ln_slchlnoninc = ', ln_slchlnoninc 
    184          WRITE(numout,*) '      Logical switch for applying schltot increments       ln_schltotinc = ', ln_schltotinc 
    185          WRITE(numout,*) '      Logical switch for applying slphytot increments     ln_slphytotinc = ', ln_slphytotinc 
    186          WRITE(numout,*) '      Logical switch for applying slphydia increments     ln_slphydiainc = ', ln_slphydiainc 
    187          WRITE(numout,*) '      Logical switch for applying slphynon increments     ln_slphynoninc = ', ln_slphynoninc 
    188          WRITE(numout,*) '      Logical switch for applying spco2 increments           ln_spco2inc = ', ln_spco2inc 
    189          WRITE(numout,*) '      Logical switch for applying sfco2 increments           ln_sfco2inc = ', ln_sfco2inc 
    190          WRITE(numout,*) '      Logical switch for applying plchltot increments     ln_plchltotinc = ', ln_plchltotinc 
    191          WRITE(numout,*) '      Logical switch for applying pchltot increments       ln_pchltotinc = ', ln_pchltotinc 
    192          WRITE(numout,*) '      Logical switch for applying pno3 increments             ln_pno3inc = ', ln_pno3inc 
    193          WRITE(numout,*) '      Logical switch for applying psi4 increments             ln_psi4inc = ', ln_psi4inc 
    194          WRITE(numout,*) '      Logical switch for applying pdic increments             ln_pdicinc = ', ln_pdicinc 
    195          WRITE(numout,*) '      Logical switch for applying palk increments             ln_palkinc = ', ln_palkinc 
    196          WRITE(numout,*) '      Logical switch for applying pph increments               ln_pphinc = ', ln_pphinc 
    197          WRITE(numout,*) '      Logical switch for applying po2 increments               ln_po2inc = ', ln_po2inc 
     173         WRITE(numout,*) '      Logical switch for Direct Initialization (DI)            ln_asmdin = ', ln_asmdin 
    198174         WRITE(numout,*) '      Logical switch for Incremental Analysis Updating (IAU)   ln_asmiau = ', ln_asmiau 
    199175         WRITE(numout,*) '      Timestep of background in [0,nitend-nit000-1]            nitbkg    = ', nitbkg 
     
    204180         WRITE(numout,*) '      Logical switch for ensuring that the sa > salfixmin      ln_salfix = ', ln_salfix 
    205181         WRITE(numout,*) '      Minimum salinity after applying the increments           salfixmin = ', salfixmin 
    206          WRITE(numout,*) '      Choice of MLD for BGC assimilation                  mld_choice_bgc = ', mld_choice_bgc 
    207          WRITE(numout,*) '      Maximum absolute chlorophyll increment (<=0 = off)    rn_maxchlinc = ', rn_maxchlinc 
    208182      ENDIF 
    209183 
     
    365339      !-------------------------------------------------------------------- 
    366340 
    367       ALLOCATE( t_bkginc(jpi,jpj,jpk) ) 
    368       ALLOCATE( s_bkginc(jpi,jpj,jpk) ) 
    369       ALLOCATE( u_bkginc(jpi,jpj,jpk) ) 
    370       ALLOCATE( v_bkginc(jpi,jpj,jpk) ) 
    371       ALLOCATE( ssh_bkginc(jpi,jpj)   ) 
    372       ALLOCATE( seaice_bkginc(jpi,jpj)) 
    373       ALLOCATE( sit_bkginc(jpi,jpj)   ) 
     341      IF ( ln_trainc ) THEN 
     342         ALLOCATE( t_bkginc(jpi,jpj,jpk) ) 
     343         ALLOCATE( s_bkginc(jpi,jpj,jpk) ) 
     344         t_bkginc(:,:,:) = 0.0 
     345         s_bkginc(:,:,:) = 0.0 
     346      ENDIF 
     347      IF ( ln_dyninc ) THEN  
     348         ALLOCATE( u_bkginc(jpi,jpj,jpk) ) 
     349         ALLOCATE( v_bkginc(jpi,jpj,jpk) ) 
     350         u_bkginc(:,:,:) = 0.0 
     351         v_bkginc(:,:,:) = 0.0 
     352      ENDIF 
     353      IF ( ln_sshinc ) THEN 
     354         ALLOCATE( ssh_bkginc(jpi,jpj)   ) 
     355         ssh_bkginc(:,:) = 0.0 
     356      ENDIF 
     357      IF ( ln_seaiceinc ) THEN  
     358         ALLOCATE( seaice_bkginc(jpi,jpj)) 
     359         seaice_bkginc(:,:) = 0.0 
     360      ENDIF 
     361      IF ( ln_sitinc ) THEN  
     362         ALLOCATE( sit_bkginc(jpi,jpj)) 
     363         sit_bkginc(:,:) = 0.0 
     364      ENDIF 
    374365#if defined key_asminc 
    375366      ALLOCATE( ssh_iau(jpi,jpj)      ) 
    376 #endif 
    377       t_bkginc(:,:,:) = 0.0 
    378       s_bkginc(:,:,:) = 0.0 
    379       u_bkginc(:,:,:) = 0.0 
    380       v_bkginc(:,:,:) = 0.0 
    381       ssh_bkginc(:,:) = 0.0 
    382       seaice_bkginc(:,:) = 0.0 
    383       sit_bkginc(:,:) = 0.0 
    384 #if defined key_asminc 
    385367      ssh_iau(:,:)    = 0.0 
    386368#endif 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r12380 r12609  
    968968            ENDIF 
    969969 
     970#if defined key_cice 
    970971            IF ( TRIM(cobstypessurf(jtype)) == 'sit' ) THEN 
    971972               CALL obs_rea_snowdepth( surfdataqc(jtype), n2dintsurf(jtype), thick_s(:,:) ) 
    972973            ENDIF  
     974#endif 
    973975 
    974976            IF ( TRIM(cobstypessurf(jtype)) == 'sst' .AND. ln_sstbias ) THEN 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r12380 r12609  
    950950         IF ( surfdataqc%lclim ) surfdataqc%rclm(jobs,1) = zclm(1) 
    951951 
     952#if defined key_cice 
    952953         IF ( TRIM(surfdataqc%cvars(1)) == 'FBD' ) THEN 
    953954            ! Convert radar freeboard to true freeboard (add 1/4 snow depth; 1/4 based on ratio of speed of light in vacuum compared to snow (3.0e8 vs 2.4e8 m/s)) 
     
    971972            surfdataqc%nqc(jobs) = 4 
    972973         ENDIF 
     974#endif defined key_cice 
    973975          
    974976      END DO 
Note: See TracChangeset for help on using the changeset viewer.