Changeset 3814 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/carbon_cycle_mod.F90

Timestamp:

Apr 15, 2015, 6:49:07 PM (10 years ago)

Author:

ymipsl

Message:

remove all dynamic dependency in LMDZ physics except for the include "dimensions.h"

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/carbon_cycle_mod.F90 (modified) (10 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/carbon_cycle_mod.F90

@@ -88 +88 @@
     USE comgeomphy
     USE mod_phys_lmdz_transfert_para
-    USE infotrac
+    USE infotrac_phy
     USE IOIPSL
     USE surface_data, ONLY : ok_veget, type_ocean
@@ -135 +135 @@
        WRITE(lunout,*) 'carbon_cycle_emis_comp = ',carbon_cycle_emis_comp
        IF (carbon_cycle_emis_comp) THEN
-          CALL abort_gcm('carbon_cycle_init', 'carbon_cycle_emis_comp option not yet implemented!!',1)
+          CALL abort_physic('carbon_cycle_init', 'carbon_cycle_emis_comp option not yet implemented!!',1)
        END IF
     END IF
@@ -225 +225 @@
           co2trac(itc)%updatefreq = 86400
           ! DOES THIS WORK ???? Problematic due to implementation of the coupled fluxes...
-          CALL abort_gcm('carbon_cycle_init','transport of total CO2 has to be implemented and tested',1)
+          CALL abort_physic('carbon_cycle_init','transport of total CO2 has to be implemented and tested',1)
        END SELECT
     END DO
@@ -247 +247 @@
     ! Allocate vector for storing fluxes to inject
     ALLOCATE(dtr_add(klon,maxco2trac), stat=ierr)
-    IF (ierr /= 0) CALL abort_gcm('carbon_cycle_init', 'pb in allocation 11',1)       
+    IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 11',1)       
     
     ! Allocate variables for cumulating fluxes from ORCHIDEE
@@ -253 +253 @@
        IF (.NOT. carbon_cycle_tr .AND. carbon_cycle_cpl) THEN
           ALLOCATE(fco2_land_day(klon), stat=ierr)
-          IF (ierr /= 0) CALL abort_gcm('carbon_cycle_init', 'pb in allocation 2',1)
+          IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 2',1)
           fco2_land_day(1:klon) = 0.
           
           ALLOCATE(fco2_lu_day(klon), stat=ierr)
-          IF (ierr /= 0) CALL abort_gcm('carbon_cycle_init', 'pb in allocation 3',1)
+          IF (ierr /= 0) CALL abort_physic('carbon_cycle_init', 'pb in allocation 3',1)
           fco2_lu_day(1:klon)   = 0.
        END IF
@@ -267 +267 @@
 !    IF (carbon_cycle_cpl .AND. type_ocean/='couple') THEN
 !       WRITE(lunout,*) 'Coupling with ocean model is needed for carbon_cycle_cpl'
-!       CALL abort_gcm('carbon_cycle_init', 'coupled ocean is needed for carbon_cycle_cpl',1)
+!       CALL abort_physic('carbon_cycle_init', 'coupled ocean is needed for carbon_cycle_cpl',1)
 !    END IF
 !
 !    IF (carbon_cycle_cpl .AND..NOT. ok_veget) THEN
 !       WRITE(lunout,*) 'Coupling with surface land model ORCHDIEE is needed for carbon_cycle_cpl'
-!       CALL abort_gcm('carbon_cycle_init', 'ok_veget is needed for carbon_cycle_cpl',1)
+!       CALL abort_physic('carbon_cycle_init', 'ok_veget is needed for carbon_cycle_cpl',1)
 !    END IF
 
@@ -278 +278 @@
     teststop=0
     DO it=1,teststop
-       CALL abort_gcm('carbon_cycle_init', 'Entering loop from 1 to 0',1)
+       CALL abort_physic('carbon_cycle_init', 'Entering loop from 1 to 0',1)
     END DO
 
@@ -284 +284 @@
        ! No carbon tracers found in tracer.def. It is not possible to do carbon cycle 
        WRITE(lunout,*) 'No carbon tracers found in tracer.def. Not ok with carbon_cycle_tr and/or carbon_cycle_cp'
-       CALL abort_gcm('carbon_cycle_init', 'No carbon tracers found in tracer.def',1)
+       CALL abort_physic('carbon_cycle_init', 'No carbon tracers found in tracer.def',1)
     END IF
     
@@ -303 +303 @@
 ! - Calculate CO2 flux to send to ocean and land models (PISCES and ORCHIDEE)
 
-    USE infotrac
+    USE infotrac_phy
     USE dimphy
     USE mod_phys_lmdz_transfert_para
@@ -370 +370 @@
              CASE DEFAULT
                 WRITE(lunout,*) 'Error with tracer ',co2trac(it)%name
-                CALL abort_gcm('carbon_cycle', 'No coupling implemented for this tracer',1)
+                CALL abort_physic('carbon_cycle', 'No coupling implemented for this tracer',1)
              END SELECT
           ELSE