--- trunk/libf/phylmd/clesphys2.f90 2008/07/25 19:59:34 13 +++ trunk/phylmd/clesphys2.f 2014/03/05 14:57:53 82 @@ -1,47 +1,42 @@ module clesphys2 - ! v 1.3 2005/06/06 13:16:33 fairhead + ! From version 1.3 2005/06/06 13:16:33 implicit none - LOGICAL:: cycle_diurne= .TRUE. + LOGICAL:: cycle_diurne = .TRUE. ! Cette option permet d'éteindre le cycle diurne. ! Peut être utile pour accélérer le code. - LOGICAL:: soil_model= .TRUE. - ! Help = 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. - ! Help = Permet de mettre en route la - ! nouvelle parametrisation de l'eau liquide + LOGICAL:: soil_model = .TRUE. + ! Choix du modele de sol (Thermique ?) - LOGICAL:: ok_orodr= .TRUE. ! orodr pour l'orographie + LOGICAL:: new_oliq = .TRUE. + ! Permet de mettre en route la nouvelle parametrisation de l'eau liquide + ! Pour l'orographie: + LOGICAL:: ok_orodr = .TRUE. LOGICAL:: ok_orolf = .TRUE. - ! Desc = Orolf ou non pour l orographie - LOGICAL:: ok_limitvrai= .FALSE. - ! Help = 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 ? + LOGICAL:: ok_limitvrai = .FALSE. + ! On peut forcer le modele a lire le fichier SST de la bonne + ! annee. - INTEGER:: nbapp_rad= 12 - ! (nombre d'appels des routines de rayonnements par jour) + INTEGER:: nbapp_rad = 12 + ! 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 + namelist /clesphys2_nml/cycle_diurne, soil_model, new_oliq, & ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con @@ -49,16 +44,12 @@ print *, "Enter namelist 'clesphys2_nml'." read(unit=*, nml=clesphys2_nml) - write(unit=*, nml=clesphys2_nml) + write(unit_nml, 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 + if (iflag_con <= 1) then + print *, "bad value for iflag_con" + stop 1 + end if end subroutine read_clesphys2