Changeset 11938


Ignore:
Timestamp:
2019-11-20T15:45:47+01:00 (9 months ago)
Author:
dcarneir
Message:

Reverting commit to version where only OBS scripts were updated.

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

Legend:

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

    r9486 r11938  
    2121   !!   dyn_asm_inc    : Apply the dynamic (u and v) increments 
    2222   !!   ssh_asm_inc    : Apply the SSH increment 
    23    !!   seaice_asm_inc : Apply the seaice increment 
     23   !!   seaice_asm_inc : Apply the sea ice increment 
    2424   !!---------------------------------------------------------------------- 
    2525   USE wrk_nemo         ! Memory Allocation 
     
    6262   LOGICAL, PUBLIC :: ln_dyninc = .FALSE.      !: No dynamics (u and v) assimilation increments 
    6363   LOGICAL, PUBLIC :: ln_sshinc = .FALSE.      !: No sea surface height assimilation increment 
    64    LOGICAL, PUBLIC :: ln_seaiceinc             !: No sea ice concentration increment 
     64   LOGICAL, PUBLIC :: ln_seaiceinc = .FALSE.   !: No sea ice concentration increment 
    6565   LOGICAL, PUBLIC :: ln_salfix = .FALSE.      !: Apply minimum salinity check 
    6666   LOGICAL, PUBLIC :: ln_temnofreeze = .FALSE. !: Don't allow the temperature to drop below freezing 
     
    160160         WRITE(numout,*) '      Logical switch for applying SSH increments               ln_sshinc = ', ln_sshinc 
    161161         WRITE(numout,*) '      Logical switch for Direct Initialization (DI)            ln_asmdin = ', ln_asmdin 
    162          WRITE(numout,*) '      Logical switch for applying sea ice increments        ln_seaiceinc = ', ln_seaiceinc 
     162         WRITE(numout,*) '      Logical switch for applying sea ice increments           ln_seaiceinc = ', ln_seaiceinc 
    163163         WRITE(numout,*) '      Logical switch for Incremental Analysis Updating (IAU)   ln_asmiau = ', ln_asmiau 
    164164         WRITE(numout,*) '      Timestep of background in [0,nitend-nit000-1]            nitbkg    = ', nitbkg 
     
    216216 
    217217      IF (      ( ( .NOT. ln_asmdin ).AND.( .NOT. ln_asmiau ) ) & 
    218            .AND.( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ) .OR. ( ln_seaiceinc) )) & 
    219          & CALL ctl_stop( ' One or more of ln_trainc, ln_dyninc, ln_sshinc and ln_seaiceinc is set to .true.', & 
     218           .AND.( ( ln_trainc ).OR.( ln_dyninc ).OR.( ln_sshinc ).OR.( ln_seaiceinc) )) & 
     219         & CALL ctl_stop( ' One or more of ln_trainc, ln_dyninc, ln_sshinc,', & 
     220         &                ' and ln_seaiceinc is set to .true.', & 
    220221         &                ' but ln_asmdin and ln_asmiau are both set to .false. :', & 
    221222         &                ' Inconsistent options') 
     
    227228      IF ( ( .NOT. ln_trainc ).AND.( .NOT. ln_dyninc ).AND.( .NOT. ln_sshinc ).AND.( .NOT. ln_seaiceinc ) & 
    228229         &                     )  & 
    229          & CALL ctl_warn( ' ln_trainc, ln_dyninc, ln_sshinc and ln_seaiceinc are set to .false. :', & 
     230         & CALL ctl_warn( ' ln_trainc, ln_dyninc, ln_sshinc', & 
     231         &                ' and ln_seaiceinc are set to .false. :', & 
    230232         &                ' The assimilation increments are not applied') 
    231233 
     
    922924      ! 
    923925   END SUBROUTINE ssh_asm_inc 
    924  
    925926 
    926927   SUBROUTINE seaice_asm_inc( kt, kindic ) 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90

    r11932 r11938  
    101101   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   fr_iu              !: ice fraction at NEMO U point 
    102102   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   fr_iv              !: ice fraction at NEMO V point 
    103    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_iu              !: ice thickness at NEMO U point 
    104    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_iv              !: ice thickness at NEMO V point 
    105    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_su              !: snow depth at NEMO U point 
    106    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_sv              !: snow depth at NEMO V point 
    107103    
    108104   ! variables used in the coupled interface 
     
    160156                wndj_ice(jpi,jpj)     , nfrzmlt(jpi,jpj)      , ss_iou(jpi,jpj)       , & 
    161157                ss_iov(jpi,jpj)       , fr_iu(jpi,jpj)        , fr_iv(jpi,jpj)        , & 
    162                 thick_iu(jpi,jpj)     , thick_iv(jpi,jpj)     ,                         & 
    163                 thick_su(jpi,jpj)     , thick_sv(jpi,jpj)     ,                         & 
    164                 ht_i(jpi,jpj,ncat)    , ht_s(jpi,jpj,ncat)    ,                         & 
    165158                a_i(jpi,jpj,ncat)     , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , & 
    166159                STAT= ierr(1) ) 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r11935 r11938  
    121121   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   sprecip           !: solid precipitation                          [Kg/m2/s] 
    122122   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   fr_i              !: ice fraction = 1 - lead fraction      (between 0 to 1) 
    123    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_i           !: ice thickness [m] 
    124    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_s           !: snow depth [m] 
    125    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vol_i             !: ice volume [m3] 
    126123#if defined key_cpl_carbon_cycle 
    127124   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   atm_co2           !: atmospheric pCO2                             [ppm] 
     
    170167         ! 
    171168      ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     & 
    172          &      thick_i(jpi,jpj) , thick_s(jpi,jpj) , vol_i(jpi,jpj),     & 
    173169#if defined key_cpl_carbon_cycle 
    174170         &      atm_co2(jpi,jpj) ,                                        & 
  • branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r11935 r11938  
    229229      CALL lbc_lnk ( fr_iu , 'U', 1. ) 
    230230      CALL lbc_lnk ( fr_iv , 'V', 1. ) 
    231  
    232 ! Snow and ice thickness 
    233 ! vice, vsno is volume per unit area of grid cell = thickness 
    234       CALL cice2nemo(vice,thick_i,'T', 1. ) 
    235       CALL cice2nemo(vsno,thick_s,'T', 1. ) 
    236        
    237       IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 
    238          DO jl = 1,ncat 
    239             CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1.) 
    240             CALL cice2nemo(vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1.) 
    241          ENDDO 
    242       ENDIF 
    243        
    244 ! T point to U point 
    245 ! T point to V point 
    246  
    247 ! Sea ice thickness 
    248       thick_iu(:,:)=0.0 
    249       thick_iv(:,:)=0.0 
    250       DO jj=1,jpjm1 
    251          DO ji=1,jpim1 
    252             thick_iu(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji+1,jj))*umask(ji,jj,1) 
    253             thick_iv(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji,jj+1))*vmask(ji,jj,1) 
    254          ENDDO 
    255       ENDDO 
    256  
    257       CALL lbc_lnk ( thick_iu , 'U', 1. ) 
    258       CALL lbc_lnk ( thick_iv , 'V', 1. ) 
    259       
    260 ! Snow depth 
    261       thick_su(:,:)=0.0 
    262       thick_sv(:,:)=0.0 
    263       DO jj=1,jpjm1 
    264          DO ji=1,jpim1 
    265             thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1) 
    266             thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1) 
    267          ENDDO 
    268       ENDDO 
    269          
    270       CALL lbc_lnk ( thick_su , 'U', 1. ) 
    271       CALL lbc_lnk ( thick_sv , 'V', 1. )      
    272231 
    273232      !                                      ! embedded sea ice 
Note: See TracChangeset for help on using the changeset viewer.