Changeset 2171 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Dec 19, 2014, 4:21:08 PM (10 years ago)
Author:
acozic
Message:

There are some commits that we must not do just before holiday .... so be back to rev 2168

Location:
LMDZ5/trunk/libf
Files:
11 edited

Legend:

Unmodified
Added
Removed

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

    r2169 r2171  
    2929
    3030  CHARACTER(len=4),SAVE :: type_trac
    31   CHARACTER(len=8),DIMENSION(:),ALLOCATABLE, SAVE :: solsym
    3231 
    3332CONTAINS
    … …  
    6362
    6463    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/='aeNP' .AND. config_inca/='chem') THEN
     96       IF (config_inca/='aero' .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), solsym(nbtr))
     174    ALLOCATE(conv_flg(nbtr), pbl_flg(nbtr), tracnam(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             solsym)
     256            tracnam)
    257257#endif
    258258       tnom_0(1)='H2Ov'
    … …  
    260260
    261261       DO iq =3,nqtrue
    262           tnom_0(iq)=solsym(iq-2)
     262          tnom_0(iq)=tracnam(iq-2)
    263263       END DO
    264264       nqo = 2
    … …  
    394394!
    395395    DEALLOCATE(tnom_0, hadv, vadv)
    396 
     396    DEALLOCATE(tracnam)
    397397
    398398  END SUBROUTINE infotrac_init

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

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

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

    r2170 r2171  
    16041604c$OMP MASTER
    16051605                 call fin_getparam
    1606 c$OMP END MASTER
    1607 
    1608 #ifdef INCA
    1609                  call finalize_inca
    1610 #endif
    1611 
    1612 c$OMP MASTER
    16131606                 call finalize_parallel
    16141607c$OMP END MASTER

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

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

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

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

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

    r2170 r2171  
    13081308             annee_ref, &
    13091309             day_ref,  &
    1310              itau_phy, &
    1311              io_lon, &
    1312              io_lat)
     1310             itau_phy)
    13131311
    13141312        CALL VTe(VTinca)

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

    r2169 r2171  
    308308!$OMP THREADPRIVATE(lessivage)
    309309
     310    CHARACTER(len=8),DIMENSION(nbtr) :: solsym
    310311    !RomP >>>
    311312    INTEGER,SAVE  :: iflag_lscav_omp,iflag_lscav
    … …  
    556557            cdragh,  coefh, yu1, yv1, ftsol, pctsrf, xlat, xlon,iflag_vdf_trac>=0,sh, &
    557558            rh, pphi, ustar, wstar, ale_bl, ale_wake,  u10m, v10m, &
    558             tr_seri, source, d_tr_cl,d_tr_dec, zmasse)               !RomP
     559            tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
    559560
    560561    CASE('inca')
    … …  
    571572            tau_aero, piz_aero, cg_aero,        ccm,       &
    572573            rfname,                                        &
    573             tr_seri,  source)     
     574            tr_seri,  source,   solsym)     
    574575
    575576    CASE('repr')
    … …  
    579580            presnivs, xlat, xlon, pphis, pphi, &
    580581            t_seri, pplay, paprs, sh , &
    581             tr_seri)
     582            tr_seri, solsym)
    582583
    583584    END SELECT

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

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

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

    r2169 r2171  
    3535       tau_aero, piz_aero, cg_aero,        ccm,       &
    3636       rfname,                                        &
    37        tr_seri,  source)     
     37       tr_seri,  source,   solsym)     
    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
    115116
    116117!=======================================================================================
    … …  
    133134       pdel(:,k) = paprs(:,k) - paprs (:,k+1)
    134135    END DO
    135  
     136   
     137    zpmfu(:,:)=pmfu(:,:)
     138
     139    IF (config_inca == 'aero') THEN
    136140#ifdef INCA
    137     IF (config_inca == 'aero') THEN
    138        zpmfu(:,:)=pmfu(:,:)       
    139     ELSE IF (config_inca == 'aeNP') THEN
     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)
     148#endif
     149    END IF
     150
     151    IF (config_inca == 'aeNP') THEN
     152#ifdef INCA
    140153       zpmfu(:,:)=upwd(:,:)
    141     ENDIF
    142 
    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)
     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)
    150160#endif
     161    END IF
    151162
    152163
    … …  
    185196         iip1,       & !nx
    186197         jjp1,       & !ny
    187          source )
     198         source,     &
     199         solsym)
    188200#endif
    189201   

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

    r2169 r2171  
    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, d_tr_cl,d_tr_dec, zmasse)               !RomP
     340!!          tr_seri, source, solsym, d_tr_cl, zmasse)                      !RomP
     341          tr_seri, source, solsym, d_tr_cl,d_tr_dec, zmasse)               !RomP
    341342   
    342343    USE dimphy
    … …  
    396397
    397398! Output argument
     399    CHARACTER(len=8),DIMENSION(nbtr), INTENT(OUT) :: solsym
    398400    REAL,DIMENSION(klon,nbtr), INTENT(OUT)        :: source  ! a voir lorsque le flux de surface est prescrit
    399401    REAL,DIMENSION(klon,klev,nbtr), INTENT(OUT)   :: d_tr_cl ! Td couche limite/traceur

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

    r2169 r2171  
    99       presnivs, xlat, xlon, pphis, pphi, &
    1010       t_seri, pplay, paprs, sh , &
    11        tr_seri)
     11       tr_seri, solsym)
    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
    4445 
    4546
Note: See TracChangeset for help on using the changeset viewer.