Changeset 3387
- Timestamp:
- Sep 4, 2018, 7:42:50 PM (6 years ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 2 edited
-
carbon_cycle_mod.F90 (modified) (8 diffs)
-
physiq_mod.F90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90
r3385 r3387 1 1 MODULE carbon_cycle_mod 2 ! Controle module for the carbon CO2 tracers : 2 !======================================================================= 3 ! Authors: Patricia Cadule and Laurent Fairhead 4 ! base sur un travail anterieur mene par Patricia Cadule et Josefine Ghattas 5 ! 6 ! Purpose and description: 7 ! ----------------------- 8 ! Control module for the carbon CO2 tracers : 3 9 ! - Identification 4 10 ! - Get concentrations comming from coupled model or read from file to tracers … … 6 12 ! - Calculate new carbon flux for sending to coupled models (PISCES and ORCHIDEE) 7 13 ! 8 ! Author : Josefine GHATTAS, Patricia CADULE 14 ! Module permettant de mettre a jour les champs (puits et sources) pour le 15 ! transport de CO2 en online (IPSL-CM et LMDZOR) et offline (lecture de carte) 16 ! 17 ! Le cas online/offline est defini par le flag carbon_cycle_cpl (y/n) 18 ! Le transport du traceur CO2 est defini par le flag carbon_cycle_tr (y/n) 19 ! la provenance des champs (termes de puits) est denini par le flag level_coupling_esm 20 ! 21 ! level_coupling_esm : level of coupling of the biogeochemical fields between 22 ! LMDZ, ORCHIDEE and NEMO 23 ! Definitions of level_coupling_esm in physiq.def 24 ! level_coupling_esm = 0 ! No field exchange between LMDZ and ORCHIDEE models 25 ! ! No field exchange between LMDZ and NEMO 26 ! level_coupling_esm = 1 ! Field exchange between LMDZ and ORCHIDEE models 27 ! ! No field exchange between LMDZ and NEMO models 28 ! level_coupling_esm = 2 ! No field exchange between LMDZ and ORCHIDEE models 29 ! ! Field exchange between LMDZ and NEMO models 30 ! level_coupling_esm = 3 ! Field exchange between LMDZ and ORCHIDEE models 31 ! ! Field exchange between LMDZ and NEMO models 32 !======================================================================= 9 33 10 34 IMPLICIT NONE 11 35 SAVE 12 36 PRIVATE 13 PUBLIC :: carbon_cycle_init, carbon_cycle 37 PUBLIC :: carbon_cycle_init, carbon_cycle, infocfields_init 14 38 15 39 ! Variables read from parmeter file physiq.def … … 18 42 LOGICAL, PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 19 43 !$OMP THREADPRIVATE(carbon_cycle_cpl) 20 INTEGER, SAVE,PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 344 INTEGER, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3 21 45 !$OMP THREADPRIVATE(level_coupling_esm) 22 46 … … 64 88 !$OMP THREADPRIVATE(co2_send) 65 89 90 ! nbfields : total number of fields 91 INTEGER, PUBLIC :: nbcf 92 !$OMP THREADPRIVATE(nbcf) 93 94 ! nbcf_in : number of fields IN 95 INTEGER, PUBLIC :: nbcf_in 96 !$OMP THREADPRIVATE(nbcf_in) 97 98 ! nbcf_in_orc : number of fields IN 99 INTEGER, PUBLIC :: nbcf_in_orc 100 !$OMP THREADPRIVATE(nbcf_in_orc) 101 102 ! nbcf_in_inca : number of fields IN (from INCA) 103 INTEGER, PUBLIC :: nbcf_in_inca 104 !$OMP THREADPRIVATE(nbcf_in_inca) 105 106 ! nbcf_in_nemo : number of fields IN (from nemo) 107 INTEGER, PUBLIC :: nbcf_in_nemo 108 !$OMP THREADPRIVATE(nbcf_in_nemo) 109 110 ! nbcf_in_ant : number of fields IN (from anthropogenic sources) 111 INTEGER, PUBLIC :: nbcf_in_ant 112 !$OMP THREADPRIVATE(nbcf_in_ant) 113 114 ! nbcf_out : number of fields OUT 115 INTEGER, PUBLIC :: nbcf_out 116 !$OMP THREADPRIVATE(nbcf_out) 117 118 ! Name of variables 119 CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname ! coupling field short name for restart (?) and diagnostics 120 !$OMP THREADPRIVATE(cfname) 121 122 CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname_in ! coupling field short name for restart (?) and diagnostics 123 !$OMP THREADPRIVATE(cfname_in) 124 125 CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname_out ! coupling field short name for restart (?) and diagnostics 126 !$OMP THREADPRIVATE(cfname_out) 127 128 CHARACTER(len=15), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfunits_in ! coupling field units for diagnostics 129 !$OMP THREADPRIVATE(cfunits_in) 130 131 CHARACTER(len=15), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfunits_out ! coupling field units for diagnostics 132 !$OMP THREADPRIVATE(cfunits_out) 133 134 CHARACTER(len=120), ALLOCATABLE, DIMENSION(:), PUBLIC :: cftext_in ! coupling field long name for diagnostics 135 !$OMP THREADPRIVATE(cftext_in) 136 137 CHARACTER(len=120), ALLOCATABLE, DIMENSION(:), PUBLIC :: cftext_out ! coupling field long name for diagnostics 138 !$OMP THREADPRIVATE(cftext_out) 139 140 CHARACTER(len=5), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfmod1 ! model 1 (rreference) : LMDz 141 !$OMP THREADPRIVATE(cfmod1) 142 143 CHARACTER(len=5), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfmod2 ! model 2 144 !$OMP THREADPRIVATE(cfmod2) 145 146 REAL, ALLOCATABLE, DIMENSION(:,:), PUBLIC :: zcfields_in 147 !$OMP THREADPRIVATE(zcfields_in) 148 149 REAL, ALLOCATABLE, DIMENSION(:,:), PUBLIC :: zcfields_out 150 !$OMP THREADPRIVATE(zcfields_out) 66 151 67 152 TYPE, PUBLIC :: co2_trac_type … … 111 196 INTEGER :: teststop 112 197 113 114 115 198 ! 1) Read controle parameters from .def input file 116 199 ! ------------------------------------------------ … … 125 208 WRITE(lunout,*) 'carbon_cycle_fos_fuel = ', fos_fuel_s 126 209 END IF 127 128 210 129 211 ! Read parmeter for calculation compatible emission … … 297 379 END SUBROUTINE carbon_cycle_init 298 380 299 300 381 SUBROUTINE carbon_cycle(nstep, pdtphys, pctsrf, tr_seri, source) 301 382 ! Subroutine for injection of co2 in the tracers … … 483 564 END SUBROUTINE carbon_cycle 484 565 566 SUBROUTINE infocfields_init 567 568 USE control_mod, ONLY: planet_type 569 USE phys_cal_mod, ONLY : mth_cur 570 USE mod_synchro_omp 571 USE mod_phys_lmdz_para, ONLY: is_mpi_root, is_omp_root 572 USE mod_phys_lmdz_transfert_para 573 USE mod_phys_lmdz_omp_transfert 574 USE dimphy, ONLY: klon 575 576 IMPLICIT NONE 577 578 !======================================================================= 579 ! 580 ! Authors: Patricia Cadule and Laurent Fairhead 581 ! ------- 582 ! 583 ! Purpose and description: 584 ! ----------------------- 585 ! 586 ! Infofields 587 ! this routine enables to define the field exchanges in both directions between 588 ! the atmospheric circulation model (LMDZ) and ORCHIDEE. In the future this 589 ! routing might apply to other models (e.g., NEMO, INCA, ...). 590 ! Therefore, currently with this routine, it is possible to define the coupling 591 ! fields only between LMDZ and ORCHIDEE. 592 ! The coupling_fields.def file enables to define the name of the exchanged 593 ! fields at the coupling interface. 594 ! field_in_names : the set of names of the exchanged fields in input to ORCHIDEE 595 ! (LMDZ to ORCHIDEE) 596 ! field_out_names : the set of names of the exchanged fields in output of 597 ! ORCHIDEE (ORCHIDEE to LMDZ) 598 ! n : the number of exchanged fields at th coupling interface 599 ! nb_fields_in : number of inputs fields to ORCHIDEE (LMDZ to ORCHIDEE) 600 ! nb_fields_out : number of ouput fields of ORCHIDEE (ORCHIDEE to LMDZ) 601 ! 602 ! The syntax for coupling_fields.def is as follows: 603 ! IMPORTANT: each column entry must be separated from the previous one by 3 604 ! spaces and only that 605 ! field name coupling model 1 model 2 long_name 606 ! direction 607 ! 10char -3spaces- 3char -3spaces- 4char -3spaces- 4char -3spaces- 30char 608 ! 609 ! n 610 ! FIELD1 IN LMDZ ORC 611 ! .... 612 ! FIELD(j) IN LMDZ ORC 613 ! FIELD(j+1) OUT LMDZ ORC 614 ! ... 615 ! FIELDn OUT LMDZ ORC 616 ! 617 !======================================================================= 618 ! ... 22/12/2017 .... 619 !----------------------------------------------------------------------- 620 ! Declarations 621 622 INCLUDE "clesphys.h" 623 INCLUDE "dimensions.h" 624 INCLUDE "iniprint.h" 625 626 ! Local variables 627 628 INTEGER :: iq, ierr, stat, error 629 630 CHARACTER(LEN=20), ALLOCATABLE, DIMENSION(:), SAVE :: cfname_root 631 CHARACTER(LEN=120), ALLOCATABLE, DIMENSION(:), SAVE :: cftext_root 632 CHARACTER(LEN=15), ALLOCATABLE, DIMENSION(:), SAVE :: cfunits_root 633 634 CHARACTER(len=3), ALLOCATABLE, DIMENSION(:) :: cfintent_root 635 CHARACTER(len=5), ALLOCATABLE, DIMENSION(:) :: cfmod1_root 636 CHARACTER(len=5), ALLOCATABLE, DIMENSION(:) :: cfmod2_root 637 638 LOGICAL, ALLOCATABLE, DIMENSION(:), SAVE :: mask_in_root 639 LOGICAL, ALLOCATABLE, DIMENSION(:), SAVE :: mask_out_root 640 641 CHARACTER(len=*),parameter :: modname="infocfields" 642 643 !----------------------------------------------------------------------- 644 645 nbcf=0 646 nbcf_in=0 647 nbcf_out=0 648 649 IF (planet_type=='earth') THEN 650 651 IF (is_mpi_root .AND. is_omp_root) THEN 652 653 IF (level_coupling_esm.GT.0) THEN 654 655 OPEN(200,file='coupling_fields.def',form='formatted',status='old', iostat=ierr) 656 657 IF (ierr.EQ.0) THEN 658 659 WRITE(lunout,*) trim(modname),': Open coupling_fields.def : ok' 660 READ(200,*) nbcf 661 WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf=',nbcf 662 IF (.NOT. ALLOCATED(cfname_root)) ALLOCATE(cfname_root(nbcf)) 663 IF (.NOT. ALLOCATED(cfintent_root)) ALLOCATE(cfintent_root(nbcf)) 664 IF (.NOT. ALLOCATED(cfmod1_root)) ALLOCATE(cfmod1_root(nbcf)) 665 IF (.NOT. ALLOCATED(cfmod2_root)) ALLOCATE(cfmod2_root(nbcf)) 666 IF (.NOT. ALLOCATED(cftext_root)) ALLOCATE(cftext_root(nbcf)) 667 IF (.NOT. ALLOCATED(cfunits_root)) ALLOCATE(cfunits_root(nbcf)) 668 IF (.NOT. ALLOCATED(mask_in_root)) ALLOCATE(mask_in_root(nbcf)) 669 IF (.NOT. ALLOCATED(mask_out_root)) ALLOCATE(mask_out_root(nbcf)) 670 671 nbcf_in=0 672 nbcf_out=0 673 DO iq=1,nbcf 674 WRITE(lunout,*) 'infofields : field=',iq 675 READ(200,'(A15,3X,A3,3X,A5,3X,A5,3X,A120,3X,A15)',IOSTAT=ierr) & 676 cfname_root(iq),cfintent_root(iq),cfmod1_root(iq),cfmod2_root(iq),cftext_root(iq),cfunits_root(iq) 677 cfname_root(iq)=TRIM(cfname_root(iq)) 678 cfintent_root(iq)=TRIM(cfintent_root(iq)) 679 cfmod1_root(iq)=TRIM(cfmod1_root(iq)) 680 cfmod2_root(iq)=TRIM(cfmod2_root(iq)) 681 cftext_root(iq)=TRIM(cftext_root(iq)) 682 cfunits_root(iq)=TRIM(cfunits_root(iq)) 683 WRITE(lunout,*) 'coupling field: ',cfname_root(iq), & 684 ', number: ',iq,', INTENT: ',cfintent_root(iq) 685 WRITE(lunout,*) 'coupling field: ',cfname_root(iq), & 686 ', number: ',iq,', model 1 (ref): ',cfmod1_root(iq),', model 2: ',cfmod2_root(iq) 687 WRITE(lunout,*) 'coupling field: ',cfname_root(iq), & 688 ', number: ',iq,', long name: ',cftext_root(iq),', units ',cfunits_root(iq) 689 IF ((TRIM(cfintent_root(iq)).NE.'OUT') .AND. (nbcf_in.LE.nbcf)) THEN 690 nbcf_in=nbcf_in+1 691 mask_in_root(iq)=.TRUE. 692 mask_out_root(iq)=.FALSE. 693 ELSE IF ((TRIM(cfintent_root(iq)).EQ.'OUT') .AND. (nbcf_out.LE.nbcf)) THEN 694 nbcf_out=nbcf_out+1 695 mask_in_root(iq)=.FALSE. 696 mask_out_root(iq)=.TRUE. 697 ELSE 698 WRITE(lunout,*) 'abort_gcm --- nbcf: ',nbcf 699 WRITE(lunout,*) 'abort_gcm --- nbcf_in: ',nbcf_in 700 WRITE(lunout,*) 'abort_gcm --- nbcf_out: ',nbcf_out 701 CALL abort_gcm('infocfields_init','Problem in the definition of the coupling fields',1) 702 ENDIF 703 ENDDO !DO iq=1,nbcf 704 ELSE 705 WRITE(lunout,*) trim(modname),': infocfields_mod.F90 --- Problem in opening coupling_fields.def' 706 WRITE(lunout,*) trim(modname),': infocfields_mod.F90 --- WARNING using defaut values' 707 ENDIF ! ierr 708 CLOSE(200) 709 710 ENDIF ! level_coupling_esm 711 712 IF (nbcf_out .EQ. 0) nbcf_out=-1 713 IF (nbcf_in .EQ. 0) nbcf_in=-1 714 715 ENDIF ! (is_mpi_root .AND. is_omp_root) 716 !$OMP BARRIER 717 718 CALL bcast(nbcf) 719 CALL bcast(nbcf_in) 720 CALL bcast(nbcf_out) 721 722 WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf=',nbcf 723 WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf_in=',nbcf_in 724 WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf_out=',nbcf_out 725 726 IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfname))) ALLOCATE(cfname(nbcf)) 727 IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cfname_in))) ALLOCATE(cfname_in(nbcf_in)) 728 IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cftext_in))) ALLOCATE(cftext_in(nbcf_in)) 729 IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cfname_out))) ALLOCATE(cfname_out(nbcf_out)) 730 IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cftext_out))) ALLOCATE(cftext_out(nbcf_out)) 731 IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfmod1))) ALLOCATE(cfmod1(nbcf)) 732 IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfmod2))) ALLOCATE(cfmod2(nbcf)) 733 IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cfunits_in))) ALLOCATE(cfunits_in(nbcf_in)) 734 IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cfunits_out))) ALLOCATE(cfunits_out(nbcf_out)) 735 736 IF (is_mpi_root .AND. is_omp_root) THEN 737 738 IF (nbcf.GT.0) cfname=cfname_root 739 IF (nbcf_in.GT.0) cfname_in=PACK(cfname_root,mask_in_root) 740 IF (nbcf_out.GT.0) cfname_out=PACK(cfname_root,mask_out_root) 741 IF (nbcf_in.GT.0) cftext_in=PACK(cftext_root,mask_in_root) 742 IF (nbcf_out.GT.0) cftext_out=PACK(cftext_root,mask_out_root) 743 IF (nbcf.GT.0) cfmod1=cfmod1_root 744 IF (nbcf.GT.0) cfmod2=cfmod2_root 745 IF (nbcf_in.GT.0) cfunits_in=PACK(cfunits_root,mask_in_root) 746 IF (nbcf_out.GT.0) cfunits_out=PACK(cfunits_root,mask_out_root) 747 748 nbcf_in_orc=0 749 nbcf_in_nemo=0 750 nbcf_in_inca=0 751 nbcf_in_ant=0 752 753 DO iq=1,nbcf 754 IF (cfmod1(iq) == "ORC") nbcf_in_orc=nbcf_in_orc+1 755 IF (cfmod1(iq) == "NEMO") nbcf_in_nemo=nbcf_in_nemo+1 756 IF (cfmod1(iq) == "INCA") nbcf_in_inca=nbcf_in_inca+1 757 IF (cfmod1(iq) == "ALL") nbcf_in_orc=nbcf_in_orc+1 ! ALL = ORC/NEMO/INCA 758 IF (cfmod1(iq) == "ALL") nbcf_in_nemo=nbcf_in_nemo+1 ! ALL = ORC/NEMO/INCA 759 IF (cfmod1(iq) == "ALL") nbcf_in_inca=nbcf_in_inca+1 ! ALL = ORC/NEMO/INCA 760 IF (cfmod1(iq) == "ANT") nbcf_in_ant=nbcf_in_ant+1 761 ENDDO 762 763 ENDIF ! (is_mpi_root .AND. is_omp_root) 764 !$OMP BARRIER 765 766 CALL bcast(nbcf_in_orc) 767 CALL bcast(nbcf_in_nemo) 768 CALL bcast(nbcf_in_inca) 769 CALL bcast(nbcf_in_ant) 770 771 WRITE(lunout,*) 'nbcf_in_orc ',nbcf_in_orc 772 WRITE(lunout,*) 'nbcf_in_nemo ',nbcf_in_nemo 773 WRITE(lunout,*) 'nbcf_in_inca ',nbcf_in_inca 774 WRITE(lunout,*) 'nbcf_in_ant ',nbcf_in_ant 775 776 IF (nbcf_in.GT.0) THEN 777 DO iq=1,nbcf_in 778 CALL bcast(cfname_in(iq)) 779 CALL bcast(cftext_in(iq)) 780 CALL bcast(cfunits_in(iq)) 781 ENDDO 782 ENDIF 783 784 IF (nbcf_out.GT.0) THEN 785 DO iq=1,nbcf_out 786 CALL bcast(cfname_out(iq)) 787 CALL bcast(cftext_out(iq)) 788 CALL bcast(cfunits_out(iq)) 789 ENDDO 790 ENDIF 791 792 IF (nbcf.GT.0) THEN 793 DO iq=1,nbcf 794 CALL bcast(cfmod1(iq)) 795 CALL bcast(cfmod2(iq)) 796 ENDDO 797 ENDIF 798 799 IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfields_mod --- cfname_in: ',cfname_in 800 IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfields_mod --- cfname_out: ',cfname_out 801 802 IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfields_mod --- cftext_in: ',cftext_in 803 IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfields_mod --- cftext_out: ',cftext_out 804 805 IF (nbcf.GT.0) WRITE(lunout,*)'infocfields_mod --- cfmod1: ',cfmod1 806 IF (nbcf.GT.0) WRITE(lunout,*)'infocfields_mod --- cfmod2: ',cfmod2 807 808 IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfunits_mod --- cfunits_in: ',cfunits_in 809 IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfunits_mod --- cfunits_out: ',cfunits_out 810 811 IF (nbcf_in.GT.0) WRITE(*,*)'infocfields_init --- number of fields in to LMDZ: ',nbcf_in 812 IF (nbcf_out.GT.0) WRITE(*,*)'infocfields_init --- number of fields out of LMDZ: ',nbcf_out 813 814 ELSE 815 ! Default values for other planets 816 nbcf=0 817 nbcf_in=0 818 nbcf_out=0 819 ENDIF ! planet_type 820 821 IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(zcfields_in))) ALLOCATE(zcfields_in(klon,nbcf_in),stat=error) 822 IF (error /= 0) CALL abort_gcm(modname,'Pb in allocation zcfields_in',1) 823 IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(zcfields_out))) ALLOCATE(zcfields_out(klon,nbcf_out),stat=error) 824 IF (error /= 0) CALL abort_gcm(modname,'Pb in allocation zcfields_out',1) 825 826 END SUBROUTINE infocfields_init 827 485 828 END MODULE carbon_cycle_mod -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r3370 r3387 16 16 d_u, d_v, d_t, d_qx, d_ps) 17 17 18 useassert_m, only: assert18 USE assert_m, only: assert 19 19 USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, & 20 20 histwrite, ju2ymds, ymds2ju, getin … … 168 168 omega, & 169 169 epmax_diag, & 170 171 170 ! Deep convective variables used in phytrac 172 171 pmflxr, pmflxs, & … … 184 183 epmlmMm, eplaMm, & 185 184 sij, & 186 185 ! 187 186 cldemi, & 188 187 cldfra, cldtau, fiwc, & … … 199 198 prfl, psfl, fraca, Vprecip, & 200 199 zw2, & 201 200 ! 202 201 fluxu, fluxv, & 203 202 fluxt, & 204 203 ! 205 204 uwriteSTD, vwriteSTD, & !pour calcul_STDlev.h 206 205 wwriteSTD, phiwriteSTD, & !pour calcul_STDlev.h 207 206 qwriteSTD, twriteSTD, rhwriteSTD, & !pour calcul_STDlev.h 208 207 ! 209 208 beta_prec, & 210 209 rneb, & … … 213 212 USE phys_state_var_mod ! Variables sauvegardees de la physique 214 213 #ifdef CPP_Dust 215 USE phys_output_write_spl_mod214 USE phys_output_write_spl_mod 216 215 #else 217 216 USE phys_output_var_mod ! Variables pour les ecritures des sorties … … 222 221 USE phys_output_mod 223 222 USE phys_output_ctrlout_mod 224 useopen_climoz_m, only: open_climoz ! ozone climatology from a file225 useregr_pr_time_av_m, only: regr_pr_time_av226 usenetcdf95, only: nf95_close223 USE open_climoz_m, only: open_climoz ! ozone climatology from a file 224 USE regr_pr_time_av_m, only: regr_pr_time_av 225 USE netcdf95, only: nf95_close 227 226 !IM for NMC files 228 ! use netcdf, only: nf90_fill_real 229 use netcdf, only: nf90_fill_real 230 use mod_phys_lmdz_mpi_data, only: is_mpi_root 227 USE netcdf, only: nf90_fill_real 228 USE mod_phys_lmdz_mpi_data, only: is_mpi_root 231 229 USE aero_mod 232 useozonecm_m, only: ozonecm ! ozone of J.-F. Royer233 useconf_phys_m, only: conf_phys234 useradlwsw_m, only: radlwsw235 usephyaqua_mod, only: zenang_an230 USE ozonecm_m, only: ozonecm ! ozone of J.-F. Royer 231 USE conf_phys_m, only: conf_phys 232 USE radlwsw_m, only: radlwsw 233 USE phyaqua_mod, only: zenang_an 236 234 USE time_phylmdz_mod, only: day_step_phy, annee_ref, day_ref, itau_phy, & 237 235 start_time, pdtphys, day_ini … … 246 244 USE indice_sol_mod 247 245 USE phytrac_mod, ONLY : phytrac 246 USE carbon_cycle_mod, ONLY : infocfields_init 248 247 249 248 #ifdef CPP_RRTM … … 263 262 !IM stations CFMIP 264 263 USE CFMIP_point_locations 265 useFLOTT_GWD_rando_m, only: FLOTT_GWD_rando266 useACAMA_GWD_rando_m, only: ACAMA_GWD_rando264 USE FLOTT_GWD_rando_m, only: FLOTT_GWD_rando 265 USE ACAMA_GWD_rando_m, only: ACAMA_GWD_rando 267 266 USE VERTICAL_LAYERS_MOD, ONLY: aps,bps 268 267 269 270 IMPLICIT none 268 IMPLICIT NONE 271 269 !>====================================================================== 272 270 !! … … 1305 1303 nvm_lmdz = 13 1306 1304 CALL getin_p('NVM',nvm_lmdz) 1305 1306 !--PC: defining fields to be exchanged between LMDz, ORCHIDEE and NEMO 1307 WRITE(lunout,*) 'Call to infocfields from physiq' 1308 CALL infocfields_init 1309 1307 1310 ENDIF 1308 1311 1309 1312 IF (prt_level.ge.1) print *,'CONVERGENCE PHYSIQUE THERM 1 ' 1310 1311 1313 1312 1314 !======================================================================
Note: See TracChangeset
for help on using the changeset viewer.
