Changeset 9259 for branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2018-01-18T13:14:23+01:00 (5 years ago)

Author:

jcastill

Message:

Changes to stop reading and using the pressure forcing file when coupling to the atmospheric mean sea level pressure

File:

• branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (4 diffs)

branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -185 +185 @@
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   albedo_oce_mix     ! ocean albedo sent to atmosphere (mix clear/overcast sky)
    
-   REAL(wp) ::   rpref = 101000._wp   ! reference atmospheric pressure[N/m2]  
-   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rau0
-
    INTEGER , ALLOCATABLE, SAVE, DIMENSION(    :) ::   nrcvinfo           ! OASIS info argument
 
@@ -523 +520 @@
       !                                                      ! 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
       
       !                                                      ! ------------------------- !
@@ -977 +978 @@
       !!----------------------------------------------------------------------
       USE sbcflx ,  ONLY : ln_shelf_flx
+      USE sbcssm ,  ONLY : sbc_ssm_cpl
+      USE lib_fortran     ! distributed memory computing library
 
       INTEGER, INTENT(in)           ::   kt          ! ocean model time step index
@@ -1175 +1178 @@
       !  
       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  
+         IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields  
       
-          r1_grau = 1.e0 / (grav * rau0)               !* 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  
+         !                                                  !* update the reference atmospheric pressure (if necessary)
+         IF( ln_ref_apr )  rn_pref = glob_sum( 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 ) ssh_ibb(:,:) = ssh_ib(:,:)  ! correct this later (read from restart if possible)  
+         !                                         ! ---------------------------------------- !
+         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
+
+         ! Update mean ssh
+         CALL sbc_ssm_cpl( kt )
       END IF  
       !