Changeset 12434


Ignore:
Timestamp:
2020-02-21T14:51:51+01:00 (12 months ago)
Author:
jcastill
Message:

Merge development branch AMM15_v3_6_STABLE_package_collate_utils325, see Met Office utils branch 325. As per reviewer's comments, indentation has been included for the added if sentence

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_collate/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r11277 r12434  
    12111211      !    
    12121212      IF( srcv(jpr_mslp)%laction ) THEN                    ! UKMO SHELF effect of atmospheric pressure on SSH    
    1213           IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields    
    1214         
    1215           !                                                  !* update the reference atmospheric pressure (if necessary)  
    1216           IF( ln_ref_apr )  rn_pref = glob_sum( frcv(jpr_mslp)%z3(:,:,1) * e1e2t(:,:) ) / tarea  
    1217          
    1218           ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rn_pref ) * r1_grau    ! equivalent ssh (inverse barometer)    
    1219           apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1) !atmospheric pressure    
    1220           !  
    1221           CALL iom_put( "ssh_ib", ssh_ib )                                    !* output the inverse barometer ssh 
    1222         
    1223           !                                         ! ---------------------------------------- !  
    1224           IF( kt == nit000 ) THEN                   !   set the forcing field at nit000 - 1    !  
    1225              !                                      ! ---------------------------------------- !  
    1226              !* Restart: read in restart file  
    1227              IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN  
    1228                 IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb read in the restart file'  
    1229                 CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb )   ! before inv. barometer ssh  
    1230              ELSE                                         !* no restart: set from nit000 values  
    1231                 IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb set to nit000 values'  
    1232                 ssh_ibb(:,:) = ssh_ib(:,:)  
     1213          IF( ln_apr_dyn ) THEN 
     1214             IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields    
     1215            
     1216             !                                                  !* update the reference atmospheric pressure (if necessary)  
     1217             IF( ln_ref_apr )  rn_pref = glob_sum( frcv(jpr_mslp)%z3(:,:,1) * e1e2t(:,:) ) / tarea  
     1218            
     1219             ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rn_pref ) * r1_grau    ! equivalent ssh (inverse barometer)    
     1220             apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1) !atmospheric pressure    
     1221             !  
     1222             CALL iom_put( "ssh_ib", ssh_ib )                                    !* output the inverse barometer ssh 
     1223            
     1224             !                                         ! ---------------------------------------- !  
     1225             IF( kt == nit000 ) THEN                   !   set the forcing field at nit000 - 1    !  
     1226                !                                      ! ---------------------------------------- !  
     1227                !* Restart: read in restart file  
     1228                IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN  
     1229                   IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb read in the restart file'  
     1230                   CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb )   ! before inv. barometer ssh  
     1231                ELSE                                         !* no restart: set from nit000 values  
     1232                   IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb set to nit000 values'  
     1233                   ssh_ibb(:,:) = ssh_ib(:,:)  
     1234                ENDIF  
    12331235             ENDIF  
    1234           ENDIF  
    1235           !                                         ! ---------------------------------------- !  
    1236           IF( lrst_oce ) THEN                       !      Write in the ocean restart file    !  
    1237              !                                      ! ---------------------------------------- !  
    1238              IF(lwp) WRITE(numout,*)  
    1239              IF(lwp) WRITE(numout,*) 'sbc_cpl : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp  
    1240              IF(lwp) WRITE(numout,*) '~~~~'  
    1241              CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib )  
     1236             !                                         ! ---------------------------------------- !  
     1237             IF( lrst_oce ) THEN                       !      Write in the ocean restart file    !  
     1238                !                                      ! ---------------------------------------- !  
     1239                IF(lwp) WRITE(numout,*)  
     1240                IF(lwp) WRITE(numout,*) 'sbc_cpl : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp  
     1241                IF(lwp) WRITE(numout,*) '~~~~'  
     1242                CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib )  
     1243             ENDIF  
    12421244          ENDIF  
    12431245         
Note: See TracChangeset for help on using the changeset viewer.