Changeset 2180 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Jan 16, 2015, 2:15:03 PM (10 years ago)
Author:
acozic
Message:

After checking I reload commit 2169 and 2170

Location:
LMDZ5/trunk/libf
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3d_common/infotrac.F90

    r2171 r2180  
    2929
    3030  CHARACTER(len=4),SAVE :: type_trac
     31  CHARACTER(len=8),DIMENSION(:),ALLOCATABLE, SAVE :: solsym
    3132 
    3233CONTAINS
    … …  
    6263
    6364    CHARACTER(len=15), ALLOCATABLE, DIMENSION(:) :: tnom_0  ! tracer short name
    64     CHARACTER(len=8), ALLOCATABLE, DIMENSION(:) :: tracnam ! name from INCA
    6565    CHARACTER(len=3), DIMENSION(30) :: descrq
    6666    CHARACTER(len=1), DIMENSION(3)  :: txts
    … …  
    9494       WRITE(lunout,*) 'You have choosen to couple with INCA chemestry model : type_trac=', &
    9595            type_trac,' config_inca=',config_inca
    96        IF (config_inca/='aero' .AND. config_inca/='chem') THEN
     96       IF (config_inca/='aero' .AND. config_inca/='aeNP' .AND. config_inca/='chem') THEN
    9797          WRITE(lunout,*) 'Incoherence between type_trac and config_inca. Model stops. Modify run.def'
    9898          CALL abort_gcm('infotrac_init','Incoherence between type_trac and config_inca',1)
    … …  
    172172!
    173173    ALLOCATE(tnom_0(nqtrue), hadv(nqtrue), vadv(nqtrue))
    174     ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(nbtr))
     174    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), solsym(nbtr))
    175175    conv_flg(:) = 1 ! convection activated for all tracers
    176176    pbl_flg(:)  = 1 ! boundary layer activated for all tracers
    … …  
    254254            conv_flg, &
    255255            pbl_flg,  &
    256             tracnam)
     256            solsym)
    257257#endif
    258258       tnom_0(1)='H2Ov'
    … …  
    260260
    261261       DO iq =3,nqtrue
    262           tnom_0(iq)=tracnam(iq-2)
     262          tnom_0(iq)=solsym(iq-2)
    263263       END DO
    264264       nqo = 2
    … …  
    394394!
    395395    DEALLOCATE(tnom_0, hadv, vadv)
    396     DEALLOCATE(tracnam)
     396
    397397
    398398  END SUBROUTINE infotrac_init

  • LMDZ5/trunk/libf/dyn3dmem/gcm.F

    r2171 r2180  
    240240     $        iphysiq,day_step,nday,
    241241     $        nbsrf, is_oce,is_sic,
    242      $        is_ter,is_lic)
     242     $        is_ter,is_lic, calend)
    243243
    244244         call init_inca_para(

  • LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F

    r2171 r2180  
    16041604c$OMP MASTER
    16051605                 call fin_getparam
     1606c$OMP END MASTER
     1607
     1608#ifdef INCA
     1609                 call finalize_inca
     1610#endif
     1611
     1612c$OMP MASTER
    16061613                 call finalize_parallel
    16071614c$OMP END MASTER

  • LMDZ5/trunk/libf/dyn3dpar/gcm.F

    r2171 r2180  
    246246     $        iphysiq,day_step,nday,
    247247     $        nbsrf, is_oce,is_sic,
    248      $        is_ter,is_lic)
     248     $        is_ter,is_lic, calend)
    249249
    250250         call init_inca_para(

  • LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F

    r2171 r2180  
    14101410c$OMP MASTER
    14111411              call fin_getparam
    1412               call finalize_parallel
     1412c$OMP END MASTER
     1413#ifdef INCA
     1414                 call finalize_inca
     1415#endif
     1416c$OMP MASTER
     1417               call finalize_parallel
    14131418c$OMP END MASTER
    14141419              abort_message = 'Simulation finished'

  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2159 r2180  
    228228    USE ocean_slab_mod, only: tslab, slab_bils
    229229    USE indice_sol_mod, only: nbsrf
    230     USE infotrac, only: nqtot, nqo
     230    USE infotrac, only: nqtot, nqo, type_trac
    231231    USE comgeomphy, only: airephy
    232232    USE surface_data, only: type_ocean, ok_veget, ok_snow
    … …  
    13231323        IF (nqtot.GE.nqo+1) THEN
    13241324            DO iq=nqo+1,nqtot
     1325              IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN
    13251326
    13261327             CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
    … …  
    13451346             ENDIF
    13461347             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
     1348             endif
    13471349          ENDDO
    13481350       ENDIF

  • LMDZ5/trunk/libf/phylmd/physiq.F90

    r2171 r2180  
    13081308             annee_ref, &
    13091309             day_ref,  &
    1310              itau_phy)
     1310             itau_phy, &
     1311             io_lon, &
     1312             io_lat)
    13111313
    13121314        CALL VTe(VTinca)

  • LMDZ5/trunk/libf/phylmd/phytrac_mod.F90

    r2171 r2180  
    308308!$OMP THREADPRIVATE(lessivage)
    309309
    310     CHARACTER(len=8),DIMENSION(nbtr) :: solsym
    311310    !RomP >>>
    312311    INTEGER,SAVE  :: iflag_lscav_omp,iflag_lscav
    … …  
    557556            cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon,iflag_vdf_trac>=0,sh, &
    558557            rh, pphi, ustar, wstar, ale_bl, ale_wake,  u10m, v10m, &
    559             tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
     558            tr_seri, source, d_tr_cl,d_tr_dec, zmasse)               !RomP
    560559
    561560    CASE('inca')
    … …  
    572571            tau_aero, piz_aero, cg_aero,        ccm,       &
    573572            rfname,                                        &
    574             tr_seri,  source,   solsym)     
     573            tr_seri,  source)     
    575574
    576575    CASE('repr')
    … …  
    580579            presnivs, xlat, xlon, pphis, pphi, &
    581580            t_seri, pplay, paprs, sh , &
    582             tr_seri, solsym)
     581            tr_seri)
    583582
    584583    END SELECT

  • LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

    r2171 r2180  
    8484  REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer
    8585  REAL, DIMENSION(klon,klev,naero_tot) :: m_allaer_pi !RAF 
    86   !  REAL, DIMENSION(klon,naero_tot)      :: fractnat_allaer !RAF delete??
    87   character(len=8), dimension(nbtr) :: tracname
     86
    8887  integer :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
    8988  integer :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
    … …  
    104103     !--convert to ug m-3 unit for consistency with offline fields
    105104     !
    106 #ifdef INCA
    107      call tracinca_name(tracname)
    108 #endif
    109 
    110105     do i=1,nbtr
    111         select case(trim(tracname(i)))
     106        select case(trim(solsym(i)))
    112107           case ("ASBCM")
    113108              id_ASBCM = i

  • LMDZ5/trunk/libf/phylmd/tracinca_mod.F90

    r2171 r2180  
    3535       tau_aero, piz_aero, cg_aero,        ccm,       &
    3636       rfname,                                        &
    37        tr_seri,  source,   solsym)     
     37       tr_seri,  source)     
    3838
    3939!========================================================
    … …  
    113113  ! Output arguments
    114114    REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit
    115     CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
    116115
    117116!=======================================================================================
    … …  
    134133       pdel(:,k) = paprs(:,k) - paprs (:,k+1)
    135134    END DO
    136    
    137     zpmfu(:,:)=pmfu(:,:)
     135 
     136#ifdef INCA
     137    IF (config_inca == 'aero') THEN
     138       zpmfu(:,:)=pmfu(:,:)       
     139    ELSE IF (config_inca == 'aeNP') THEN
     140       zpmfu(:,:)=upwd(:,:)
     141    ENDIF
    138142
    139     IF (config_inca == 'aero') THEN
    140 #ifdef INCA
    141        CALL aerosolmain(                    &
    142             aerosol_couple,tr_seri,pdtphys, &
    143             pplay,pdel,prfl,pmflxr,psfl,    &
    144             pmflxs,zpmfu,itop_con,ibas_con,  &
    145             pphi,airephy,nstep,rneb,t_seri, &     
    146             rh,tau_aero,piz_aero,cg_aero,   &
    147             rfname,ccm,lafin)
     143    CALL aerosolmain(                    &
     144         aerosol_couple,tr_seri,pdtphys, &
     145         pplay,pdel,prfl,pmflxr,psfl,    &
     146         pmflxs,zpmfu,itop_con,ibas_con,  &
     147         pphi,airephy,nstep,rneb,t_seri, &     
     148         rh,tau_aero,piz_aero,cg_aero,   &
     149         rfname,ccm,lafin, config_inca)
    148150#endif
    149     END IF
    150 
    151     IF (config_inca == 'aeNP') THEN
    152 #ifdef INCA
    153        zpmfu(:,:)=upwd(:,:)
    154        CALL aerosolmainNP(                  &
    155             aerosol_couple,tr_seri,pdtphys, &
    156             pplay,pdel,prfl,pmflxr,psfl,    &
    157             pmflxs,zpmfu,itop_con,ibas_con,  &
    158             pphi,airephy,nstep,rneb,t_seri, &     
    159             rh,lafin)
    160 #endif
    161     END IF
    162151
    163152
    … …  
    196185         iip1,       & !nx
    197186         jjp1,       & !ny
    198          source,     &
    199          solsym)
     187         source )
    200188#endif
    201189   

  • LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

    r2171 r2180  
    338338       cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon, couchelimite, sh, &
    339339       rh, pphi, ustar, wstar, ale_bl, ale_wake,  zu10m, zv10m, &
    340 !!          tr_seri, source, solsym, d_tr_cl, zmasse)                      !RomP
    341           tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
     340       tr_seri, source, d_tr_cl,d_tr_dec, zmasse)               !RomP
    342341   
    343342    USE dimphy
    … …  
    397396
    398397! Output argument
    399     CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
    400398    REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit
    401399    REAL,DIMENSION(klon,klev,nbtr), INTENT(OUT)   :: d_tr_cl ! Td couche limite/traceur

  • LMDZ5/trunk/libf/phylmd/tracreprobus_mod.F90

    r2171 r2180  
    99       presnivs, xlat, xlon, pphis, pphi, &
    1010       t_seri, pplay, paprs, sh , &
    11        tr_seri, solsym)
     11       tr_seri)
    1212
    1313    USE dimphy
    … …  
    4242!----------------
    4343    REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT)  :: tr_seri ! Concentration Traceur [U/KgA] 
    44     CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
    4544 
    4645
Note: See TracChangeset for help on using the changeset viewer.