Ignore:
Timestamp:
Jul 23, 2024, 3:29:36 PM (8 weeks ago)
Author:
abarral
Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_flux2.F90

    r5101 r5103  
    8686! Verification de la nullite des entrainement et detrainement au dessus
    8787! de lmax(ig)
    88 ! Active uniquement si check_debug=.true. ou prt_level>=10
     88! Active uniquement si check_debug=.TRUE. ou prt_level>=10
    8989!-------------------------------------------------------------------------
    9090
    … …  
    9696            if (l<=lmax(ig)) then
    9797               if (entr_star(ig,l)>1.) then
    98                     print*,'WARNING thermcell_flux 1 ig,l,lmax(ig)',ig,l,lmax(ig)
    99                     print*,'entr_star(ig,l)',entr_star(ig,l)
    100                     print*,'alim_star(ig,l)',alim_star(ig,l)
    101                     print*,'detr_star(ig,l)',detr_star(ig,l)
     98                    PRINT*,'WARNING thermcell_flux 1 ig,l,lmax(ig)',ig,l,lmax(ig)
     99                    PRINT*,'entr_star(ig,l)',entr_star(ig,l)
     100                    PRINT*,'alim_star(ig,l)',alim_star(ig,l)
     101                    PRINT*,'detr_star(ig,l)',detr_star(ig,l)
    102102               endif
    103103            else
    104104               if (abs(entr_star(ig,l))+abs(alim_star(ig,l))+abs(detr_star(ig,l))>0.) then
    105                     print*,'cas 1 : ig,l,lmax(ig)',ig,l,lmax(ig)
    106                     print*,'entr_star(ig,l)',entr_star(ig,l)
    107                     print*,'alim_star(ig,l)',alim_star(ig,l)
    108                     print*,'detr_star(ig,l)',detr_star(ig,l)
     105                    PRINT*,'cas 1 : ig,l,lmax(ig)',ig,l,lmax(ig)
     106                    PRINT*,'entr_star(ig,l)',entr_star(ig,l)
     107                    PRINT*,'alim_star(ig,l)',alim_star(ig,l)
     108                    PRINT*,'detr_star(ig,l)',detr_star(ig,l)
    109109                    abort_message = ''
    110                     labort_physic=.true.
     110                    labort_physic=.TRUE.
    111111                    CALL abort_physic (modname,abort_message,1)
    112112               endif
    … …  
    166166
    167167!      if (prt_level.ge.10) &
    168 !    &    call printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
     168!    &    CALL printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
    169169!    &    ptimestep,masse,entr,detr,fm,'2  ')
    170170
    … …  
    201201         do ig=1,ngrid
    202202            if (fm(ig,l+1)<0.) then
    203 !              print*,'fm1<0',l+1,lmax(ig),fm(ig,l+1)
     203!              PRINT*,'fm1<0',l+1,lmax(ig),fm(ig,l+1)
    204204                ncorecfm1=ncorecfm1+1
    205205               fm(ig,l+1)=fm(ig,l)
    … …  
    270270
    271271
    272          labort_physic=.false.
     272         labort_physic=.FALSE.
    273273         do ig=1,ngrid
    274274            if (entr(ig,l)<0.) then
    275                labort_physic=.true.
     275               labort_physic=.TRUE.
    276276               igout=ig
    277277               lout=l
    … …  
    280280
    281281         if (labort_physic) then
    282             print*,'N1 ig,l,entr',igout,lout,entr(igout,lout)
     282            PRINT*,'N1 ig,l,entr',igout,lout,entr(igout,lout)
    283283            abort_message = 'entr negatif'
    284284            CALL abort_physic (modname,abort_message,1)
    … …  
    310310         enddo
    311311
    312          labort_physic=.false.
     312         labort_physic=.FALSE.
    313313         do ig=1,ngrid
    314314            if (entr(ig,l)<0.) then
    315                labort_physic=.true.
     315               labort_physic=.TRUE.
    316316               igout=ig
    317317            endif
    … …  
    319319         if (labort_physic) then
    320320            ig=igout
    321             print*,'ig,l,lmax(ig)',ig,l,lmax(ig)
    322             print*,'entr(ig,l)',entr(ig,l)
    323             print*,'fm(ig,l)',fm(ig,l)
     321            PRINT*,'ig,l,lmax(ig)',ig,l,lmax(ig)
     322            PRINT*,'entr(ig,l)',entr(ig,l)
     323            PRINT*,'fm(ig,l)',fm(ig,l)
    324324            abort_message = 'probleme dans thermcell flux'
    325325            CALL abort_physic (modname,abort_message,1)
    … …  
    348348         enddo
    349349
    350          labort_physic=.false.
     350         labort_physic=.FALSE.
    351351         do ig=1,ngrid
    352352            if (detr(ig,l)<0.) then
    353                labort_physic=.true.
     353               labort_physic=.TRUE.
    354354               igout=ig
    355355            endif
    … …  
    357357        if (labort_physic) then
    358358               ig=igout
    359                print*,'cas 2 : ig,l,lmax(ig)',ig,l,lmax(ig)
    360                print*,'detr(ig,l)',detr(ig,l)
    361                print*,'fm(ig,l)',fm(ig,l)
     359               PRINT*,'cas 2 : ig,l,lmax(ig)',ig,l,lmax(ig)
     360               PRINT*,'detr(ig,l)',detr(ig,l)
     361               PRINT*,'fm(ig,l)',fm(ig,l)
    362362               abort_message = 'probleme dans thermcell flux'
    363363               CALL abort_physic (modname,abort_message,1)
    … …  
    398398!             lmax(ig)=l+1
    399399!             zmax(ig)=zlev(ig,lmax(ig))
    400 !             print*,'alpha>1',l+1,lmax(ig)
     400!             PRINT*,'alpha>1',l+1,lmax(ig)
    401401              ncorecalpha=ncorecalpha+1
    402402           endif
    … …  
    413413
    414414!      if (prt_level.ge.10) &
    415 !    &    call printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
     415!    &    CALL printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
    416416!    &    ptimestep,masse,entr,detr,fm,'8  ')
    417417
    … …  
    423423
    424424      if (1==1) then
    425       labort_physic=.false.
     425      labort_physic=.FALSE.
    426426      do l=1,nlay-1
    427427         do ig=1,ngrid
    … …  
    446446                         igout=ig
    447447                         lout=l
    448                          labort_physic=.true.
     448                         labort_physic=.TRUE.
    449449                      endif
    450450                      entr(ig,l+1)=entr(ig,l+1)-ddd
    … …  
    460460                         ig=igout
    461461                         l=lout
    462                          print*,'ig,l',ig,l
    463                          print*,'eee0',eee0
    464                          print*,'ddd0',ddd0
    465                          print*,'eee',eee
    466                          print*,'ddd',ddd
    467                          print*,'entr',entr(ig,l)
    468                          print*,'detr',detr(ig,l)
    469                          print*,'masse',masse(ig,l)
    470                          print*,'fomass_max',fomass_max
    471                          print*,'masse(ig,l)*fomass_max/ptimestep',masse(ig,l)*fomass_max/ptimestep
    472                          print*,'ptimestep',ptimestep
    473                          print*,'lmax(ig)',lmax(ig)
    474                          print*,'fm(ig,l+1)',fm(ig,l+1)
    475                          print*,'fm(ig,l)',fm(ig,l)
     462                         PRINT*,'ig,l',ig,l
     463                         PRINT*,'eee0',eee0
     464                         PRINT*,'ddd0',ddd0
     465                         PRINT*,'eee',eee
     466                         PRINT*,'ddd',ddd
     467                         PRINT*,'entr',entr(ig,l)
     468                         PRINT*,'detr',detr(ig,l)
     469                         PRINT*,'masse',masse(ig,l)
     470                         PRINT*,'fomass_max',fomass_max
     471                         PRINT*,'masse(ig,l)*fomass_max/ptimestep',masse(ig,l)*fomass_max/ptimestep
     472                         PRINT*,'ptimestep',ptimestep
     473                         PRINT*,'lmax(ig)',lmax(ig)
     474                         PRINT*,'fm(ig,l+1)',fm(ig,l+1)
     475                         PRINT*,'fm(ig,l)',fm(ig,l)
    476476                         abort_message = 'probleme dans thermcell_flux'
    477477                         CALL abort_physic (modname,abort_message,1)
    … …  
    494494!     if (ncorecfm1+ncorecfm2+ncorecfm3+ncorecfm4+ncorecfm5+ncorecalpha > 0 ) then
    495495
    496 !         print*,'PB thermcell : on a du coriger ',ncorecfm1,'x fm1',&
     496!         PRINT*,'PB thermcell : on a du coriger ',ncorecfm1,'x fm1',&
    497497!   &     ncorecfm2,'x fm2',ncorecfm3,'x fm3 et', &
    498498!   &     ncorecfm4,'x fm4',ncorecfm5,'x fm5 et', &
    … …  
    505505
    506506!      if (prt_level.ge.10) &
    507 !    &    call printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
     507!    &    CALL printflux(ngrid,nlay,lunout1,igout,f,lmax,lalim, &
    508508!    &    ptimestep,masse,entr,detr,fm,'fin')
    509509
Note: See TracChangeset for help on using the changeset viewer.