Changeset 7162 for branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2016-11-01T14:23:51+01:00 (7 years ago)

Author:

cetlod

Message:

new top interface : Add PISCES quota model

File:

• branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (5 diffs)

branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -20 +20 @@
    USE dom_oce        ! ocean space and time domain
    USE sbc_oce        ! Surface boundary condition: ocean fields
+   USE trc_oce        ! share SMS/Ocean variables
    USE sbc_ice        ! Surface boundary condition: ice fields
    USE sbcapr         ! Stochastic param. : ???
@@ -36 +37 @@
    USE albedo         ! 
    USE eosbn2         ! 
-   USE sbcrnf  , ONLY : l_rnfcpl
-#if defined key_cpl_carbon_cycle
-   USE p4zflx, ONLY : oce_co2
-#endif
+   USE sbcrnf, ONLY : l_rnfcpl
 #if defined key_cice
    USE ice_domain_size, only: ncat
@@ -467 +465 @@
       !                                                      !      Atmospheric CO2      !
       !                                                      ! ------------------------- !
-      srcv(jpr_co2 )%clname = 'O_AtmCO2'   ;   IF( TRIM(sn_rcv_co2%cldes   ) == 'coupled' )    srcv(jpr_co2 )%laction = .TRUE.
+      srcv(jpr_co2 )%clname = 'O_AtmCO2'   
+      IF( TRIM(sn_rcv_co2%cldes   ) == 'coupled' )  THEN
+         srcv(jpr_co2 )%laction = .TRUE.
+         l_co2cpl = .TRUE.
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) '   Atmospheric pco2 received from oasis '
+         IF(lwp) WRITE(numout,*)
+      ENDIF
       !                                                      ! ------------------------- !
       !                                                      !   topmelt and botmelt     !   
@@ -984 +989 @@
       ENDIF
 
-#if defined key_cpl_carbon_cycle
       !                                                      ! ================== !
       !                                                      ! atmosph. CO2 (ppm) !
       !                                                      ! ================== !
       IF( srcv(jpr_co2)%laction )   atm_co2(:,:) = frcv(jpr_co2)%z3(:,:,1)
-#endif
 
       !  Fields received by SAS when OASIS coupling
@@ -1919 +1922 @@
          IF( ssnd(jps_hsnw)%laction )   CALL cpl_snd( jps_hsnw, isec, ztmp4, info )
       ENDIF
-      !
-#if defined key_cpl_carbon_cycle
       !                                                      ! ------------------------- !
       !                                                      !  CO2 flux from PISCES     ! 
       !                                                      ! ------------------------- !
-      IF( ssnd(jps_co2)%laction )   CALL cpl_snd( jps_co2, isec, RESHAPE ( oce_co2, (/jpi,jpj,1/) ) , info )
-      !
-#endif
+      IF( ssnd(jps_co2)%laction .AND. l_co2cpl )   CALL cpl_snd( jps_co2, isec, RESHAPE ( oce_co2, (/jpi,jpj,1/) ) , info )
+      !
       !                                                      ! ------------------------- !
       IF( ssnd(jps_ocx1)%laction ) THEN                      !      Surface current      !