Changeset 3557 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

Timestamp:

2012-11-14T18:04:49+01:00 (11 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB minor modifications for PISCES when using AGRIF

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90 (modified) (3 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

@@ -274 +274 @@
          IF( sn_dust%ln_tint )     ALLOCATE( sf_dust(1)%fdta(jpi,jpj,1,2) )
          !
-         ! Get total input dust ; need to compute total atmospheric supply of Si in a year
-         CALL iom_open (  TRIM( sn_dust%clname ) , numdust )
-         CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust)  ! get number of record in file
-         ALLOCATE( zdust(jpi,jpj,ntimes_dust) )
-         DO jm = 1, ntimes_dust
-            CALL iom_get( numdust, jpdom_data, TRIM( sn_dust%clvar ), zdust(:,:,jm), jm )
-         END DO
-         CALL iom_close( numdust )
-         ztimes_dust = 1._wp / FLOAT( ntimes_dust ) 
-         sumdepsi = 0.e0
-         DO jm = 1, ntimes_dust
-            sumdepsi = sumdepsi + glob_sum( zdust(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_dust )
-         ENDDO
-         sumdepsi = sumdepsi / ( nyear_len(1) * rday ) * 12. * 8.8 * 0.075 / 28.1 
-         DEALLOCATE( zdust)
+         IF( Agrif_Root() ) THEN   !  Only on the master grid
+            ! Get total input dust ; need to compute total atmospheric supply of Si in a year
+            CALL iom_open (  TRIM( sn_dust%clname ) , numdust )
+            CALL iom_gettime( numdust, zsteps, kntime=ntimes_dust)  ! get number of record in file
+            ALLOCATE( zdust(jpi,jpj,ntimes_dust) )
+            DO jm = 1, ntimes_dust
+               CALL iom_get( numdust, jpdom_data, TRIM( sn_dust%clvar ), zdust(:,:,jm), jm )
+            END DO
+            CALL iom_close( numdust )
+            ztimes_dust = 1._wp / FLOAT( ntimes_dust ) 
+            sumdepsi = 0.e0
+            DO jm = 1, ntimes_dust
+               sumdepsi = sumdepsi + glob_sum( zdust(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_dust )
+            ENDDO
+            sumdepsi = sumdepsi / ( nyear_len(1) * rday ) * 12. * 8.8 * 0.075 / 28.1 
+            DEALLOCATE( zdust)
+         ENDIF
       ELSE
          sumdepsi  = 0._wp
@@ -335 +337 @@
             IF( slf_river(ifpr)%ln_tint ) ALLOCATE( sf_river(ifpr)%fdta(jpi,jpj,1,2) )
          END DO
-         ! Get total input rivers ; need to compute total river supply in a year
-         DO ifpr = 1, jpriv
-            CALL iom_open ( TRIM( slf_river(ifpr)%clname ), numriv )
-            CALL iom_gettime( numriv, zsteps, kntime=ntimes_riv)
-            ALLOCATE( zriver(jpi,jpj,ntimes_riv) )
-            DO jm = 1, ntimes_riv
-               CALL iom_get( numriv, jpdom_data, TRIM( slf_river(ifpr)%clvar ), zriver(:,:,jm), jm )
-            END DO
-            CALL iom_close( numriv )
-            ztimes_riv = 1._wp / FLOAT(ntimes_riv) 
-            DO jm = 1, ntimes_riv
-               rivinput(ifpr) = rivinput(ifpr) + glob_sum( zriver(:,:,jm) * tmask(:,:,1) * ztimes_riv ) 
-            END DO
-            DEALLOCATE( zriver)
-         END DO
-         ! N/P and Si releases due to coastal rivers
-         ! -----------------------------------------
-         rivdicinput = (rivinput(jr_dic) + rivinput(jr_doc) ) * 1E3 / 12._wp
-         rivdininput = (rivinput(jr_din) + rivinput(jr_don) ) * 1E3 / rno3 / 14._wp
-         rivdipinput = (rivinput(jr_dip) + rivinput(jr_dop) ) * 1E3 / po4r / 31._wp
-         rivdsiinput = rivinput(jr_dsi) * 1E3 / 28.1_wp
-         rivalkinput = rivinput(jr_dic) * 1E3 / 12._wp
-         !
+         IF( Agrif_Root() ) THEN   !  Only on the master grid
+            ! Get total input rivers ; need to compute total river supply in a year
+            DO ifpr = 1, jpriv
+               CALL iom_open ( TRIM( slf_river(ifpr)%clname ), numriv )
+               CALL iom_gettime( numriv, zsteps, kntime=ntimes_riv)
+               ALLOCATE( zriver(jpi,jpj,ntimes_riv) )
+               DO jm = 1, ntimes_riv
+                  CALL iom_get( numriv, jpdom_data, TRIM( slf_river(ifpr)%clvar ), zriver(:,:,jm), jm )
+               END DO
+               CALL iom_close( numriv )
+               ztimes_riv = 1._wp / FLOAT(ntimes_riv) 
+               DO jm = 1, ntimes_riv
+                  rivinput(ifpr) = rivinput(ifpr) + glob_sum( zriver(:,:,jm) * tmask(:,:,1) * ztimes_riv ) 
+               END DO
+               DEALLOCATE( zriver)
+            END DO
+            ! N/P and Si releases due to coastal rivers
+            ! -----------------------------------------
+            rivdicinput = (rivinput(jr_dic) + rivinput(jr_doc) ) * 1E3 / 12._wp
+            rivdininput = (rivinput(jr_din) + rivinput(jr_don) ) * 1E3 / rno3 / 14._wp
+            rivdipinput = (rivinput(jr_dip) + rivinput(jr_dop) ) * 1E3 / po4r / 31._wp
+            rivdsiinput = rivinput(jr_dsi) * 1E3 / 28.1_wp
+            rivalkinput = rivinput(jr_dic) * 1E3 / 12._wp
+            !
+         ENDIF
       ELSE
          rivdicinput = 0._wp
@@ -382 +386 @@
          IF( sn_ndepo%ln_tint )    ALLOCATE( sf_ndepo(1)%fdta(jpi,jpj,1,2) )
          !
-         ! Get total input dust ; need to compute total atmospheric supply of N in a year
-         CALL iom_open ( TRIM( sn_ndepo%clname ), numdepo )
-         CALL iom_gettime( numdepo, zsteps, kntime=ntimes_ndep)
-         ALLOCATE( zndepo(jpi,jpj,ntimes_ndep) )
-         DO jm = 1, ntimes_ndep
-            CALL iom_get( numdepo, jpdom_data, TRIM( sn_ndepo%clvar ), zndepo(:,:,jm), jm )
-         END DO
-         CALL iom_close( numdepo )
-         ztimes_ndep = 1._wp / FLOAT( ntimes_ndep ) 
-         nitdepinput = 0._wp
-         DO jm = 1, ntimes_ndep
-           nitdepinput = nitdepinput + glob_sum( zndepo(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_ndep )
-         ENDDO
-         nitdepinput = nitdepinput / rno3 / 14E6 
-         DEALLOCATE( zndepo)
+         IF( Agrif_Root() ) THEN   !  Only on the master grid
+            ! Get total input dust ; need to compute total atmospheric supply of N in a year
+            CALL iom_open ( TRIM( sn_ndepo%clname ), numdepo )
+            CALL iom_gettime( numdepo, zsteps, kntime=ntimes_ndep)
+            ALLOCATE( zndepo(jpi,jpj,ntimes_ndep) )
+            DO jm = 1, ntimes_ndep
+               CALL iom_get( numdepo, jpdom_data, TRIM( sn_ndepo%clvar ), zndepo(:,:,jm), jm )
+            END DO
+            CALL iom_close( numdepo )
+            ztimes_ndep = 1._wp / FLOAT( ntimes_ndep ) 
+            nitdepinput = 0._wp
+            DO jm = 1, ntimes_ndep
+              nitdepinput = nitdepinput + glob_sum( zndepo(:,:,jm) * e1e2t(:,:) * tmask(:,:,1) * ztimes_ndep )
+            ENDDO
+            nitdepinput = nitdepinput / rno3 / 14E6 
+            DEALLOCATE( zndepo)
+         ENDIF
       ELSE
          nitdepinput = 0._wp