--- trunk/libf/phylmd/clesphys2.f90 2011/09/20 09:14:34 51 +++ trunk/phylmd/clesphys2.f 2014/07/02 18:39:15 99 @@ -10,35 +10,34 @@ LOGICAL:: soil_model = .TRUE. ! Choix du modele de sol (Thermique ?) - ! Option qui pourait un string afin de pouvoir - ! plus de choix ! Ou meme une liste d'options LOGICAL:: new_oliq = .TRUE. ! Permet de mettre en route la nouvelle parametrisation de l'eau liquide - LOGICAL:: ok_orodr = .TRUE. ! pour l'orographie - LOGICAL:: ok_orolf = .TRUE. ! pour l'orographie + ! Pour l'orographie: + LOGICAL:: ok_orodr = .TRUE. + LOGICAL:: ok_orolf = .TRUE. LOGICAL:: ok_limitvrai = .FALSE. ! On peut forcer le modele a lire le fichier SST de la bonne - ! annee. C'est une tres bonne idee, pourquoi ne pas mettre toujours - ! a y ? + ! annee. INTEGER:: nbapp_rad = 12 - ! (nombre d'appels des routines de rayonnements par jour) + ! nombre d'appels des routines de rayonnements par jour INTEGER:: iflag_con = 2 - ! Flag pour la convection : - ! 1 LMD, - ! 2 Tiedtke, - ! 3 CCM(NCAR) - ! 3 KE - ! 4 KE vect + ! Convection scheme: + ! 2 Tiedtke + ! 3 Emanuel + ! 4 Emanuel vect contains subroutine read_clesphys2 + use unit_nml_m, only: unit_nml + use nr_util, only: assert + namelist /clesphys2_nml/cycle_diurne, soil_model, new_oliq, & ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con @@ -46,16 +45,8 @@ print *, "Enter namelist 'clesphys2_nml'." read(unit=*, nml=clesphys2_nml) - write(unit=*, nml=clesphys2_nml) - - select case (iflag_con) - case (1) - PRINT *, 'Schéma convection LMD' - case (2) - PRINT *, 'Schéma convection Tiedtke' - case (3) - PRINT *, 'Schéma convection CCM' - END select + write(unit_nml, nml=clesphys2_nml) + call assert(iflag_con >= 2, "read_clesphys2 iflag_con") end subroutine read_clesphys2