Changeset 15537 for NEMO/branches/UKMO/v15531_cpl-mslp-apr/src/OCE/SBC/sbccpl.F90

Timestamp:

2021-11-25T12:47:25+01:00 (3 years ago)

Author:

jcastill

Message:

Changes compatible with the ticket (not tested in this revision)

File:

• NEMO/branches/UKMO/v15531_cpl-mslp-apr/src/OCE/SBC/sbccpl.F90 (modified) (4 diffs)

NEMO/branches/UKMO/v15531_cpl-mslp-apr/src/OCE/SBC/sbccpl.F90

@@ -220 +220 @@
 #endif
 
-   REAL(wp) ::   rpref = 101000._wp   ! reference atmospheric pressure[N/m2]
-   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rho0)
-
    INTEGER , ALLOCATABLE, SAVE, DIMENSION(:) ::   nrcvinfo           ! OASIS info argument
 
@@ -601 +598 @@
       !                                                      ! Mean Sea Level Pressure   !
       !                                                      ! ------------------------- !
-      srcv(jpr_mslp)%clname = 'O_MSLP'     ;   IF( TRIM(sn_rcv_mslp%cldes  ) == 'coupled' )    srcv(jpr_mslp)%laction = .TRUE.
+      srcv(jpr_mslp)%clname = 'O_MSLP'
+      IF( TRIM(sn_rcv_mslp%cldes  ) == 'coupled' ) THEN
+         srcv(jpr_mslp)%laction = .TRUE.
+         cpl_mslp = .TRUE.
+      ENDIF
       !
       !                                                      ! --------------------------------- !
@@ -1175 +1176 @@
       !!----------------------------------------------------------------------
       USE zdf_oce,  ONLY :   ln_zdfswm
+      USE sbcssm ,  ONLY :   sbc_ssm_cpl   
+      USE lib_fortran     ! distributed memory computing library
       !
       INTEGER, INTENT(in) ::   kt          ! ocean model time step index
@@ -1348 +1351 @@
       !                                                      ! ========================= !
       IF( srcv(jpr_mslp)%laction ) THEN                    ! UKMO SHELF effect of atmospheric pressure on SSH
-          IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields
-
-          r1_grau = 1.e0 / (grav * rho0)               !* constant for optimization
-          ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_grau    ! equivalent ssh (inverse barometer)
-          apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1)                         !atmospheric pressure
-
-          IF( kt == nit000 ) ssh_ibb(:,:) = ssh_ib(:,:)  ! correct this later (read from restart if possible)
+         IF( ln_apr_dyn ) THEN  
+            IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)     !* Swap of ssh_ib fields
+               !                                  !* update the reference atmospheric pressure (if necessary)   
+               IF( ln_ref_apr )  rn_pref = glob_sum( 'sbccpl', frcv(jpr_mslp)%z3(:,:,1) * e1e2t(:,:) ) / tarea
+       
+               ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rn_pref ) * r1_grau  ! equivalent ssh (inverse barometer)   
+               apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1)                        !atmospheric pressure   
+               !  
+               CALL iom_put( "ssh_ib", ssh_ib )   !* output the inverse barometer ssh  
+               !                                         ! ---------------------------------------- !   
+               IF( kt == nit000 ) THEN                   !   set the forcing field at nit000 - 1    !   
+                  !                                      ! ---------------------------------------- !   
+                  !* Restart: read in restart file   
+                  IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN   
+                     IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb read in the restart file'   
+                     CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb )   ! before inv. barometer ssh   
+                  ELSE                                         !* no restart: set from nit000 values   
+                     IF(lwp) WRITE(numout,*) 'sbc_cpl:   ssh_ibb set to nit000 values'   
+                     ssh_ibb(:,:) = ssh_ib(:,:)   
+                  ENDIF   
+               ENDIF   
+               !                                         ! ---------------------------------------- !   
+               IF( lrst_oce ) THEN                       !      Write in the ocean restart file     !   
+                  !                                      ! ---------------------------------------- !   
+                  IF(lwp) WRITE(numout,*)   
+                  IF(lwp) WRITE(numout,*) 'sbc_cpl : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp   
+                  IF(lwp) WRITE(numout,*) '~~~~'   
+                  CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib )   
+               ENDIF   
+            ENDIF  
+        
+            ! Update mean ssh   
+            IF( nn_components /= jp_iam_sas ) CALL sbc_ssm_cpl( kt )
+         ENDIF
       ENDIF
       !